Hypervelocity Drag

Drag on an accelerating launch loop vehicle with a hemispherical nose

I am not an aeronautical engineer and probably misunderstand the sources. In any case, the numbers are approximate, and should be treated skeptically.

The drag and heating is acceptable at all altitudes. However, incoming debris impactors in decaying orbits will be more unpredictable when the drag is higher, and more difficult to shield or dodge. Perhaps most of them may be intercepted a few orbit earlier, reducing flux at launch loop track level.

Based on Trajectory Optimization for an Apollo-type Vehicle under Entry Conditions Encountered During Lunar Returm by John W. Young (famous astronaut) and Robert E. Smith Jr., May 1967, NASA TR-R-258, Langley Research Center.

Equations on Page 5 in Foot-second-slug-BTU :

Density in slugs/ft3: multiply kg/m3 by 1.9403203e-3

Power in Btu/ft2-s: multiply by 11350.54 to get W/m2

Velocity in ft/s: divide m/s by 0.3048

Metric equations:

These are for a 1 foot diameter nose, and scale by {r_n}^{-1/2} according to equation 4B-4 on page 520 of Part 4B (Entry Heat Transfer) of the SAE Aerospace Applied Thermodynamics Manual. That sites reference 1, A study of the motion and aerodynamic heating of missiles entering the earth's atmosphere at high supersonic speeds, H. Julian Allen and A. J. Eggers, Jr, NACA TN 4047, 1957. If r_n is in meters, scale by 0.552 {r_n}^{-1/2} .

If we scale these for a half-spherical nose, area \pi {r_n}^2 , we get:

Effective time:

Assume constant acceleration for the vehicle, v = a t , to a maximum velocity V = a T .

define t_{eff} = {\Large { T \over { n+1 } } } = { \Large { V \over { a ( n+1) } } }

If the drag power \dot Q = k v^n = k a^n t^n , then the time integrated power:

Q=k a^n{\Large {T^{n+1}\over {n+1}}}=k a^n T^n{\Large {T\over{n+1}}} = k V^n t_{eff} = \dot Q_{max} t_{eff}

There will also be additional exit or climb-out time for the launch loop added to t_{eff} , TBD. This additional time will be proportionally larger for the radiation fraction, but that will remain small, especially in thinner, higher altitude atmosphere.

The drag losses are much higher; most of the lost energy ends up heating the upper atmosphere (where it radiates efficiently into space, not to the ground). The drag power is P = C_D \rho Area V^3 and the drag loss is P = C_D \rho Area V^3 T/4

Examples

For a 1 meter diameter nose, V=11 km/s, a=3*9.8m/s, T=374 s, CD = 2.0 and density at 80, 100, and 120 km:

altitude km

80

100

120

density kg/m3

1.85e-5

5.60e-7

2.22e-8

\dot Q_c W

3.5e-6

6.1e+5

1.2e+5

\dot Q_r W

1.2e+5

6.2e+2

4.9e+0

exponent

t_{eff}

Q_c J

3.3e+8

5.7e+7

1.1e+7

3

94

Q_r J

2.1e+6

1.1e+4

8.7e+1

20

18

Q_{total} J

3.3e+8

5.7e+7

1.1e+7

heat fraction

1.1e-3

1.9e-4

3.7e-5

drag loss J

1.4e+10

4.4e+8

1.7e+7

drag fraction

4.8e-2

1.5e-3

5.8e-5

heat/drag

0.023

0.13

0.65

Add spreadsheet here.

Validity?

The Young and Smith paper was for Apollo lunar reentry, and the equations may not generalize to lower drag regions. Exit trajectory must still be computed, keeping in mind that Earth's rotation velocity (470 m/s) should be added to the outgoing orbit vector.