HOMESPACE ELEVATOR$4M CHALLENGEABOUT US



Space Elevator Power Analysis and Optimization

Summary

This paper lays out the basics constraints for a Space Elevator power system, performs parameter optimization, and compares the results with real-life technology parameters. The paper also considers the special case of solar climbers that have the additional constraint of a once-per-day launch rate.

Motivation

The Space Elevator is a transportation system, and as these systems go, we’re interested in moving as much payload as possible to orbit for the lowest overall cost. Because the equations governing the Space Elevator are simple, it is possible to make quite a bit of headway using simple mathematical analysis.

Since the Space Elevator is linearly scalable, we normalize the calculations by the maximum mass that is allowed to hang from the bottom of the tether (mmax). Thus a “20-ton” Elevator is one that can support a single 20 ton climber at ground level. Typically, this means that the tether weight 4000 – 6000 tons, and the climbers will actually weigh around 15 tons (since we have multiple climbers simultaneously on the ribbon). Using mmax normalized mass units, the tether weighs 200-300, the climber weighs 0.75, etc. We’ll define a “Standard throughput unit” (STU) as being able to launch one mmax per year. Unless specified otherwise, we’ll be talking about the payload throughput.

To simplify matters, we divide the mass of the climber into payload and power system, assuming the “dead structure” is small in comparison to either of them. (Any structure that scales with the power system (e.g. motors) is incorporated into the overall power density of the power system).

The parameters determined by technology are the power density of the power system, and the maximum speed we can handle the tether. The variables we can tune as part of the design are the climber mass, the power-payload ratio division within the climber, and the time between climber launches.

 

The Scoop:

 

The Scoop:

Ascent Power Profile

The power required to move a mass at a certain velocity is a function of the effective gravity at the altitude that the climber is at: P = m•geff•v, where geff(r) = g(re/r)2 - ω2rm (re=6400 km is the radius of the earth, and ωe = 7.3E-5 rad/sec, the angular velocity of the Earth).

For r < 2re, we can approximate this very well as simply geff(r) = g(re/r)2, which means that for a specific power system, the climber’s velocity will increase according to a square law v(r) = P/(m•g) • (r/re)2 as it moves out, until it reaches some maximum terminal velocity vT determined by the tether handling system. From that point onwards, the climber moves at vT, and the power system is under-utilized.

The formula for the time it takes a climber that is following a constant power velocity profile to cover the distance between re and r is:

Where t0=(r-re)/ve is the time it would have taken the climber to cover the distance if it were moving at a constant velocity, and Q=r/re is the radius ratio.

The distance traveled by a constant power velocity climber (relative to re) is: d/re = (r-re)/re = t/(re/ve-t).

Handoff and Throughput

The spacing between climbers can be characterized by a handoff fraction kH, so that a new climber is launched when the old climber reached geff/g=kH. The handoff altitude rH = re/kH0.5 is the location where this happens, and the time since launch that this happens at is tH)

The mass of each climber can only be (1-kH)mmax, so that the geometrical series (1-kH)+(1-kH)kH+(1-kH)kH2… =1. (This is a slightly conservative, since the spacing between the climbers does not remain constant, and there are only a finite number of climbers).

If the initial parameters are such that rT>rH, (or (vT/ve)0.5 > kH-0.5 ) then the climber will follow a constant-power velocity profile all the way out to the handoff point. We call this a “power limited” profile. Otherwise, the climber will “max-out” on the way to the handoff point, and the profile is called “speed limited”. There is also the possibility the rT < re, which means that the climber power system starts out under-utilized, even at ground level. This is clearly a non-optimal case.

The payload per climber is therefore mP = (1-β)(1-kH) and the mass throughput of the system is P=(1-β)(1-kH)/tH. The more frequently we launch climbers, the smaller each one can be, but the larger the throughput. This trend continues to the limiting case of a continuous (variable speed) belt of cargo, though we see no practical way of doing that. Similarly, the faster we move, the more climbers we can launch, but the larger power systems leave less room for payload.

The parameters dictated by technology are the power density ρP and the terminal velocity vT. The variables we can tune are the handoff constant kH, the power system mass fraction β, and the handoff time tH. The variables represent only 2 degrees of freedom, since clearly once we set β and kH, tH is already determined.

When power limited scenarios, while vT can still be tuned, it does not affect the throughput, since it only alters the behavior of the climber beyond the handoff point.

 

For Example:

 

The Scoop:

Optimization

To find the maximum throughput:
 
The relative payload ismP = (1-kH)(1-β)
And the throughput isP = (1-kH)(1-β)/tH
 

We can now link β and kH through tH

 
Speed-limited:
Initial velocity:ve = β•ρP/g
Terminal point:rT = re(vT/ve)0.5 = re•QT
Terminal altitude:aT = rT-re
Time to terminal point:tT = (re/rT)(rT-re)/ve = (re/ve)(QT-1)/QT
Handoff point:rH = re/kH0.5
Handoff altitude:aH = rH-re
Time to handoff point:tH = tT + (rH-rT)/vT = (re/ve)(QT-1)/QT + (rek-0.5-reQT)/vT =
(re/vT)•[QT(QT-1)-QT + kH-0.5] = (re/vT)[QT(QT-2) + kH-0.5]
Extracting kH:kH = [(tHvT/re) - QT(QT-2)]-2
 
Power-limited:
Initial velocity:ve = β•ρP/g
Handoff point:rH = re(vH/ve)0.5 = re•QH = re•(kH)-0.5
Handoff altitude:aH = rH-re
Velocity at handoff point:vH = veQH2
Time to handoff point:tH = (re/rH)(rH-re)/ve = (re/ve)(QH-1)/QH = (re/vH)QH(QH-1)
Extracting β:β = g•veP = (g/ρP)•(re/tH)•(QH-1)/QH = cH•(QH-1)/QH
Solving for kH:kH = QH-2 = [cH / (cH-β)]-2

It is possible to express dP/dβ is closed form, but the resulting expression is only solvable numerically. For the same effort, it is more interesting to directly optimize P(β).

Results

The worksheet is implemented in MS Excel, and works with the built-in numerical solver to yield optimal values for P in respect to β. Below is one instance of the optimization, for ρP = 1500 and vT = 80. The source worksheet is available online.

Table 1: Optimization Worksheet

The parameter space is 3-dimensional, and we are interested in quite a few of the resulting quantities. The approach taken for aggregating the data is to hold tH constant and plot one table per observed quantity, then experiment with other values of tH.

This process requires a considerable amount of manual work, but gives the experimenter a good insight into the behavior of the system.

Below are the results for daily-cycle operations (tH=86400). Power limited scenarios are shaded. Our focus is on the payload throughput (P). mP and kH are shown for “situational awareness”.

Table 2: Results for daily operations (tH=86400)

The first observation is that once the system becomes power-limited, vT (as expected) no longer influences the result. The second observation is that even before the system becomes power-limited, the performance only advances slowly. If we stay in the “reasonable” vT range of 60-120 m/s, the throughput values are mostly a function of the power density.

As estimated before (using the rH=2re point), the weaker power systems run with mP ≈ 0.25, but we find out that the stronger ones reach much higher, into mP ≈ 0.6. Since even with mP = 0.6 we only have P = 219, let’s look at what can be gained by increasing the launch rates.

Looking at bi-daily operations (tH=43200) and keeping in mind that for the same P, mP will be half its previous value, we get:

Table 3: Results for bi-daily operations (tH=43200)

The shaded region is larger since the smaller climbers need to get out of the way faster and so carry larger power systems, thus maxing out sooner. For this reason while the higher performing system gain up to 50% in throughput, the lower performing systems gain only about 10%. This is to be expected, since there’s little point expediting the launch rate if the system is not capable of getting the climbers far enough out of the way by in half a day.

Conclusions

We can draw the following table, to be used as a rough guide for the throughput available from a power system: (Throughput again is in units of mmax/yr, or STU)

Table 4: Approximate throughput rates as a funciton of power density and launch rate

When compared to the requirement 200-250 STU imposed by the Space Elevator Feasibility condition for 30 MYuri tethers, these results show very little salvation for the low power density systems, since their performance can’t be improved by either a higher vT, or a lower tH. For higher performing systems, as long as we stay with daily operations, even with the benefit of high power density systems and faster travel speeds, throughput remains just under 300 STU.

If we need to get to the 300-400 STU throughput range, the only way to get there is to have a 2500 – 3500 kWatt/kg power system, be able to travel at vT>100 m/s, and get into Bi-Daily operational tempo. Faster-than-daily launches, however, open an entire Pandora’s Box of operational issues pertaining to the day-night cycle and power beaming, very fast atmospheric crossings, etc.

In summary:
- The universe has again conspired to make the Space Elevator feasible, but make us work very hard at it.
- High intensity power beaming which requires cooling radiators is probably not a viable power source.
- 1-sun thin film solar technology is adequate from a power density perspective.
- 1-2 sun power beaming is viable and might be needed to augment solar operations.

References

[1] Artsutanov, Y., "Into the Cosmos by Electric Rocket", Komsomolskaya Pravda, 31 July 1960. (The contents are described in English by Lvov in Science, 158, 946-947, 1967)
[2] Artsutanov, Y., "Into the Cosmos without Rockets", Znanije-Sila 7, 25, 1969
[3] Pearson, J., "The Orbital Tower: A Spacecraft Launcher Using the Earth's Rotational Energy", Acta Astronautica 2, 785-799, 1975
[4] Edwards, B. C., and Westling, E. A., "The Space Elevator: A Revolutionary Earth-to-Space Transportation System", published by the authors, January 2003
[5] B. Shelef, "Space Elevator Power System Analysis and Optimization", The Spaceward Foundation [2008]
[6] B. Shelef, "Space Elevator Calculation Scrapbook", The Spaceward Foundation [2008]
Optimization Worksheet [XLS]
Back to Library

© The Spaceward Foundation 2008 - www.spaceward.org - Mountain View, CA