45381
Comment:
|
45967
|
Deletions are marked like this. | Additions are marked like this. |
Line 230: | Line 230: |
$ A ~=~ { \Large { { ( 1 - {e_L}^2 ) ~ a_L } \over { ( 1 - {e_c}^2 ) ~ a_c } } } $ So that $ A ( 1 + e_c \cos( \theta ) ~=~ 1 + e_L \cos( \theta + \theta_L ) $ Rearranging: $ ( A - 1 ) /e_L + ( A e_c / e_L ) \cos( \theta ) ~=~ \cos( \theta + \theta_L ) $ |
$ K ~=~ { \Large { { ( 1 - {e_L}^2 ) ~ a_L } \over { ( 1 - {e_c}^2 ) ~ a_c } } } $ So that $ K ( 1 + e_c \cos( \theta ) ~=~ 1 + e_L \cos( \theta + \theta_L ) $ Rearranging: $ ( K - 1 ) /e_L + ( K e_c / e_L ) \cos( \theta ) ~=~ \cos( \theta + \theta_L ) $ |
Line 238: | Line 238: |
$ ( A - 1 ) / e_L + ( A e_c / e_L ) \cos( \theta ) ~=~ \cos( \theta ) \cos( \theta_L ) - \sin( \theta ) \sin( \theta_L ) $ We know all the terms except $ \cos( \theta ) $ and $ \sin( \theta ) $; how do we solve for those? Let's solve the equation for $ \cos( theta ) $ in terms of the constants and $ \sin ( theta ) $: $ $ |
$ ( K - 1 ) / e_L + ( K e_c / e_L ) \cos( \theta ) ~=~ \cos( \theta ) \cos( \theta_L ) - \sin( \theta ) \sin( \theta_L ) $ We know all the terms except $ \cos( \theta ) $ and $ \sin( \theta ) $; how do we solve for those? Let's solve the equation for $ \cos( \theta ) $ in terms of the constants and $ \sin ( \theta ) $: $ \cos( \theta_L ) \cos( \theta ) - ( K e_c / e_L ) \cos( \theta ) ~=~ ( K - 1 ) / e_L + \sin( \theta_L ) \sin( \theta ) $ $~~~~ ( \cos( \theta_L ) - ( K e_c / e_L ) ) \cos( \theta ) ~=~ ( K - 1 ) / e_L + \sin( \theta_L ) \sin( \theta ) $ $~~~~ \cos( \theta ) ~=~ $ $ \Large { { ( K - 1 ) / e_L + \sin( \theta_L ) \sin( \theta ) } \over { \cos( \theta_L ) - ( K e_c / e_L ) } } $ Define two more intermediate variables: $ M = { \Large { { K - 1 } \over { e_L \cos( \theta_L ) - K e_c } } } $ and $ N = { \Large { \sin( \theta_L ) \over { e_L \cos( \theta_L ) - K e_c } } } $ |
Construction1
This will be merged with ConstructionOrbits real soon now.
This analysis ignores second-order effects like lunar/solar/Jupiter tides, apsidal precession from Earth oblateness, light pressure, and outgassing, which will add seconds and many meters of perturbations to the precision calculations needed to deliver a vehicle to a rendezvous at an exact position and time. With years of experience and calibration, the launch loop will evolve into a precision instrument for millimeter-accurate delivery of vehicles to destinations.
However, the following calculations will be better than 0.1% accurate, useful for estimating delivery rates and infrastructure requirements. It is better to be roughly right than exactly wrong.
Construction orbit, Act 1
Begin with a construction orbit, with a period of one sidereal day. There are 366.2422 sidereal days per 365.2422 solar days of 86400 seconds, so the construction orbit period is ~~~~ T_c = 86164.0905 seconds
The angular frequency of the construction orbit is ~~~~ \omega_c = 2 \pi / T_c = 7.2921158e-5 radians/second
The Earth's standard gravitational constant is μ = 398600.4418 km³/s². The semimajor axis of a one sidereal day orbit is ~~~ a_c = ( \mu ( T_c / 2\pi )^2 )^{1/3} = 42164.1696 km
The construction orbit perigee should be above the LEO space debris, collision, and litigation belt. Arbitrarily choose a perigee radius 2000 km above Earth equatorial radius: ~~ r_{cp} = 2000 + 6378 km = 8378 km
The construction orbit apogee is twice the semimajor axis on the opposite side of the Earth from the (negative) perigee, so the apogee radius is ~~~ r_{ca} = 2 a_c - r_{cp} = 2 * 42164.1696 - 8378 = 75950.3392 km
Deviating from "standard" treatments, we will choose our reference angle, time, and other orbit parameters in relation to apogee, not the usual perigee. The eccentricity of this orbit is negative: ~~~ e_c = r_{cp} / a_c - 1 = 8378.0 / 42164.1696 - 1.0 = -0.8013005 (no units)
The construction orbit perigee occurs one half a period before time zero: T_{cp} = 0.5 T_c = -43082.0452 seconds
The construction orbit perigee velocity is:
v_{cp} = \sqrt{ ( \mu / a_c ) ( r_{ca} / r_{cp} ) } = \sqrt{ ( { \small 398600.4418 } / { \small 42164.1696 } ) \times ( { \small 75950.3392 } / \small { 8378 } ) ~ ~ } = 9.25746 km/s
The construction orbit apogee velocity is:
v_{ca} = \sqrt{ ( \mu / a_c ) ( r_{cp} / r_{ca} ) } = \sqrt{ ( { \small 398600.4418 } / { \small 42164.1696 } ) \times ( { \small 8378} / \small { 75950.3392 } ) ~ ~ } = 1.02118 km/s
One sidereal day construction orbit |
|||||||||||
Apogee |
Perigee |
Semimajor |
eccentricity |
Apogee |
Perigee |
period |
apogee |
perigee |
|||
km |
km |
Axis, km |
|
m/s |
m/s |
sec |
time, sec |
time, sec |
|||
75950.34 |
8378.00 |
42164.17 |
-0.8013005 |
1021.28 |
9257.46 |
86164.09 |
0.0 |
-43082.05 |
There can be thousands of different one day construction orbits with the same apogee, perigee, and period, but arriving at apogee at different times throughout the sidereal day. The constellation will resemble the petals of a flower; while the orbital tracks intersect, two different stations associated with the same launch loop will not pass through the same intersection simultaneously.
Unless the entire construction station is well shielded, the crew must be confined to a shielded core surrounded by lots of hydrogen (food, waste, water, fuel) that can attenuate the dose. If the orbit passes through the radiation belt only once per day, this can be the rigidly scheduled sleeping period, every 23 hours and 56 minutes.
Astronauts are superbly competent people, so they can get their day's work done 4 minutes faster than the rest of us.
If the construction orbit perigee is "sidereal midnight", apogee is "sidereal noon". As we will see, "noon" is when passengers and supplies arrive from Earth, and when passengers and products return for reentry. Launch and entry orbits will have smaller perigees, a smaller semimajor axis, and thus somewhat shorter orbit times than the construction orbit.
Multiple Day Construction Orbits
The period of the construction orbit need not be one sidereal day, but it should be an integer ratio fraction of a sidereal day.
For example, a 5/3 sidereal day orbit will have a larger semimajor axis than the one day orbit, proportional to the 2/3 power of the period, or 59271.063 km. The perigee might still be 8378 km, yielding a 110164.127 km apogee. Since the launch loop rotates around the Earth once per sidereal day, the apogee of this construction orbit is accessable only once every three orbits, or every five sidereal days. Such orbits are interesting, but not as practical as simple integer ratios. A 5 sidereal day orbit will actually be more useful - one third of the passes through the van Allen belt, and a higher perigee velocity for the start (and end) of interplanetary missions.
A more useful example is a 1/2 sidereal day orbit, with a semimajor axis of 26561.762 km. If perigee is 8378 km as before, apogee is 44745.524 km. This orbit can intersect with loop launches twice per day. However, this accentuates a problem shared by all orbits with "low" perigees; they pass through the van Allen belt during descent and ascent more often and more slowly than higher orbits, so more crew time spent in the shielded core. Indeed, the peak of the outer belt is at 5 Re, about 32000 km, Such orbits are more practical for crew-less automated stations; if the crew cannot work outside shielding, then telepresence is a more practical way to work there.
High integer multiple orbits: 2, 3, 4 etc. sidereal days, may be useful, but orbits that with perigees near lunar radius may be (WAG) more sensitive to lunar gravitational perturbations. The Moon's semimajor axis is 384399 km, though the perigee and apogee drift due to solar tidal forces - a similar drift will perturb construction orbits and launches to them.
A perigee of 8378 km and an apogee of 384400 km is 196389 km, with an orbital period of 866137 seconds, or 10.05 sidereal days.
So, the most practical construction orbits are likely to be 1, 2, 3 ... 6 sidereal days. For the rest of this discussion, we will focus on launches to the 1 sidereal day construction orbit; this may offer the best compromise between rapid assembly and a higher v_{\infty} .
We will continue the construction orbit description later, after considering how to get there from a launch loop.
Prime orbit, Act 1
From here on out, we will round distance to the nearest kilometer and the time to the nearest second.
The prime orbit is an unmodified launch orbit, with the same apogee time ( 0.0 seconds ) and the same apogee radius as the one sidereal day construction orbit: r_{pa} = 75950 km.
The launch orbit perigee is the altitude of the launch loop (assumed to be 80 km) added to the radius of the Earth: r_{pp} = 6378 + 80 km = 6458 km. The prime orbit semimajor axis is: a_p = 0.5 * ( 6458 + 75950 ) = 42104 km The period of the prime orbit is: T_p = 2 \pi \sqrt{ a_p / \mu } = 85980 seconds.
The trajectory time is half the full orbit, or 42990 seconds, or 11 hours 56 minutes and 30 seconds. Prime orbit perigee time is before zero time, or t_{pp} = -42990 seconds.
That is 92 seconds after the construction orbit passed overhead. The prime orbit perigee velocity is 10551 m/s. 471 m/s is provided by earth rotation velocity, and an additional 30 m/s WAG) is needed for air friction during exit. The estimate loop-related launch velocity is 10110 m/s, or 337 seconds of acceleration at 30 m/s . The acceleration run begins 337-92 = 245 seconds before construction orbit perigee, and 1700 km to the west of the release point.
The construction orbit station is moving at 9.257 km/s, so it will be aproximately 2200 km west of the release point, and 2000 km overhead, or 500 km west of launch loop west station, at an elevation of atan( 500/2000 ) or 14 degrees from zenith. There will be plenty of warning time to abort the launch from the loop, though that can lead to logistic complications for fast assembly.
Prime Launch
The launch loop launch track can be at altitudes between 70 km and 120 km. The drag and heating at 120 km is very low, but that means space debris deorbits more slowly, greatly increasing the debris flux and impact hazard. At 70 km loop altitude, the drag is 4.7x higher and the temperature (in Kelvins) is 11% higher, so the vehicle drag and heating is 5x higher than the slightly higher 80 km, making the vehicle nose shroud heavier and more expensive. We will choose 80 km for this exercise.
The earth rotates one sidereal day in 86141.0905 seconds. At 6458 km radius, that is a rotation velocity of ( 2 π 6458000 / 86141.0905 ) = 471.05 m/s . The launch velocity is the apogee velocity minus the rotation velocity. More precisely, we must add a little bit of extra launch velocity to counteract drag leaving the atmosphere; the amount of drag loss will depend on vehicle aerodynamics and atmospheric conditions, but will probably be less than 30 m/s. Production launches will use realtime atmospheric measurements and exact vehicle measurements to model this effect, resulting in exact trajectories trimmed to millimeters per second.
Prime launch delivers vehicles to the apogee of the construction orbit from the launch site at 6458 km radius (80 km altitude). Before a more versatile capture system is assembled, this is the easiest way to get initial materials to the construction orbit station, without a plane change, once per sidereal day.
The arrival velocity is the prime orbit apogee velocity minus the 1021.28 m/s construction station apogee velocity. That velocity will be faster, and some onboard vehicle thrust will be needed to bring incoming vehicles up to speed. Later, we will discuss cable capture, a way to minimize vehicle cost and complexity by using construction station momentum and makeup thrust.
Orbital characteristics referenced from apogee:
Prime Launch Orbit |
|||||||||||
Apogee |
Perigee |
Semimajor |
eccentricity |
Apogee |
Perigee |
period |
arrival |
launch |
launch |
arrival |
|
km |
km |
Axis, km |
|
m/s |
m/s |
sec |
time, sec |
time, sec |
m/s |
m/s |
|
75950 |
6458 |
42104 |
-0.84326 |
897.2 |
10551.7 |
85979.72 |
0.0 |
-42989.86 |
10081 |
142.1 |
A 142 m/s apogee boost is much smaller than the injection velocity for a circular orbit at 75950 km radius; indeed, it is smaller than the injection velocity into any circular orbit above 275 km altitude. This is why construction orbits are the best destinations from the launch loop.
Reentry after Capture/Rendezvous failure
If the vehicle fails to rendezvous with the station, it cannot quickly resynchronize for a second pass. The prime orbit period is shorter than the one day sidereal orbit, so the vehicle will arrive after another orbit 161.37 seconds ( 86141.09 - 85979.72 ) early. It may be best to re-enter (assuming a heat shield) and re-launch. If a low-cost cargo vehicle lacks a heat shield, let it re-enter and burn up. A 40 km altitude ( 6418 km radius) is assumed sufficient to lose enough velocity to descend, so a launch orbit abort reentry orbit will look like so:
75950 km Apogee Prime Reentry Orbit |
|||||||||||
Apogee |
Perigee |
Semimajor |
eccentricity |
Apogee |
Perigee |
period |
apogee |
entry |
slowdown |
entry |
|
km |
km |
Axis, km |
|
m/s |
m/s |
sec |
time, sec |
time, sec |
m/s |
m/s |
|
75950 |
6418 |
41184 |
-0.84416 |
0.90435 |
10.702 |
83177.6 |
0.0 |
41588.8 |
2.59 |
10234.1 |
This delta V can be accomplished with laser ablation. In addition, the vehicle should have a rocket capable of 3 m/s retrograde thrust.
A much larger retrograde thrust at can re-enter from the construction orbit, starting at apogee:
75950 km Apogee Construction Station Reentry Orbit |
|||||||||||
Apogee |
Perigee |
Semimajor |
eccentricity |
Apogee |
Perigee |
period |
apogee |
entry |
slowdown |
entry |
|
km |
km |
Axis, km |
|
m/s |
m/s |
sec |
time, sec |
time, sec |
m/s |
m/s |
|
75950 |
6418 |
41184 |
-0.84416 |
0.90435 |
10.702 |
83177.6 |
0.0 |
41588.8 |
116.83 |
10234.1 |
Prime Capture at Docking Station
Approaches to ISS involve a complicated set of manuevers that take up to 8 hours. ISS is a bag of air wrapped in aluminum foil, easily damaged by a collision. A capture station will be larger, heavier and far more robust, and approach speeds from launch orbits will be slower. An amusing alternative is an armored lower docking station suspended on stout cables below the main station, plus deployable capture lassos that grab the tailhooks of approaching vehicles. This is an incredible kludge, but it means that system propulsion can be a fixed rocket on the docking station, shared among hundreds of vehicles a year (some of them propellant tankers).
A single construction station may deploy dozens of docking stations, capturing many incoming vehicles sequentially, perhaps 100 tonnes per day. This is no job for pilots; computers will measure, simulate, and adjust positions and velocities to millimeter accuracy. There are no waves or wind in orbit, and computers and sensors are very fast. In zero gee and high vacuum, with laser interferometry, real-time computer simulation, and high precision actuators, we will learn to capture incoming vehicles faster and more accurately than a tailhook landing on a aircraft carrier.
An incoming vehicle with 200 m/s retrograde velocity (the station moves faster than a vehicle at apogee) moves only 200 micrometers in a microsecond, which is ample time to make an interferometric measurement, update a simulation model, and apply power to an actuator. With limited power, mechanical movements will be slow, but the system has performed nearly identical operations many times before; the control system will be correcting unexpected errors, not discovering de novo what to do. At the speed of light in optical fiber, a sensor reading can travel 200 meters in a microsecond. Compared to "electronics world", "mechanical rendezvous world" is very slow and predictable.
5000 kg vehicles launch from the loop with an acceleration of three+ gees, 30 m/s²; if vehicle contents remained well secured during the 12 hour zero gee trajectory, a 10 m/s² deceleration along the same vehicle axis is tolerable. That is 20 seconds of deceleration, and a 2000 meter deceleration run. One possibility for hookup is to launch the tethered deceleration probe slowly along the trajectory of the incoming vehicle direction, stop it perhaps 1200 meters away, then pulling it towards the docking station at 10 gees, accelerating to 200 m/s inwards in 2 seconds, exactly matching the incoming vehicle's velocity vector and position 1000 meters out.
This allows 5 seconds to manuever and latch onto the vehicle, while the deceleration probe tether is reeled in at the station. The vehicle whips by at 200 m/s, and the pull motor reverses in a couple of seconds, becoming a generator. The cable pulls on the generator; If the vehicle and cable mass 6000 kg, and the combination decelerates at 10 m/s². The force transmitted to the station will be 6000 kg × 10 m/s² = 60 kN. This decelerates the vehicle and cable to zero velocity 20 seconds, and adds a net 1000 kN-s momentum impulse to the station as vehicle, cable, and station match velocities.
The peak generated power is 6000×10×200 = 12 MW, dropping to zero as the vehicle slows to a stop.
Restoring Docking Station Momentum
This power might help expel propellant at 1600 m/s retrograde, which will drop the propellant exhaust plume down to the earth perhaps 11 hours later (and before the loop and its next launches rotate into the tenuous plume). This initially generates 14.4 kN of thrust, reducing the propellant feed as power diminishes, a total impulse of 144 kN-s.
We must generate 856 kN-s more thrust and momentum by other means during vehicle capture, requiring a total energy of perhaps 1.8 GJ. At 100 Wh/kg, that is 5000 kg of fast-discharge lithium-ion battery. Total propellant used is 1000 kN-s / 1600 m/s or 625 kg of inert propellant; perhaps water, 12.5% of incoming vehicle mass. Perhaps 1 in 6 of the incoming vehicles are water tankers.
2000 km/s exhaust molecules would be launched into a retrograde Earth orbit, and plague LEO satellites and subsequent loop launches. A plume launched retrograde with significant additional earthwards velocity will have the same angular momentum and will probably still miss the atmosphere, merely doing so more quickly.
A much faster retrograde exhaust plume velocity, such that the slowest molecules are moving faster than 4.3 km/s relative to the station, will launch the propellant entirely away from the Earth faster than escape velocity. That may require a 5000 m/s or better electric thruster and perhaps 6 GJ of stored battery energy and 17000 kg of battery. It is cheaper to launch more water with the launch loop.
Docking Station to Construction Station
Multiple docking stations may be attached by transfer elevator tethers to the main Construction Station, allowing multiple captures per sidereal day. One way to arrange the tethers is as a rotating "pinwheel" around the construction station, connected by fairly long tethers, helping protect the construction station from errant incoming vehicles. The rotation period may be fairly slow, perhaps a few revolutions per sidereal day, but fast enough to keep them taut as the capture station travels around its highly elliptical orbit. Without rotation, the docking station cables will slacken, making transfer elevator operations more difficult.
The construction station itself may include a one gee centrifugal wheel habitat, encased in a very slowly rotating shielding sheath that attaches to the docking stations; axis oriented north-south. A 50 meter radius, 0.9 gee habitat might turn in 15 seconds; 4 RPM, or 5744 turns per sidereal day. The outer shell might turn every 900 seconds, about 96 times per sidereal day. The centrifugal gravity at the end of 300 meter radial tethers would be 0.015 m/s², or 1.5 milli-gees, and a rotation velocity of 2 meters. This is enough for orientation, but not enough to fling careless workers far from construction station. This rotates the docking stations through a 900 second group of arriving vehicles, discussed in the next section. A 1 m/s² freight elevator can lift captured vehicles to the construction station hub in one minute, including low-jerk starts and stops.
Early and Late Launch
Either Before or After may not work, finish analysis soon.
The prime orbit is easiest to compute, and is coplanar with the construction orbit, but there are other trajectories from the launch loop to the construction station, launching earlier or later from prime launch time. A fully-powered launch loop can launch a 5 tonne vehicle every 45 seconds; over a +/- 900 second window, that is 40 additional vehicles and 200 extra tonnes to the construction station every sidereal day.
A launch loop oriented along the 8 degree south latitude line launches into a plane defined by that velocity vector and the center of the Earth. The launch plane rotates, and orbits launched into one plane will need some north/south delta V to transition to a different plane; for example from an early or late launch into the prime launch plane that also defines the construction orbit.
Over 24 sidereal hours (86164 seconds), the normal vector rotates around the pole with a circle with a radius of sin(8°) on the unit sphere. The angular spacing on that circle from the "prime" launch orbit is the plane change required from a launch at an earlier or later time.
For example, for a launch 900 seconds later, the Earth has turned 360*900/86164 = 3.76°. That is an angle of sin(8°) × sin(3.76°) or 9.13 milliradians around that circle. The plane change velocity near apogee is 2 x sin( 0.5 x 9.13 mrad ) x Vapogee. Vapogee is approximately 900 m/s for a one day construction launch orbit, so the plane change velocity would be approximately 8.2 m/s northward. Similarly, for a launch only 300 seconds after the prime launch orbit, the angle is 3.04 mrad and the plane change velocity is 2.7 m/s northward.
The apogee of an early launch occurs east of the construction orbit perigee, and so does the plane change. The plane change for and apogee of an early launch occurs west of the construction orbit perigee. Assuming a north-south thrust for plane change, and no radial or in-plane thrust, we can remap the problem onto one plane, and compute the intersection of the launch orbit and the construction orbit.
We can change the apogee of the launch orbit by adjusting the loop exit velocity. A slightly higher apogee and velocity requires exiting the loop slightly further east; for now we will ignore this complication and assume all launches occur at the same longitude on the Earth as the prime launch.
The Earth rotates 360 degrees in a 86164.09 second sidereal day, or 3.7603 degrees in 900 seconds. Since C programs and spreadsheets prefer radians, that is 6.5629E-02 radians in 900 seconds, or 7.292116e-5 radians per second. Launching after prime orbit launch results in an apogee to the east of construction and prime orbit apogees (0 degrees or 0 radians); counterclockwise when viewed from the north. Launching before the prime orbit results in an apogee to the west (clockwise when viewed from the north). If the apogee of the candidate orbit is above the radius of the construction orbit at that angle, the two ellipses will intersect at two angles, west and east of that angle.
The candidate launch orbit (hereafter, for this section, launch orbit ) will have less angular momentum than the construction orbit, so at the same radius, the westward orbital velocity will be somewhat less (roughly 10% less, approximately 100 m/s). That means that where the orbits cross, the curvatures differ; the launch orbit is more curved, and the radial accelerations are higher for the launch orbit. Near the west crossing, the vehicle in the launch orbit will be descending (towards the Earth) relative to the construction orbit; near the east crossing, the vehicle ascends.
The construction orbit perigee is (8378-6458) = 1920 km higher than the candidate launch orbit perigee (at launch loop exit), so if the launch orbit is less than 1920 km higher than the construction orbit, 77870 km, it will have a smaller semimajor axis and less than a 86164 second orbital period. By choosing a higher or lower launch apogee, we can choose a longer or shorter delay to apogee. A higher launch velocity results in a higher apogee, and (paradoxically) a longer time delay to apogee.
The western crossing is after launch orbit apogee; given the slower launch orbit velocity, it will cross the construction orbit much later. The eastern crossing is before launch orbit apogee, and will cross the construction orbit sooner.
Our goal is to adjust the launch velocity (and thus the apogee radius) so that the launch vehicle crosses the construction orbit at the same time as the construction station.
More or less; we will actually aim at a capture device a few kilometers from a docking station attached to the construction station. The actual positions will be tweaked for apsidal precession, lunar and solar and Jupiter tides, and a host of other small perturbations. Those tweaks will be computed and optimized for large ensembles of launched vehicles during launch scheduling; this keeps trajectory design engineers employed, and gives the programmers something to miscalculate . "Orbital accountant" will be a job title, and they will not graduate from orbital accountancy school until they risk their lives travelling an outbound and homewards trajectory they've designed for themselves.
So, given the parameters for a launch orbit:
perigee time t_{pl} (the launch time, approximately 12 hours before construction orbit apogee time t=0
apogee angle \theta_l (determined by launch or perigee time)
apogee radius r_{al} (determined by launch velocity)
perigee radius r_{pl} = 6458 km, (launch loop exit radius)
semimajor axis a_{l} = ( r_{pl} + r_{al} ) / 2
eccentricity e_l = ( r_{pl} - r_{al} ) / ( r_{pl} + r_{al} ) ~~ negative, relative to apogee
... we have fully specified the launch orbit, and can compute orbit period, apogee time, and the position and velocity vectors as a function of time.
The construction orbit is specified as before; perigee radius r_{pc} = 8378 km, apogee radius r_{ac} = 75950 km, apogee angle \theta_c = 0°, apogee time t_{ac} = 0 seconds. As before, the calculated period is t_c = 86164.0905 seconds (one sidereal day), the semimajor axis is a_c = 42164 km (the same as a circular geostationary orbit), and the eccentricity is e_c = -0.8013 .
Given the orbit apogee angle \theta_c, we can compute the minimum launch apogee radius from the construction orbit radius at that angle:
r_{al-min} = (1 - {e_c}^2 ) r_{ac} / ( 1 + e_c \cos( \theta_c ) )
The minimum radius is r_{ac} \approx 75950 km for t_{ac} = 0 seconds. For a launch time 900 seconds before prime launch, the launch angle \theta_L is ( 2 \pi 900 / 86164 ) = -0.065629 radians = +3.76° (positive is westward in the eastern Pacific), and for 900 seconds after prime launch, the launch angle \theta_L is -0.065629 radians = -3.76°. In either case, the minimum launch apogee radius must be greater than 75296 km to insure two intersection points.
Given the coplanar launch and construction ellipses, we can:
- 1 compute the two angles at which they intersect
- 2 For each intersection angle, we can compute the arrival time of both the construction station and the launch vehicle
- 3 Then we raise or lower apogee and recompute.
Repeat until one of the time differences is zero; that defines a rendezvous in time and space.
Let's guess a candidate launch orbit apogee r_{al} . We will find the angles where the radii of launch and construction orbits are the same by equating their radii and computing the angles \theta_1, \theta_2 where this happens.
Starting with r_{pc} = 8378 km, r_{ac}=75950 km, r_{pl} = 6458 km, and r_{al} TBD, compute the semimajor axes and eccentricities of the two orbits:
a_c = 42164 km |
e_c = -0.8013 |
a_l = 3229 km + 0.5*r_{al} |
e_l = ( 6458 km - r_{al} ) / ( 6458 km + r_{al} ) |
For the radius to be the same, the two orbit equations are equal:
r(\theta) ~=~ { \Large { { ( 1 - {e_L}^2 ) ~ a_L } \over { 1 + e_L \cos( \theta + \theta_L ) } } } ~=~ { \Large { { ( 1 - {e_c}^2 ) ~ a_c } \over { 1 + e_c \cos( \theta ) } } }
Solving for \theta is easier if we define some intermediate variables:
K ~=~ { \Large { { ( 1 - {e_L}^2 ) ~ a_L } \over { ( 1 - {e_c}^2 ) ~ a_c } } }
So that K ( 1 + e_c \cos( \theta ) ~=~ 1 + e_L \cos( \theta + \theta_L )
Rearranging: ( K - 1 ) /e_L + ( K e_c / e_L ) \cos( \theta ) ~=~ \cos( \theta + \theta_L )
"Take apart" the second cosine term using the trig identity \cos( x + y ) = \cos(x) \cos(y) - \sin(x) \sin(y) :
( K - 1 ) / e_L + ( K e_c / e_L ) \cos( \theta ) ~=~ \cos( \theta ) \cos( \theta_L ) - \sin( \theta ) \sin( \theta_L )
We know all the terms except \cos( \theta ) and \sin( \theta ) ; how do we solve for those? Let's solve the equation for \cos( \theta ) in terms of the constants and \sin ( \theta ) :
\cos( \theta_L ) \cos( \theta ) - ( K e_c / e_L ) \cos( \theta ) ~=~ ( K - 1 ) / e_L + \sin( \theta_L ) \sin( \theta )
~~~~ ( \cos( \theta_L ) - ( K e_c / e_L ) ) \cos( \theta ) ~=~ ( K - 1 ) / e_L + \sin( \theta_L ) \sin( \theta )
~~~~ \cos( \theta ) ~=~ \Large { { ( K - 1 ) / e_L + \sin( \theta_L ) \sin( \theta ) } \over { \cos( \theta_L ) - ( K e_c / e_L ) } }
Define two more intermediate variables:
M = { \Large { { K - 1 } \over { e_L \cos( \theta_L ) - K e_c } } } and N = { \Large { \sin( \theta_L ) \over { e_L \cos( \theta_L ) - K e_c } } }
The below may be invalid ...
The first table describes a series of increasingly higher altitude construction orbits, with periods that are multiples of sidereal days, synchronizing the orbit with the launch loop as it rotates below.
The first "one sidereal day" orbit will be convenient for the construction of space solar power satellites in synchronous orbits. There may be as many as 96 construction station orbits, spaced around the "sidereal clock", fed by separately synchronous vehicle streams. The quickest return to Earth (say, to return a stabilized accident victim to a hospital on earth) will be from the one day orbit, with a 116 m/s retrograde delta V to drop perigee to 50 km reentry altitude in the atmosphere.
If the launch is not captured by the construction station, the vehicle will be in a shorter period orbit and not synchronize with the station on it's next pass. A 2 m/s retrograde delta V from the launch orbit will also drop perigee to 50 km reentry altitude.
The larger long period orbits will be suitable for the launch of interplanetary missions. They will suffer from larger tidal effects from the moon. Still, if the 5 day orbit at the bottom can be made to work, then the capture delta V will be a mere 40 m/s. After an interplanetary vehicle is assembled in this long period construction orbit and the interplanetary trajectory window opens, the perigee is lowered to perhaps 222 km altitude (6600 km radius) with a 37 m/s thrust at apogee. 2.5 days later at perigee, a 1.1 km/s delta V rocket burn can launch the assembled vehicle with an escape velocity excess vinf = 5.5 km/s for rendezvous with an Aldrin Mars cycler orbit. The referenced paper is for a spartan 75 metric tonne cycler; a vehicle built with 20 4-tonne additions every 5 days over a two year period could mass more than 10,000 tonnes, adding to a vast, shielded wheel suitable for centuries of continuous occupation.
construction orbit perigee radius = 8378 km launch orbit perigee radius = 6458 km |
||||||||
sid |
period sec |
apogee |
apogee velocity m/s |
perigee velocity m/s |
||||
day |
constr. |
launch |
radius km |
constr. |
launch |
diff. |
constr. |
loop |
1 |
86164 |
83238 |
75950.3 |
1021.18 |
906.95 |
114.23 |
9257.46 |
10665.8 |
2 |
172328 |
168634 |
125484.9 |
630.56 |
557.63 |
72.94 |
9444.51 |
10834.8 |
3 |
258492 |
254260 |
167032.0 |
477.45 |
421.50 |
55.95 |
9518.89 |
10901.3 |
4 |
344656 |
339996 |
204116.1 |
392.41 |
346.09 |
46.32 |
9560.46 |
10938.4 |
5 |
430820 |
425798 |
238199.6 |
337.21 |
297.22 |
39.99 |
9587.54 |
10962.5 |
6 |
516985 |
511647 |
270068.0 |
298.02 |
262.56 |
35.46 |
9606.82 |
10979.6 |
7 |
603149 |
597528 |
300205.2 |
268.51 |
236.48 |
32.03 |
9621.36 |
10992.5 |
8 |
689313 |
683436 |
328935.4 |
245.35 |
216.02 |
29.32 |
9632.79 |
11002.6 |
9 |
775477 |
769364 |
356489.5 |
226.60 |
199.47 |
27.13 |
9642.05 |
11010.8 |
10 |
861641 |
855309 |
383039.5 |
211.06 |
185.76 |
25.30 |
9649.73 |
11017.6 |
Multiple vehicles per day to a one sidereal day construction orbit
Vehicles can be launched from the launch loop into higher apogee orbits over a ±15 minute window around the prime orbit; they will arrive with a bit more tangential velocity, and far more radial velocity. Vehicles launched before the prime orbit time will arrive with downward radial velocity; vehicles launched after the prime orbit time will arrive after.
I presume the vehicles will be as cheap and as close to passive and uncomplicated as possible, and will arrive near the station to be captured by an active maneuvering tether. They will be perturbed somewhat by the turbulent passage out of the thin remaining atmosphere after a high precision launch by the launch loop, and some way of thrusting them into an exactly precise (to millimeters absolute position and micrometers per second relative velocity) is needed.
- My guess is that a reusable "thrust panel" consisting of a "black" laser-absorbent refractory metal sponge, saturated (somehow) with lots of hydrogen, can be designed to eject that hydrogen at very high velocity when it is smacked by a laser pulse from an orbiting station with very big optics. Wild handwaving, but a laser propulsion genius like Leik Myrabo or Jordin Kare might have some good ideas about how to do it (is Jordin's collaborator Tom Nugent still active?). Depleted panels will be returned to Earth for recharge and resuse, or used for bulk shielding.
An ablative rubber panel might also work, but might scatter too much material into persistent retrograde orbits and pollute low-Earth orbital space with ram-surface-eroding material. A '''No Gram Left Behind''' policy will be necessary for a permanent gigatonne/year spacefaring civilization.
When vehicles arrive, they will be "lassoed" by a velocity and position-matched loop on a deployed cable. They will pull the cable against a drum and a generator, producing electricity to drive some form of propulsion attached to the station itself. We do not need high ISP, but propellant plumes with tightly constrained velocity profiles will be designed to launch all of the propellant into a retrograde orbit with a perigee below the top of the atmosphere. An inert material like argon might be best, so it does not upset upper atmosphere chemistry too much.
Vehicles arriving at ± 900 seconds will have radial velocities around 660 m/s, or excess kinetic energies of 220 kJ/kg; if that was converted 50% efficiently to propellant kinetic energy, a 10% propellant fraction could be launched retrograde at ≈1500 m/s, an impulse of 150 kg-m/s per vehicle kg, more than enough to restore station momentum "lost" to a vehicle capture. A 670 m/s combined radial and tangential velocity capture is frightening - or perhaps "a mere engineering detail", as antimatter propulsion advocate Bob Forward was fond of saying. I don't know how to do it safely and reliably, but someone reading this may be inspired to learn how.
If the expelled propellant fraction is 10%, we can guess that perhaps one in eight of the incoming vehicles are tankers delivering liquid (argon?) propellant.
Note: There may be a computation error in this table, and half the values invalid; MoreLater
capture8 construction perigee = 8378 km launch perigee = 6458 km |
||||||
sidereal period 1 days 86164.091 seconds |
||||||
|
period |
arrival |
apogee |
velocity change m/s |
||
|
sec |
sec |
km |
tangent |
radial |
plane |
construction |
86164.091 |
0.000 |
75950.339 |
0.00 |
0.00 |
0.00 |
prime cargo |
83238.210 |
0.000 |
75950.339 |
114.23 |
0.00 |
0.00 |
|
||||||
-900s cargo |
84605.521 |
1954.218 |
76850.339 |
114.20 |
-652.63 |
-8.32 |
-870s cargo |
84559.824 |
1888.892 |
76820.339 |
114.58 |
-635.58 |
-8.04 |
-840s cargo |
84514.136 |
1823.650 |
76790.339 |
114.92 |
-618.17 |
-7.76 |
-810s cargo |
84468.456 |
1758.485 |
76760.339 |
115.24 |
-600.39 |
-7.47 |
-780s cargo |
84422.784 |
1693.397 |
76730.339 |
115.52 |
-582.27 |
-7.19 |
-750s cargo |
84377.120 |
1628.381 |
76700.339 |
115.78 |
-563.81 |
-6.91 |
-720s cargo |
84331.465 |
1563.434 |
76670.339 |
116.01 |
-545.01 |
-6.63 |
-690s cargo |
84285.818 |
1498.552 |
76640.339 |
116.21 |
-525.89 |
-6.35 |
-660s cargo |
84240.179 |
1433.733 |
76610.339 |
116.38 |
-506.46 |
-6.08 |
-630s cargo |
84194.548 |
1368.972 |
76580.339 |
116.53 |
-486.72 |
-5.80 |
-600s cargo |
84148.925 |
1304.267 |
76550.339 |
116.65 |
-466.69 |
-5.52 |
-570s cargo |
84103.311 |
1239.613 |
76520.339 |
116.75 |
-446.37 |
-5.24 |
-540s cargo |
84057.705 |
1175.006 |
76490.339 |
116.82 |
-425.77 |
-4.96 |
-510s cargo |
84012.108 |
1110.444 |
76460.339 |
116.87 |
-404.91 |
-4.69 |
-480s cargo |
83966.518 |
1045.920 |
76430.339 |
116.89 |
-383.78 |
-4.41 |
-450s cargo |
83920.937 |
981.432 |
76400.339 |
116.89 |
-362.41 |
-4.13 |
-420s cargo |
83875.364 |
916.973 |
76370.339 |
116.87 |
-340.80 |
-3.86 |
-390s cargo |
83829.799 |
852.539 |
76340.339 |
116.82 |
-318.96 |
-3.58 |
-360s cargo |
83784.243 |
788.123 |
76310.339 |
116.76 |
-296.90 |
-3.31 |
-330s cargo |
83738.695 |
723.717 |
76280.339 |
116.67 |
-274.62 |
-3.03 |
-300s cargo |
83693.155 |
659.313 |
76250.339 |
116.55 |
-252.13 |
-2.76 |
-270s cargo |
83647.623 |
594.901 |
76220.339 |
116.42 |
-229.43 |
-2.48 |
-240s cargo |
83602.099 |
530.467 |
76190.339 |
116.26 |
-206.53 |
-2.20 |
-210s cargo |
83556.584 |
465.991 |
76160.339 |
116.09 |
-183.41 |
-1.93 |
-180s cargo |
83511.077 |
401.450 |
76130.339 |
115.89 |
-160.06 |
-1.65 |
-150s cargo |
83465.579 |
336.806 |
76100.339 |
115.67 |
-136.45 |
-1.38 |
-120s cargo |
83420.088 |
271.998 |
76070.339 |
115.43 |
-112.53 |
-1.10 |
-90s cargo |
83374.606 |
206.920 |
76040.339 |
115.16 |
-88.17 |
-0.83 |
-60s cargo |
83329.132 |
141.354 |
76010.339 |
114.88 |
-63.12 |
-0.55 |
-30s cargo |
83283.667 |
74.699 |
75980.339 |
114.57 |
-36.64 |
-0.28 |
prime cargo |
83238.210 |
0.000 |
75950.339 |
114.23 |
0.00 |
0.00 |
30s cargo |
83320.865 |
38.227 |
76004.884 |
114.85 |
42.83 |
0.28 |
60s cargo |
83391.496 |
58.300 |
76051.481 |
115.35 |
70.36 |
0.55 |
90s cargo |
83460.787 |
76.359 |
76097.180 |
115.82 |
96.13 |
0.83 |
120s cargo |
83529.441 |
93.472 |
76142.446 |
116.25 |
121.03 |
1.10 |
150s cargo |
83597.718 |
110.036 |
76187.451 |
116.67 |
145.38 |
1.38 |
180s cargo |
83665.745 |
126.248 |
76232.280 |
117.05 |
169.32 |
1.65 |
210s cargo |
83733.595 |
142.218 |
76276.980 |
117.42 |
192.94 |
1.93 |
240s cargo |
83801.311 |
158.017 |
76321.580 |
117.76 |
216.28 |
2.20 |
270s cargo |
83868.922 |
173.693 |
76366.098 |
118.08 |
239.37 |
2.48 |
300s cargo |
83936.447 |
189.278 |
76410.548 |
118.38 |
262.21 |
2.75 |
330s cargo |
84003.901 |
204.798 |
76454.939 |
118.65 |
284.81 |
3.02 |
360s cargo |
84071.292 |
220.271 |
76499.278 |
118.91 |
307.17 |
3.30 |
390s cargo |
84138.630 |
235.712 |
76543.568 |
119.14 |
329.30 |
3.57 |
420s cargo |
84205.918 |
251.132 |
76587.815 |
119.35 |
351.18 |
3.85 |
450s cargo |
84273.161 |
266.542 |
76632.020 |
119.53 |
372.81 |
4.12 |
480s cargo |
84340.361 |
281.950 |
76676.185 |
119.69 |
394.19 |
4.40 |
510s cargo |
84407.520 |
297.362 |
76720.312 |
119.83 |
415.30 |
4.67 |
540s cargo |
84474.641 |
312.786 |
76764.401 |
119.95 |
436.14 |
4.95 |
570s cargo |
84541.722 |
328.226 |
76808.454 |
120.04 |
456.71 |
5.22 |
600s cargo |
84608.766 |
343.687 |
76852.469 |
120.11 |
476.98 |
5.50 |
630s cargo |
84675.771 |
359.174 |
76896.449 |
120.15 |
496.96 |
5.78 |
660s cargo |
84742.738 |
374.691 |
76940.391 |
120.17 |
516.63 |
6.05 |
690s cargo |
84809.667 |
390.241 |
76984.296 |
120.16 |
535.98 |
6.33 |
720s cargo |
84876.555 |
405.828 |
77028.164 |
120.12 |
555.01 |
6.61 |
750s cargo |
84943.403 |
421.455 |
77071.993 |
120.06 |
573.71 |
6.88 |
780s cargo |
85010.210 |
437.125 |
77115.784 |
119.97 |
592.08 |
7.16 |
810s cargo |
85076.973 |
452.842 |
77159.535 |
119.85 |
610.09 |
7.44 |
840s cargo |
85143.693 |
468.608 |
77203.246 |
119.71 |
627.75 |
7.72 |
870s cargo |
85210.366 |
484.425 |
77246.915 |
119.53 |
645.04 |
8.00 |
900s cargo |
85276.992 |
500.297 |
77290.542 |
119.33 |
661.97 |
8.28 |
Launching from a Construction Orbit to Mars
If an interplanetary vehicle accelerates to a velocity greater than escape velocity at that orbital radius, the excess velocity results in a velocity distant from the Earth of v_{\infty} . ( "Vee infinity" ). For interplanetary launches, higher v_{\infty} reduces travel time, though the arrival velocity can be quite high, which leads to high gee loads entering atmospheres like Mars. The escape velocity at a radius r from the Earth is v^2_{esc} = 2 \mu / r . Given a starting velocity v_r at radius r , v_{\infty}^2 = v_r^2 - 2 \mu / r A small change in v_r can lead to a sizable v_{\infty} . v_{\infty}
Interplanetary launches will begin at the construction orbit apogee, with a small retrograde thrust to drop perigee close to the Earth, maximizing perigee velocity. That maximizes the boost effect of additional delta V at perigee; think of this as adding gravitational energy to the expended propellant, and maximizing the effect of thrust, thus achieving the desired v_{\infty} with a minimum of propellant.
A Hohmann trajectory to Mars results from a v_{\infty} of 2.9 km/s. Consider an "injection orbit" starting at the one day construction orbit apogee radius r_{ca} = r_{ia} = 75950 km. with an injection orbit perigee radius r_{ip} of 6500 km ( 122 km altitude ), for a semimajor axis s_a = 41225 km. The The injection orbit apogee velocity v_{ia} is:
v_{ia} = \sqrt{ ( \mu / s_i ) ( r_{cp} / r_{ca} ) } = \sqrt{ ( 398600.4418 / 41225 ) ( 6500 / 75950 ) } = 0.90966 km/s
The retrograde injection delta V is 1021.18-909.66 m/s = -112 m/s. Hopefully, the propellant plume will have a velocity on the order of -1600 m/s, so that most of the exhaust plume molecules re-enter, that is, have a perigee within the upper atmosphere. 1600 m/s is a low Isp < 180 seconds, so this should use a low energy propellant mixture, perhaps LH/LOX with injected water.
The insertion orbit perigee velocity is:
v_{ip} = \sqrt{ ( \mu / s_c ) ( r_{ca} / r_{cp} ) } = \sqrt{ ( 398600.4418 / 41225 ) ( 75950 / 6500 ) } = 10.629 km/s
for a v_{\infty} = 2.9 km/s, the interplanetary transfer orbit perigee velocity v_{tp} must be:
v_{tp} = \sqrt{ 2.9² + 2 × 398600.4418 / 6500 } = 11.448 km/s
So, the delta V from insertion orbit perigee to an interplanetary Hohmann to Mars is ( 11448 - 10629 ) = 819 m/s, and the total delta V starting at construction orbit perigee is 931 m/s, far less than the ≈11000 m/s direct rocket launch velocity from the Earth's surface.
A more likely mission launched from a construction is a massive 300,000 tonne "Island One" scale permanent habitat, constructed over years in the construction orbit, launched with a crew of thousands into a permanent Aldrin Mars Cycler orbit. This crew will finish the construction and preparation of a Deimos colonization mission, and deploy hundreds of entry vehicles to slow down in the thin Mars atmosphere, and orbit Mars until they are sequentially landed on Deimos to construct a large Deimos scientific colony. A subsequent visit years later by the Mars cycler will deliver the crew for the colony and their supplies.
Mars and Earth have a 782.0 sidereal day (779.9 solar day) synodic period; that enough time to deliver 32,000 vehicles, or 160,000 vehicle-tonnes to more than 40 construction stations from a minimum-sized fully-powered (6GW) launch loop. That is many HUGE missions. A Deimos Synchrotron Light Source ("DSLS") facility for advanced bio-molecule characterization might weigh 30,000 tonnes. A 6 GW space solar power satellite might mass 60,000 tonnes at 10 kg/kW, and require an additional 60,000 tonnes for deployment systems, stationkeeping, and propellant.