Size: 2021
Comment:
|
Size: 2156
Comment:
|
Deletions are marked like this. | Additions are marked like this. |
Line 49: | Line 49: |
This "sticks" at R=1, so the iteration should start with a "known in-between" like $ R = 0.5 * ( 1.0 + r_{GEO}/r_d ) $ MoreLater |
Orbit Circularization
What is the ΔV needed for apogee insertion into a circular equatorial orbit from a launch loop transfer orbit, and for perigee insertion from a space elevator transfer orbit?
The destination orbit has a radius of r_d and a velocity of v_d = \sqrt{ \mu / r_d } were \mu = 398600.4418 km3 / s2.
The launch loop calculation is fairly simple.
An 80 kilometer breech altitude launch loop defines a transfer orbit with a perigee r_p = 6378 + 80 km = 6458 km . The semimajor axis is a = 0.5 * ( r_p + r_d ) , the eccentricity e = ( r_d - r_p ) / ( r_d + r_p ) , the characteristic velocity is v_0 = \sqrt{ \mu / ( a * ( 1 - e^2 ) ) } , and the apogee velocity is v_a = ( 1 - e ) v_0 . Combining and simplifying:
{v_a}^2 = ( 1 - e )^2 {v_0}^2 = ( \mu / a ) ( 1 - e )^2 / ( 1 - e^2 )
{v_a}^2 = ( 2 \mu r_p ) / ( r_d ( r_d + r_p ) )
v_a = \sqrt{ ( 2 \mu r_p ) / ( r_d ( r_d + r_p ) ) }
\Delta V = v_d - v_a
The space elevator calculation is more complicated.
For transfer orbits below geosynchronous, the apogee of the orbit is the release radius r_a , and the angular velocity is \omega_E = 2 \pi / P_{sidereal} , where P_{sidereal} is the sidereal day, 86164.0905 seconds. We want to find r_a and the perigee velocity v_p .
v_a = \omega_E r_a
v_p = r_a v_a / r_d = \omega_E { r_a }^2 / r_d
Orbital energy:
{v_a}^2 - 2 \mu / r_a = {v_p}^2 - 2 \mu / r_d
Solve for r_a :
( \omega_E r_a )^2 - 2 \mu /r_a = ( \omega_E {r_a}^2 / r_d )^2 - 2 \mu / r_d
Multiply by r_a / {\omega_E}^2 {r_d}^3 :
( r_a / r_d )^3 - ( 2 \mu / {\omega_E}^2 {r_d}^3 ) = ( r_a / r_d )^5 - ( 2 \mu / {\omega_E}^2 {r_d}^3 ) ( r_a / r_d )
Normalize for R = r_a / r_d and M = 2 \mu / ( {\omega_E}^2 {r_d}^3 ) .
R^3 - M = R^5 - M * R
R^5 = R^3 + M * (R - 1) ... given M, solve iteratively for R
This "sticks" at R=1, so the iteration should start with a "known in-between" like R = 0.5 * ( 1.0 + r_{GEO}/r_d )