Reentry Heat Flux⚠ unverified
Aerospace / Structures · Compute the approximate stagnation-point convective heat flux on reentry
Parameters
| Input | Symbol | Unit | Default | Description |
|---|---|---|---|---|
| rho | ρ | kg/m^3 | 1.0 | Free-stream atmospheric density |
| V | V | m/s | 1.0 | Vehicle velocity |
| R_n | Rn | m | 1.0 | Nose radius |
| Output | Symbol | Unit | Description |
|---|---|---|---|
| result | q | W/m^2 | Stagnation-point heat flux, in watts per square metre (W/m^2) |
The science & history
Understanding the Parameters
-
Velocity $V$ — the overwhelming driver, entering as $V^{3.05}$. Re-entry heating is fundamentally the conversion of the vehicle's vast kinetic energy into heat in the shock layer; because kinetic energy goes as $V^2$ and the rate of energy deposition adds roughly another power of $V$, the flux scales as $\sim V^3$. Halving entry speed cuts peak heating nearly eight-fold — which is why lunar-return and interplanetary entries are so much harsher than low-orbit returns.
-
Free-stream density $\rho$ — enters as $\sqrt\rho$. Heating is negligible in the thin upper atmosphere and peaks lower down as density rises — but by then the vehicle has also decelerated, so peak heating occurs at an intermediate altitude, a balance between rising $\rho$ and falling $V$. The entry trajectory is shaped to keep this peak survivable (Density Altitude, ISA Density).
-
Nose radius $R_n$ — appears as $1/\sqrt{R_n}$: a blunter nose (larger $R_n$) sees lower stagnation heating. This is the famous blunt-body insight — counter-intuitively, a rounded, draggy shape is cooler than a sharp one, because it pushes the bow shock forward and spreads the heat over a larger area while radiating more energy back into the flow. Every crewed re-entry capsule is blunt for this reason.
-
The output $q$ — the convective heat flux at the stagnation point (the hottest spot), in W/m². For orbital entry it reaches megawatts per square metre; the TPS must absorb, ablate, or re-radiate it. At the highest speeds a separate radiative heating term (from the glowing shock layer) adds to this convective estimate.
Derivation (Approaching a Proof)
Stagnation-point convective heating comes from boundary-layer theory in the shock layer at the nose. The seminal Fay–Riddell analysis (1958) solved the stagnation-point boundary layer to give the heat flux in terms of the enthalpy difference across it and the velocity gradient at the wall. The key results are:
-
The wall heat flux scales with the stagnation enthalpy $h_0 \approx \tfrac12 V^2$ (the energy available to be dumped into the surface), contributing the dominant velocity dependence.
-
The stagnation-point velocity gradient $du_e/dx \propto \sqrt{2(p_0-p_\infty)/\rho}/R_n \propto V/R_n$, whose square root introduces the $1/\sqrt{R_n}$ blunt-body factor and part of the density dependence.
Combining these boundary-layer scalings gives $q \propto \sqrt{\rho/R_n}\,\cdot\,(\text{function of }h_0)$. Sutton and Graves (1971) reduced the Fay–Riddell physics to a compact engineering power-law fit for air,
$$q = k\,\sqrt{\frac{\rho}{R_n}}\;V^{N}, \qquad k = 1.83\times10^{-4}\ \tfrac{\text{kg}^{1/2}}{\text{m}},\quad N \approx 3.05,$$
with the exponent $3.05$ (rather than exactly $3$) capturing the mild departure from the ideal enthalpy scaling across the real-gas regime. $\blacksquare$ The constant $k$ is specific to Earth air; other atmospheres (Mars CO₂, Venus, the gas giants) use different Sutton–Graves constants.
Dimensional check. The constant $k$ carries the units needed to make $\sqrt{\rho/R_n}\,V^{3.05} = \sqrt{(\text{kg}/\text{m}^3)/\text{m}}\cdot(\text{m}/\text{s})^{3.05}$ resolve to $\text{W}/\text{m}^2$; $k$ is an empirical dimensional constant ($1.83\times10^{-4}$ in SI). $\checkmark$
History and Development
-
The blunt-body breakthrough. In 1951 H. Julian Allen and Alfred Eggers at NACA showed that a blunt body, not a sharp one, minimises re-entry heating by throwing most of the entry energy into the air rather than the vehicle. The $1/\sqrt{R_n}$ term is the quantitative fingerprint of their insight, and it made crewed re-entry possible.
-
Fay–Riddell and Sutton–Graves. Fay and Riddell (1958) provided the rigorous stagnation-point boundary-layer theory; Sutton and Graves (1971) turned it into the compact engineering correlation used here, still a standard first estimate in entry-vehicle design.
-
Thermal protection systems. The heat flux (integrated over the trajectory to give total heat load) sizes the TPS: ablative heatshields (Apollo, Orion, Mars landers) that char and erode to carry heat away; reusable ceramic tiles and blankets (Space Shuttle); and reinforced carbon–carbon leading edges. Peak flux sets the material; total load sets the thickness.
Related Concepts: ISA Density, Density Altitude, Mach Number, Speed Of Sound, Thermal Stress, Combustion Chamber Temperature
Notes: Registry calculator reentry-heat-flux (unverified; lives in the structures module but is an
aerothermal/reentry calc). Sutton–Graves stagnation-point convective heating — correct as shipped; the
constant $1.83\times10^{-4}$ is Earth-air specific. Heating $\sim V^3$ (velocity-dominated); $1/\sqrt{R_n}$ is the
blunt-body effect (blunter = cooler). A radiative term adds at the highest entry speeds. All defaults $1.0$.