models.parkinresearch.com

Ballistic Reentry Calculator

This model represents an object that reenters the Earth’s atmosphere. It calculates peak deceleration, peak heating, and the altitudes and times at which that happens. It is an implementation of the first order ballistic reentry model published by C. Gazley Jr of RAND Corporation in 1960, which itself traces back to the Allen-Eggars model published in 1958.

For this model, key approximations simplify the underlying trajectory equations, allowing them to be solved in closed form: The atmosphere’s density is approximated to fall off exponentially with increasing altitude, gravity is approximated as constant regardless of altitude, and the entering body is approximated to have zero lift and a constant coefficient of drag.

A consequence of these approximations is that the entry flight path angle of the body does not subsequently change; this can be the case provided that the angle is 5 degrees or more below horizontal.

Press «enter» to see the model.  To infer results, try one of the examples below, or your own inputs. To clear results, refresh your browser.

Allen and Eggers’ Iron Sphere

On page 2 of their seminal 1958 paper, Allen and Eggers consider the vertical descent of a solid iron sphere having a diameter of 1 foot.

  1.  First, let’s use the same atmospheric parameters that they do by specifying our sea-level atmospheric density as 0.0034 slugs/ft³ and atmospheric scale height as 22,000 ft by typing:
    (CTRL+g)r s.s = 0.0034 «enter»  (CTRL+g)b .f = 22000 «enter»
  2. The body’s ballistic coefficient is implied by the information that it is a 1-foot diameter iron (7,874 kg/m³) sphere with unity drag coefficient:
    B.Sp.D.f = 1 «enter»  B.Sp.(CTRL+g)r = 7874 «enter»  B.C = 1 «enter»
  3. Turning to page 3, let’s take the case where the sphere enters the atmosphere at 40 miles altitude and 30,000 ft/s, vertically (γₑ=-90°):
    he.mi = 40 «enter»  ve.f = 30000 «enter»  g = -90 «enter»
    The model infers that the sphere decelerates at a peak of 234 g’s (peak_decn.-dV/dt) at 43,854 ft altitude (peak_decn.h.ft), and this corresponds to the peak in figure 3 of the paper.
  4. To compare with figure 2 of the paper as well, let’s choose 50,000 ft altitude:
    h.f = 50000 «enter»
    The model infers that velocity at this altitude (V.ft/s) is 20,555 ft, which visually matches the figure.
  5. To estimate heating the model needs to know that the radius of curvature is 0.5 ft:
    R.f = 0.5 «enter»
    The model infers that peak heating (peak_heat.qw) of 66 MW/m² (6,600 W/cm² or 45E5 ft-lbf/s/ft²) occurs at 68,000 ft altitude (peak_heat.h.ft). This enormous heat input would occur if the walls of the sphere were cold, and it would take a blackbody at 5,840 K (peak_heat.Wall.T) to produce the same radiant exitance. Note that this radiant exitance will not exactly match the values in the paper because this model instead uses the stagnation point heating correlation of Kemp & Riddell (1957) to estimate radiant exitance.
run this example

So what about at ground level? Setting altitude to zero:

h.f = 0 «enter»

The radiant sphere of melting iron hits the ground at 233 m/s (V) only 17 seconds after atmospheric entry (t.s)! Note that when the ballistic coefficient is lower, setting an altitude that is too low may result in a consistency error. This is a sign that the model is out of its regime, and instead the body has slowed to terminal velocity. This model does not yet calculate terminal velocity or detect when it is reached, so it is the user’s responsibility to notice when the velocity gets unrealistically low.

run this example

One can also find the altitude at which a velocity (such as terminal velocity) occurs. Clearing the altitude and setting a velocity, for example 300 m/s, does this:

h = «enter»  v = 0.3 «enter»

In this example, the model infers that iron sphere decelerates through 300 m/s at 482 meters altitude.

run this example