This document contains equations, constants, and other goodies for world-building calculations. They're useful in doing any kind of hard science story, novel, or campaign. Many of these are lifted from the non-copyrighted version of the sci.space FAQ. Others are taken from Stephen Gillett's excellent book, World Building. See the biblography for more details.
Some of these equations and technologies involved are discussed in a sister page, Interplanetary Travel: An Intro (with Equations), though that's more about rocketry. Another good educational link is NASA
Where d is distance, v is velocity, and t is time.
For constant acceleration:
Distance | d = d0 + vt + 0.5at2 |
Velocity | v = v0 + at |
Velocity squared | v2 = 2ad |
Acceleration on a cylinder of radius r and rotation period t: | a = 4 p2r/t2 |
Time to travel distance d at acceleration d, given constant acceleration half-way and constant deceleration half-way | t = 2*(da)1/2 |
Surface gravity | g = GM/r2 |
Surface gravity of planet p as fraction of earth's | g = rp/re * densityp/densitye |
Escape velocity | vesc = 20.5 * vc = (2GM/r)1/2 |
Orbital velocity | vorbital = (GM/a)1/2 |
Tides (in earth units)
| T = M/R3 |
Absolute magnitude | M |
Apparent magnitude | m |
Distance in parsecs | p |
Luminosity (in solar units) | L |
intensity (solar constant = 1) | I |
Distance of planet (in AU) | R |
Diameter of star (sol = 1) | D |
Effective temperature (sol) | T (5770 K) |
Effective temperature (star) | t |
Size in degrees | S |
Absolute magnitude | M = m+5 - 5log p |
Luminosity | L = 2.52(4.85 - M) |
Apparent brightness | I = L/R2 |
Stellar diameter | D = L(T2/t2) |
Size in sky: (for sizes ~20 degrees) | S = 57.3D/R |
If you're using stars that are somewhat more extreme, you might want to calculate the bolometric magnitude instead. (Bolometric magnitude is the total amount of radiation put out by the star, not just visible radiation.) Add the correction values from this table (stolen from Gilett) to the magnitude of the star. Gillet says, "For a more complete table showing additional classes, refer to Kaler 1997, p. 263." I'll try and dig up the Kaler reference.
You might also go see the On Creating an Earthlike Planet, which is excellent (and I'm not saying that because he links back here; I only just discovered his link).
Class | Main Sequence | Giants | Supergiants |
---|---|---|---|
O3 | -4.3 | -4.2 | -4.0 |
B0 | -3.0 | -2.9 | -2.7 |
A0 | -0.15 | -0.24 | -0.3 |
F0 | -0.01 | 0.01 | 0.14 |
G0 | -0.10 | 00.13 | -0.1 |
K0 | -0.24 | -0.42 | -0.38 |
M0 | -1.21 | -1.28 | -1.3 |
M8 | -4.0 |
If you really need to calculate it, there's an empirical formula and a calculator at http://www.go.ednet.ns.ca/~larry/astro/HR_diag.html.
If you’re interested in learning about orbital mechanics, try the Orbital Mechanics website, http://www.braeunig.us/space/orbmech.htm.
T = 2p/(GM/a3)0.5
where
G | Gravitational constant |
M | Mass of both bodies |
r | Radius of orbit |
a | Semimajor axis of orbit |
Normally used for LEO to GEO calculations, this is T.N. Edelbaum's equation. Unless there are simplifications I’m not aware of, it should be valid for differences between any two circular orbits around the same primary:
Delta-Vee = (V12 - 2 V1 V2 cos (pi/2 alpha ) + V22)1/2
Where:
V1 | circular velocity initial orbit |
V2 | circular velocity final orbit |
alpha | plane change in degrees. |
Roche's Limit | L = 2.44 r (densityp/densitys)1/3 |
where
Density of planet | densityp |
Density of satellite | densitys |
Radius | r |
Orbital velocities for orbits at a distance r:
Semimajor axis | a |
G(m1 + m2) | m |
Circular orbit | v = [m/r]1/2 |
Elliptical orbit | v = [m((2/r) - (1/a))]1/2 |
Parabolic orbit | v = [m(2/r)]1/2 |
Hyperbolic orbit | v = [m( (2/r)+(1/a))]1/2 |
Energy of object in orbit | E = -Gm1m2/2a |
Eccentricities of orbits depending on orbit type, with semimajor axis a and semiminor axis b:
Circular orbit | e = 0 |
Elliptical orbit | e < 1 |
Parabolic | e = 1 |
Hyperbolic orbit | e > 1 |
Point of periapsis | Rp = a(1-e) |
Point of apoapsis | Ra = a(1+e) |
Note: | 2a = Rp + Ra |
Eccentricity of orbit | e = Rp * Vp2 / GM |
Eccentricity of orbit | e = (a2 + b2)1/2/a |
Period of orbit | P2 =4p2/ma3 |
P = 2p/[ma3]1/2 |
Energy of an object of mass m in an orbit around the sun (mass M) with semimajor axis a:
Orbital energy | E = -G*M*m/(2a) |
Insolation of a planet determines approximately how much light it gets, and (in solar units) depends on the luminosity of the star and its distance. Dole suggests the acceptable values are between 0.65 and 1.5; see the "fudged temperature" for a more recent measurement.
Insolation (relative) | I = L/D2 |
Luminosity of star | L |
Distance from star | D |
A Hohman transfer orbit is the minimum energy orbit to get from planet A to planet B, assuming they have circular Keplerian orbits. The orbit is circular, with a tangent at the perihelion of one planet and another tangent at the aphelion of the other.
Semimajor axis of planet 1 | R1 |
Semimajor axis of planet 2 | R2 |
Semimajor axis of the transfer orbit | a=(R1+R2)/2 |
Once you have the semimajor axis, you know transfer time: it's half the orbital period for a circular Keplerian orbit of that radius (use equation above).
To calculate required DV, you need to know the orbital velocity for your transfer orbit at the points where it's tangential to the orbits of the departure and destination planets:
Orbital velocity V = ( (2GM * [1/r - 1/2a])1/2
Ignoring for now the problems of calculating the angle that the destination planet needs to subtend and calculating the launch date; sample calculations for Earth to Mars can be found at:
http://www.marsacademy.com/text/angplan.htm
http://www.marsacademy.com/text/ladate.htm
There's a second kind of easily-calculated, efficient orbit, one that assumes a constant low acceleration (the sort you'd expect from an ion drive or a solar sail).
The acceleration must be very much lower than R/P2, where R is the distance from the sun and P is the period of the outermost planet. (Note that this is extremelylow; the value of R/P2 for Earth is 0.015 m/s2; for Mars, it is 0.0065 m/s2, or less than 7 ten-thousandths of a G.)
An acceptable approximation of the travel time is:
2pR1/(aP1) * (1 - R1/R2)0.5
where:
R1 | Distance from sun of inner planet |
R2 | Distance from sun of outer planet |
P1 | Period of the inner planet |
a | Aceleration of spacecraft |
Be consistent with your units! If a is in m/s2, then P must be in seconds.
You can get a value good enough for story or RPG purposes by doubling t=(2d/a)0.5, where d is half the distance to the other planet. For example, say that Mars to Earth is 7.893E10 meters, its closest approach (2.279E11 - 1.496E11 meters). You have a solar sail that gives you 0.001 G acceleration, or 0.01 m/s2. The time to accelerate half-way there is:
(7.83E10/0.01)0.5 = (7.83E12)0.5 = 2.8E6 seconds
Or a little over 32 days. Assume the same time tod ecelerate for a total Earth-to-Mars time of about 65 days.
A note from the website http://dutlsisa.lr.tudelft.nl/Propulsion/Data/V_increment_requirements.htm says:
"Transfer or trip time for constant thrust spiral is calculated by dividing total propellant mass by mass flow. Total propellant mass is calculated using the rocket equation also known as Tsiolkowky's equation. In case of negligible propellant mass (constant acceleration), transfer time can be calculated by dividing the velocity change by the acceleration."
Of course, it doesn't look to me that you get to stop at the other end, and you're not necessarily going at the same velocity as your target...but good enough for back-of-the-envelope calculations.
Temperature of a blackbody:
Albedo | A |
Incident light (sun=1) | I |
Temp in degrees Kelvin | T |
Temperature | T = 374(1-A)I1/4 |
To allow for greenhouse gases, Gillett suggests a fudge factor of about 1.1 for habitable planets:
Fudged temperature | T = 374 * 1.1(1-A)I1/4 |
Intensity of blackbody per unit area:
Stephann-Boltzmann constant | s |
Temperature, degrees K | T |
Intensity | I = sT4 |
GURPS Space uses a variant on this "law" (discovered by Titius, popularized by Bode) for placing planets.
Gillett says that current thinking is this is an example of tidal separations in the protocloud; it holds to lesser extents for moon systems as well, but with different parameters.
The classical formula for our solar system, where rn is the orbital distance for planet n:
rn = (0.3 * 2n) + 0.4 AU
A more general form, suitable for moons around planets, for planet n:
Pn = P0An
Where:
Pn | Period of orbit of nth planet (traditionally in days) |
P0 | Period of primary's rotation |
A | Semimajor axis of the orbit |
Warning: I no longer remember where I got the general form, and frankly, looking at it I no longer understand it. So use with caution (or point me to the derivation, so I can do the same for others).
Where dv is the change in velocity, Isp is the specific impulse of the engine, ve is the exhaust velocity, x is the reaction mass, m1 is the rocket mass excluding reaction mass, g is acceleration due to gravity on earth:
Exhaust velocity | ve = g Isp |
Change in velocity | DV = ve * ln((m1 +x)/m1) |
Or: | |
Ratio of masses | (m1+x)/m1) = e(d/v) |
Note that (m1+x)/m1) is the ratio of the initial mass to the final mass.
The exponent d/v is change in velocity over exhaust velocity.
For a staged rocket where each stage has the same ratio R of initial to final mass and with n stages:
Final delta-vee | D |
V = n [veln(R)] |
You may notice that's the same as the single stage orbit multiplied by n. Essentially, two stages give you twice the final velocity of a single stage rocket with the same mass ratio, and so on.
For constant acceleration:
Time (unaccel.) | tu = (c/a) * sinh(at/c) |
Distance | d = (c2/a) * (cosh(at/c)-1) |
Velocity | v = c * tanh(at/c) |
For a black hole of mass M:
Schwartzchild radius | 2GM/c2 |
Some useful constants. Since it's sometimes easier to work things out in solar or terran equivalents, some physical data for our solar system is also included.
For game purposes, one or two significant digits is all you need, but I've gone to four here.
G (gravitational constant) | 6.673E-11 Nm2/kg2 |
c (speed of light in vacuum) | 2.998E8 m/s |
Luminosity of sun | 3.827E26 W |
Solar constant (intensity@1 AU) | 1370 W/m2 |
Planck's constant h | 6.6262E-34 J-s |
"h bar" h/(2p) | 1.055E34 J-s |
Boltzmann's Constant k | 1.381E-23 J/K |
Stephann-Boltzman Constant s | 5.670E-8 W/m2/K |
Earth gravity acceleration | 9.80665 m/s2 |
One light year (meters) | 9.461E15 m |
One light year (AU) | 2.063E5 AU |
One parsec (light years) | 3.262 ly |
One parsec (meters) | 3.086E16 m |
Mean earth-moon distance | 3.844E8 m |
Mean earth-sun distance (1 AU) | 1.496E11 m |
Mean radius of earth | 1.371E6 m |
Equatorial radius of earth | 6.378E6 m |
Mass of Sun | 1.989E30 kg |
Mass of earth | 5.974E24 kg |
Mass of moon | 7.348E22 kg |
Average density of Earth | 5.5 g/cm3 |
VDS (from BTRC) has the power from a solar sail constant per square kilometer of sail (10 w at Earth orbit). Fiddling with the acceleration equation gives us these two versions:
a2 = 2P/(5M) |
M=2P/(5a2) |
Where a is the acceleration in meters/s2, M is the mass of the vehicle in kilogarms, and P is the power in watts. If you calculate the acceleration for a given vehicle at earth orbit (10 watts), the acceleration at other orbits is proportional to the distance in AU. (Calculating for another star is a different matter.) Sails for different TLs have the following "characteristic acceleration" (acceleration with no payload):
TL | Power | Mass/km2 | Characteristic acceleration |
---|---|---|---|
11 | 10w | 413 | 0.089 |
12 | 10w | 347 | 0.107 |
13 | 10w | 296 | 0.116 |
14 | 10w | 255 | 0.125 |
15 | 10w | 222 | 0.134 |
Equations and data were taken from the following references:
Get your own Free Home Page