Re-entry Heating Calculator
Compute stagnation point heat flux and surface temperature during atmospheric entry for Earth, Mars, or Venus using the Sutton-Graves correlation.
🔥 What is Re-entry Heating?
Re-entry heating is the aerodynamic heat generated when a spacecraft descends through a planetary atmosphere at hypersonic speed. As the vehicle compresses the atmosphere ahead of it, kinetic energy converts to thermal energy in the shock layer, raising gas temperatures to thousands of Kelvin. A fraction of this heat is transferred to the vehicle surface by convection (and, at very high speeds, by radiation from the glowing plasma). The rate of heat transfer per unit area is the stagnation heat flux, measured in kilowatts per square meter or watts per square centimeter. Designing a thermal protection system (TPS) that can survive peak heating is one of the most demanding engineering challenges in spacecraft development.
The Sutton-Graves correlation is the standard formula for estimating stagnation point heat flux at hypersonic speeds: q_s = k times sqrt(rho/R_n) times V^3. It was derived empirically from wind tunnel data at Mach 10 to 20 in dissociating air and validated against flight measurements from Mercury, Gemini, and Apollo. The formula shows that heat flux increases as the cube of entry velocity: doubling velocity increases stagnation heating eightfold. It also decreases inversely with the square root of nose radius, which is why blunt capsule designs are used instead of sharp cones for maximum-heating entries.
Practical re-entry vehicles range from the blunt Apollo Command Module (nose radius 4.7 m, peak heating about 50 W/cm²) to the Space Shuttle orbiter (nose cap radius 0.3 m, wing leading edges 0.04 m, peak leading edge heating 60-80 W/cm²) to Mars EDL aeroshells (2.25 m radius on 4.5 m shell, peak heating 150-200 W/cm²). On the extreme end, the Galileo probe entered Jupiter at 47 km/s and experienced peak stagnation heating of 200,000 W/cm², requiring a carbon phenolic heat shield that ablated away 80 kg of material during the entry.
This calculator implements the Sutton-Graves formula with correct atmospheric models for three planetary bodies. Earth uses the full 7-layer ISA 1976 standard atmosphere with temperature gradients and pressure derived from hydrostatics. Mars uses an exponential model with surface density 0.020 kg/m³ and scale height 11.1 km, consistent with Viking Lander and MER entry measurements. Venus uses an exponential model with surface density 64 kg/m³ and an 17 km scale height. The Wall Temperature mode solves the radiation equilibrium to estimate how hot the heat shield surface gets and recommends a TPS material class.
📐 Formula
qs = k × √(ρ∞ / Rn) × V∞3
qs = stagnation heat flux (W/m²)
k = Sutton-Graves constant: 1.7415 × 10-4 for air (Earth), 1.9027 × 10-4 for CO2 (Mars, Venus)
ρ∞ = freestream atmospheric density at entry altitude (kg/m³)
Rn = nose radius of the vehicle (m). Larger Rn = lower heating
V∞ = freestream entry velocity (m/s)
Wall temperature: Tw = (qs / (ε × σ))0.25
ε = surface emissivity (0 to 1; 0.85 for ceramic TPS)
σ = 5.67 × 10-8 W/(m²·K4) Stefan-Boltzmann constant
Example: At 75 km Earth altitude (rho = 3.5e-5 kg/m³), V = 7.8 km/s, Rn = 0.5 m: q = 1.74e-4 × 0.00836 × 4.75e11 = 690 kW/m² = 69 W/cm²
📖 How to Use This Calculator
Steps
1
Select the entry body and set trajectory - Choose Earth, Mars, or Venus from the body menu. The calculator automatically sets default altitude and velocity for that planet. Adjust altitude in km and entry velocity in km/s for your specific trajectory.
2
Set the nose radius - Enter the vehicle nose radius in meters. Blunt capsules like Apollo use 4-5 m; generic spacecraft use 0.3-1 m; missile warheads use 0.05-0.2 m. Larger radius reduces heating.
3
Read stagnation heat flux results - In Stagnation Heating mode, read the heat flux in kW/m² and W/cm², the dynamic pressure, the local atmospheric density, and the qualitative heating context comparing to known vehicles.
4
Switch to Wall Temperature mode - In Wall Temperature mode, set the surface emissivity (0.85 for typical ceramics) and read the equilibrium surface temperature in Kelvin and Celsius, along with a TPS material recommendation.
💡 Example Calculations
Example 1 — Earth LEO Re-entry (ISS-like Return)
Earth, 75 km altitude, 7.8 km/s, nose radius 0.5 m
1
ISA 1976 density at 75 km: using layer 7 (71 km, T=214.65 K, L=-0.002 K/m, P=3.956 Pa). At 75 km: T=206.6 K, P=2.07 Pa, rho = 2.07/(287.05 x 206.6) = 3.49e-5 kg/m³.
2
Stagnation flux: q = 1.7415e-4 x sqrt(3.49e-5/0.5) x 7800^3 = 1.7415e-4 x 0.00835 x 4.747e11 = 6.90e5 W/m².
3
q = 690 kW/m² = 69.0 W/cm². This falls in the "high heating" Orion/Apollo CM nose range. Dynamic pressure = 0.5 x 3.49e-5 x 7800² = 1.06 kPa.
Stagnation heat flux = 690 kW/m² (69 W/cm²), Dynamic pressure = 1.06 kPa
Try this example →Example 2 — Apollo Lunar Return Entry
Earth, 75 km altitude, 11.0 km/s, nose radius 4.7 m (Apollo CM)
1
Density at 75 km: 3.49e-5 kg/m³ (same as example 1). Lunar return entry at 11 km/s = 11,000 m/s.
2
q = 1.7415e-4 x sqrt(3.49e-5/4.7) x 11000^3 = 1.7415e-4 x 2.73e-3 x 1.331e12 = 6.32e5 W/m².
3
q = 632 kW/m² = 63.2 W/cm². The large 4.7 m Apollo nose radius substantially reduces heating despite the much higher velocity. Historical Apollo peak was 40-50 W/cm² (this estimate is conservative for the actual re-entry altitude).
Stagnation heat flux = 632 kW/m² (63.2 W/cm²), consistent with Apollo AVCOAT TPS design
Try this example →Example 3 — Mars Science Laboratory EDL
Mars, 15 km altitude, 5.9 km/s, nose radius 2.25 m (Curiosity aeroshell)
1
Mars density at 15 km: rho = 0.020 x exp(-15000/11100) = 0.020 x 0.259 = 5.18e-3 kg/m³. Much denser than at comparable Earth altitudes despite thin surface air.
2
q = 1.9027e-4 x sqrt(5.18e-3/2.25) x 5900^3 = 1.9027e-4 x 0.0479 x 2.054e11 = 1.872e6 W/m².
3
q = 1872 kW/m² = 187 W/cm². Historical MSL peak was about 200 W/cm², confirming the formula. This required the PICA heat shield on Curiosity and a large PICA-X shield on Perseverance.
Stagnation heat flux = 1872 kW/m² (187 W/cm²), matching historical MSL/Perseverance measurements
Try this example →Example 4 — Earth Re-entry Wall Temperature
Earth, 75 km, 7.8 km/s, nose 0.5 m, emissivity 0.85 - equilibrium wall temperature
1
Stagnation heat flux from Example 1: q_s = 690 kW/m² = 690,000 W/m².
2
Equilibrium temperature: T_w = (690,000 / (0.85 x 5.67e-8))^0.25 = (690,000 / 4.82e-8)^0.25 = (1.432e13)^0.25 = 1945 K.
3
T_w = 1945 K = 1672 degrees C. This exceeds the limit of reusable ceramic tiles (about 1870 K) and calls for ablative PICA or phenolic-carbon TPS. Radiation cooling is 100% effective at equilibrium.
Equilibrium wall temperature = 1945 K (1672 C), requires ablative PICA TPS
Try this example →❓ Frequently Asked Questions
What is stagnation point heat flux during re-entry?
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 times sqrt(rho/R_n) times 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 heat flux at the same nose size.
What is the Sutton-Graves formula and when is it valid?
The Sutton-Graves correlation q_s = k times sqrt(rho/R_n) times V^3 was derived empirically from hypersonic wind tunnel data and validated against flight measurements from Mercury, Gemini, and Apollo. It is most accurate for blunt bodies (R_n greater than 0.1 m) at Mach 10 to 30 in dissociating air or CO2. It overestimates heating for slender bodies and does not include radiation heating, which becomes significant above about 10-12 km/s. For Earth, k = 1.7415e-4; for CO2 atmospheres like Mars and Venus, k = 1.9027e-4. Results agree with full Navier-Stokes simulations to within 10-15% for typical capsule geometries.
Why does a larger nose radius reduce re-entry heating?
A larger nose radius spreads the stagnation region over a wider area and reduces the velocity gradient, lowering the convective heat transfer coefficient. The formula shows heat flux varies as 1/sqrt(R_n): doubling the nose radius reduces stagnation heating by 29%. Apollo used a 4.7 m radius heat shield, Orion uses 5.0 m, and the Galileo Jupiter probe used a 0.63 m nose on a 1.27 m base. The engineering trade-off is that blunter noses produce more drag, but for entry vehicles drag is beneficial. Sharp-nosed hypersonic gliders (like HTV-2) experience much higher localized heating at the leading edge despite their lower ballistic coefficient.
What is the equilibrium wall temperature of a heat shield?
The equilibrium wall temperature is where the heat shield radiates away exactly as much power as it absorbs. Setting incident flux equal to emitted radiation: q_s = epsilon times sigma times T_w^4. Solving: T_w = (q_s / (epsilon times sigma))^0.25. For 690 kW/m^2 at emissivity 0.85, T_w = 1945 K (1672 C). Higher emissivity lowers T_w because more radiation is emitted at any given temperature. Glazed or coated ceramics have emissivity 0.85-0.92; polished metals can be below 0.2, which traps heat and requires active cooling or ablation.
What TPS materials are used for different heating levels?
Below 750 K wall temperature: reusable ceramic tiles (LI-900, TUFI), used on Space Shuttle windward underbelly (1-30 W/cm^2 range). From 750 to 1200 K: advanced blankets (AFRSI, FRCI-12), used on Shuttle fuselage sides and top. From 1200 to 2000 K: ablative PICA (Phenolic Impregnated Carbon Ablator), used on Stardust capsule (60 W/cm^2), Dragon (40-50 W/cm^2), and Mars landers. Above 2000 K: reinforced carbon-carbon (RCC, used on Shuttle nose and leading edges, to 1870 K surface), AVCOAT (Apollo, to 2500 K ablated), and ultra-high-temperature ceramics (ZrB2, HfB2) for hypersonic aircraft.
How does Mars EDL heating compare to Earth re-entry heating?
Mars EDL produces peak stagnation heating of 150-200 W/cm^2, comparable to Earth lunar return entries, despite Mars atmosphere being 100 times thinner at the surface. This is because entry vehicles reach peak heating at 10-20 km altitude where Mars atmospheric density is 0.003-0.006 kg/m^3, and they arrive at 5.5-7 km/s. The cube dependence on velocity amplifies even moderate speeds. Earth LEO re-entry (7.8 km/s) at 75 km altitude (rho = 3.5e-5 kg/m^3) produces 69 W/cm^2 for a 0.5 m nose and around 12 W/cm^2 for a 5 m nose like Orion, so Mars aeroshells with 2 m noses see 10 to 20 times higher heating than Earth capsules at comparable nose sizes.
What altitude produces peak heating during Earth re-entry?
Peak heating occurs where the product rho times V^2 is maximum during the descent, which for Earth LEO re-entry is around 70-80 km altitude. Below that altitude, the vehicle has decelerated significantly, reducing V^3 faster than density increases. Above that altitude, density is too low to transfer much heat despite high velocity. For ballistic trajectories, peak heating is sharper and at slightly lower altitudes than for lifting trajectories. The Space Shuttle, with its lifting re-entry, spread peak heating over a broader altitude band (50-80 km) at moderate angles of attack, keeping surface temperatures manageable for its tile system.
How is the ISA 1976 standard atmosphere used in this calculator?
The ISA 1976 divides Earth's atmosphere into 7 layers from sea level to 86 km, each with a constant temperature lapse rate. Pressure is computed hydrostatically and density follows from the ideal gas law: rho = P/(R times T). Layer boundaries are at 0, 11, 20, 32, 47, 51, and 71 km. The troposphere (0-11 km) cools at 6.5 K/km; the lower stratosphere (11-20 km) is isothermal at 216.65 K; higher layers have varying gradients. Above 86 km, the ISA is officially undefined; this calculator extrapolates exponentially from the boundary density at 86 km (5.77e-6 kg/m³) with a 7 km scale height.
What is ablation and why is it used in heat shields?
Ablation is the controlled evaporation or charring of the outer heat shield layer during entry. When the surface material vaporizes, it absorbs latent heat equal to its heat of vaporization, removing thermal energy from the surface. The outgassed vapor forms a cool boundary layer that also blocks heat transfer from the hot shock layer. This mechanism allows ablative heat shields to survive heat fluxes far beyond what purely insulative materials can handle. PICA can absorb 1,200 W/cm^2 without burning through; AVCOAT survived 40-50 W/cm^2 on Apollo for the 2,000+ K peak. Once ablated, the material does not regenerate, making ablative TPS single-use.
What is dynamic pressure at entry and why does it matter?
Dynamic pressure q_dyn = 0.5 times rho times V^2 measures the aerodynamic force per unit area acting on the vehicle. At the same altitude and velocity, dynamic pressure drives structural loads on the heat shield, aeroshell ribs, and attachment hardware. For LEO Earth re-entry at 75 km: q_dyn = 0.5 times 3.5e-5 times 7800^2 = 1.06 kPa, a modest load. For Mars EDL at 15 km altitude: q_dyn = 0.5 times 5.2e-3 times 5900^2 = 90.5 kPa, much higher. The entry corridor must be designed so that the maximum aerodynamic load occurs before structural limits are reached, while still decelerating enough to survive.
Why is Venus entry so extreme despite the thick atmosphere?
Venus probes enter at 10-12 km/s (typical interplanetary transfer velocity) into an atmosphere 90 times denser at the surface than Earth and composed of CO2, which has different thermochemical properties. Peak heating occurs at 60-80 km altitude where density is 0.1-0.5 kg/m^3 combined with 10+ km/s velocity. Stagnation heat fluxes of 3,000-10,000 W/cm^2 were experienced by Venera probes, requiring very thick titanium aeroshells and high-performance ablators. The Soviet Venera probes successfully landed using these techniques, surviving on the surface for 23-127 minutes before being destroyed by 462 C temperatures and 90 bar pressure.
How accurate is the re-entry heating calculator?
The Sutton-Graves formula with ISA 1976 density gives stagnation heat flux estimates accurate to within 10-20% for blunt-body capsules in the Mach 10-25 range. Major sources of error: the formula does not include radiation heating (5-20% of total at LEO re-entry, dominant above 12 km/s); the ISA model deviates from actual atmospheric conditions by 10-30% at high altitudes; real vehicles do not fly purely ballistic trajectories; and wake heating, shoulder heating, and after-body heating are not captured. Use these results as first-order estimates for feasibility studies and educational analysis. Production re-entry TPS sizing requires full coupled aeroheating and ablation simulations.