Difference between revisions of "Atmosphere"

Revision as of 21:22, 2 September 2013

The pressures for all atmospheres

Planets				Moons
	Eve		Kerbin		Laythe
	Duna		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.

Drag

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

${\displaystyle F_{D}=0.5\,\rho \,v^{2}\,d\,A}$

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

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 m²/kg, so:

${\displaystyle 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 (p₀), and scale height (H):

${\displaystyle p=p_{0}\cdot e^{\frac {-altitude}{H}}}$

${\displaystyle \rho =1.2230948554874\cdot p}$

where p here is in units atm, and ρ in kg/m³. The 1.2230948554874 kg/m³ value for density at 1 atmosphere (sea level on Kerbin) is given by FlightGlobals.getAtmDensity(1.0).

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:

${\displaystyle {\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 (F_G) is:

${\displaystyle 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 F_G equal to F_D:

${\displaystyle {\begin{aligned}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{aligned}}}$

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

${\displaystyle 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:

${\displaystyle v_{T}={\sqrt {\frac {250{\frac {\text{kg}}{{\text{m}}^{2}}}\cdot GM}{r^{2}\,\rho \cdot 35.07}}}}$

${\displaystyle \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}}}}}$

${\displaystyle 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

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:

Atmospheric height

The atmospheric height depend on the scale height of the celestial body and is where 0.000001^th (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).

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

${\displaystyle 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:

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

Notes

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