Track Slope

Assume a maximum exit velocity a little larger than Earth escape (11.2 km/s) at 6458 km radius, 80 km altitude, near the equator. At that radius, the Earth's rotation velocity is 0.47 km/s ( 2π × 6458 km / 86414s ), so the maximum atmosphere-relative velocity (ignoring wind and adding drag loss) is 10.8 km/s. We will compute track slope backwards from that.

Loop inclination - the angle to the equator - will probably be between 10 and 30 degrees, To Be Determined.

Atmosphere Density Model

Altitude

Density

Density Scale

Average

km

kg/m3

Height km

km to 80

80 km

1.85e-5

5.89

6.33

70 km

8.28e-5

7.43

6.7

60 km

3.10e-4

8.25

7.1

50 km

1.03e-3

8.03

7.5

40 km

4.00e-3

6.86

7.4

30 km

1.84e-2

6.45

7.24

Linear Heating Profile

a

vehicle acceleration

x

vehicle distance along track

x_e

vehicle exit distance

v

vehicle velocity along track

v_e

vehicle exit velocity

H

vehicle heating (relative)

H_e

vehicle exit heating

t

time from start of launch run

t_e

exit time

\rho

atmospheric density

\rho_e

exit atmospheric density

z

vehicle/track altitude

z_e

exit altitude

An escape velocity run to v_e = 10800 m/s will last t_x = 10800/29.4 = 367 seconds. Most launches will be to high earth orbits, a slightly shorter launch run.

v^2 ~=~ 2 a x

Heating rates, during acceleration and at exit:

H ~~=~ { 1 \over 2 } \rho(z) ~ v^3

H_e ~=~ { 1 \over 2 } \rho(z_e) ~ {v_e}^3

Arbitrarily chosen heating rate over the launch run, ramping up proportional velocity:

H ~~=~ H_e ~ v / v_e

H ~~=~ { 1 \over 2 } \rho(z) ~ v^3

H_e ~=~ { 1 \over 2 } \rho(z_e) ~ {v_e}^3

Thus

H ~=~ { 1 \over 2 } \rho(z) ~ v^3 ~=~ { 1 \over 2 } \rho(z_e) ~ {v_e}^3 ~ ( v / v_e ) ~=~ { 1 \over 2 } \rho(z_e) ~ {v_e}^2 ~ v

Striking common factors from both sides:

\rho(z) ~=~ \rho(z_e) ~ {v_e}^2 / v^2

Substituting v^2 ~=~ 2 a x and v_e^2 ~=~ 2 a x_e :

\rho(z) ~=~ \rho(z_e) ~ 2 ~ a x_e / 2 ~ a ~ x ~=~ \rho(z_e) ~ ( x_e / x )

\rho(z) = \rho(z_e) {v_e}^2 / 2 a x

So, given the exit velocity v_e , the exit density \rho(z_e) , and the distance x along the track from the starting point, we can compute the (higher) density at distance x , and use that to determine the lower altitude at that point.

This oversimplified model blows up as x approaches zero, to infinite \rho(z) density. At some point (relatively close to the start) the track angle dz / dx starts increasing rapidly, where it makes sense to limit the track angle and run at a linear path down to west station, at a lower altitude. Also, the actual distance travelled is larger, (dx^2 + dy^2 )^{1 \over 2} , or dx / cos( i ) where i is the inclination of the track from horizontal.

Over a small range of z around z_0, the function \rho(z) will approximate a simple exponential function:

\rho(z) \approx \rho_0 e^{-(z-z_0)/l_0} ~~~~~ where \rho_0 = \rho(z_0) and l_0 = { { \Large \rho(z_0) } \over { { { \LARGE \partial\rho(z) } \over { \LARGE z } } } { \large ~at~ z~=~ z_0 } } , a log-linear fit, also called the '''density scale height''' around a z_0 .

Density scale height varies between 5.89 km and 8.25 km between 30 and 80 kilometers altitude; logarithmically average that to 7.24 km, spanning the range.

Using the simplified density model, a pure exponential with a H = 7.24 km density scale height, we can compute the change in altitude from the velocity relative to the exit velocity:

\Delta z ~=~ - H ~ ln( \rho(z) / \rho(z_e) ) ~=~ - H ~ ln( {v_e}^2 / 2 a x ) ~=~ - H ~ ln( x_e / x )

More later

Slope:

{ \Large { 1 \over \sin( \theta ) } } ~=~ { \Large { { \partial x } \over { \partial z } } } ~=~ { { \Large \left( { - x_e } \over H \right) } ~ exp { \Large { \left( { { z_e - z } \over H } \right) } } } ~=~ { \Large { x \over H } }

{ \Large \theta } ~=~ arcsin { \Large { \left( { H \over x } \right) } }

It's best not to "waste too much deflection angle" on the ramped portion of the loop between west station and the main launch run; as a wild guess, i = 5 degrees above west station seems sufficient. 1/cos( 5° ) = 1.0038 = dl / dx , good enough for rough estimates. The lower ramped section will have a fairly thick aluminum outer sheath in order to conduct extreme lightning strikes without too much heating or voltage drop.