Re-entry Heating Calculator

Stagnation point heat flux is the peak rate of heat transfer per unit area at the nose or leading edge of a spacecraft, where the flow stagnates to zero velocity. The kinetic energy of the hypersonic flow is converted to thermal energy, raising local temperature to several thousand Kelvin. The Sutton-Graves formula q_s = k*sqrt(rho/R_n)*V^3 captures this: heat flux grows as the cube of velocity, meaning doubling entry speed increases heating eightfold. LEO re-entry at 7.8 km/s generates moderate heating (10-100 W/cm^2 depending on nose radius), while lunar return at 11 km/s generates roughly three times more.

The Sutton-Graves correlation is q_s = k*sqrt(rho/R_n)*V^3, where q_s is stagnation heat flux (W/m^2), k is a gas-specific constant (1.74e-4 for air, 1.90e-4 for CO2), rho is freestream density (kg/m^3), R_n is nose radius (m), and V is entry velocity (m/s). It was derived empirically from hypersonic wind tunnel data and validated against flight measurements. The formula is most accurate for blunt bodies at Mach 10-30 in dissociating flow. It overestimates heating for slender bodies and underestimates radiation heating for entries faster than about 12 km/s.

A larger nose radius spreads the stagnation region over a larger area and reduces the velocity gradient at the nose, which lowers the convective heat transfer coefficient. The Sutton-Graves formula shows heat flux varies as 1/sqrt(R_n): doubling the nose radius reduces stagnation heating by about 29%. This is why all entry vehicles use blunt rather than sharp noses: Apollo used a 4.7 m radius, Orion uses a 5.0 m radius, and the Mars Science Laboratory aeroshell had a 2.25 m nose radius on a 4.5 m diameter shell. The penalty is higher drag, but for entry vehicles that drag is beneficial for deceleration.

The equilibrium wall temperature is reached when the heat shield radiates away exactly as much energy as it absorbs: q_s = epsilon*sigma*T_w^4. Solving for T_w: T_w = (q_s / (epsilon*sigma))^0.25. For a stagnation heat flux of 700 kW/m^2 and emissivity of 0.85, this gives about 1945 K (1672 C). This is the theoretical steady-state surface temperature; real TPS materials absorb some heat during the transient entry phase, and ablative materials cool themselves by evaporating. High emissivity (close to 1.0) is crucial for minimizing T_w.

Material selection depends on surface temperature: up to 750 K (477 C), reusable ceramic tiles like Space Shuttle TUFI (Toughened Uni-piece Fibrous Insulation) or LI-900 suffice. From 750 K to 1200 K, advanced blanket systems like AFRSI (Advanced Flexible Reusable Surface Insulation) are used. From 1200 K to 2000 K, ablative materials like PICA (Phenolic Impregnated Carbon Ablator, used on Stardust and Dragon) or SLA-561V are required. Above 2000 K, only carbon-carbon (RCC, used on Space Shuttle leading edges and nose cap) or bulk ablatives (AVCOAT, used on Apollo) can survive.

Earth LEO re-entry at 7.8 km/s produces moderate heating concentrated at 60-80 km altitude where the atmosphere is dense enough to transfer heat but the vehicle is still hypersonic. Mars EDL (Entry, Descent, Landing) at 5.5-7 km/s occurs in an atmosphere 100 times thinner than Earth, but peak heating still reaches 150-200 W/cm^2 because the vehicle enters much deeper (15-25 km) before decelerating significantly. The thin atmosphere provides less drag, so the vehicle stays supersonic longer, spreading heat load over more time. Mars EDL requires robust ablative TPS despite the thin atmosphere.

Stagnation heat flux grows as V^3, so faster is always hotter at the same altitude and nose geometry. Lunar return entries at 11 km/s are roughly 2.7 times hotter than LEO re-entries at 7.8 km/s at the same altitude and nose radius. Planetary probe entries produce extreme heating: Galileo entered Jupiter at 47 km/s and experienced 200,000 W/cm^2 peak heating, requiring a 152 kg carbon phenolic heat shield that ablated 80 kg during entry. Venus probes enter at 10-12 km/s and face a thick atmosphere that decelerates them quickly, but peak heating still reaches 3,000-5,000 W/cm^2.

Ablation is the thermal protection mechanism where the heat shield material intentionally evaporates, sublimes, or chars on the outer surface. The phase change (solid to gas) absorbs large amounts of latent heat, and the outgassed material forms a cool boundary layer that blocks heat transfer to the structure. This is far more effective per kilogram than purely insulative approaches. AVCOAT (used on Apollo and Orion) is a phenolic-filled glass micro-balloons epoxy that chars and ablates. PICA (used on Stardust, MSL, SpaceX Dragon) uses phenolic resin in a carbon fiber preform and can survive 1,200 W/cm^2 without active cooling.

Dynamic pressure q_dyn = 0.5*rho*V^2 measures the aerodynamic force per unit area on the vehicle. Peak dynamic pressure (Max-Q) during entry determines structural loads on the aeroshell, heat shield, and payload. For LEO re-entry at 7.8 km/s and 75 km altitude (rho = 3.5e-5 kg/m^3): q_dyn = 0.5*3.5e-5*7800^2 = 1064 Pa = 1.06 kPa. Much lower than launch Max-Q (which can exceed 20 kPa). For ballistic entries, peak dynamic pressure drives aeroshell structural design and determines terminal velocity before parachute deployment.

The International Standard Atmosphere (ISA 1976) defines temperature, pressure, and density as a function of altitude through 7 layers from sea level to 86 km. Each layer has a linear temperature lapse rate: the troposphere (0-11 km) cools at 6.5 K/km, the lower stratosphere (11-20 km) is isothermal at 216.65 K, and higher layers have varying gradients. Density is derived from pressure and temperature using the ideal gas law: rho = P/(R*T). Above 86 km, the ISA is undefined; this calculator uses an exponential extrapolation from the known 86 km boundary conditions. The ISA accurately predicts Earth entry heating to within 5-10%.

At entry velocities above about 10 km/s, the shock layer plasma ahead of the spacecraft becomes so hot that it radiates significant energy directly to the heat shield surface (radiative heating), in addition to convective heat transfer. The Sutton-Graves formula captures only convective heating. For lunar returns at 11 km/s, radiation adds about 5-15% to total heat flux. For Mars direct entries at 7-8 km/s, it is negligible. For Jupiter probe entries at 47 km/s, radiation dominated: over 90% of the 200,000 W/cm^2 peak flux was radiative. Purpose-built radiation heating models (such as NEQAIR or HARA) are needed for accurate analysis above 12 km/s.

Yes, with caution. The Sutton-Graves formula was derived for blunt bodies during atmospheric entry but gives reasonable estimates for any hypersonic blunt stagnation point above Mach 5. For slender bodies or leading edges (very small R_n), the formula underestimates heating because real small-radius leading edges are closer to the sharp-body limit where heating scales differently. For flat windward surfaces (like the Space Shuttle belly), use the Eckert reference enthalpy method instead. For missiles and re-entry vehicles at Mach 5-15, the Sutton-Graves formula with appropriate nose radius is accurate to within 20-30%.

Convective heating transfers energy through molecule-to-molecule collisions between the hot shock layer gas and the heat shield surface. It dominates below about 10-12 km/s entry velocity and follows the Sutton-Graves V^3 relationship. Radiative heating transfers energy through photon emission from the hot plasma; it scales roughly as V^8-10 and dominates at very high velocities. For planetary probes (Galileo at Jupiter: 47 km/s), radiation overwhelms convection by a factor of 10 or more. For interplanetary returns such as Genesis (returning solar wind samples at 11 km/s), radiation heating was about 20% of the total. Most Earth LEO re-entries are purely convective.