Hand Calculations logo Hand Calculations All help pages ▾

Peak Heating Altitude⚠ unverified

Aerospace / Reentry · Compute the altitude at which peak heating occurs

Parameters

InputSymbolUnitDefaultDescription
betaβkg/m**21.0Ballistic coefficient
HHm7000.0Atmospheric scale height, in metres (m). Default is 7000
OutputSymbolUnitDescription
resulthmAltitude of peak heating, in metres (m)

The science & history

Understanding the Parameters

Derivation (Approaching a Proof)

Stagnation-point convective heating follows $q \propto \sqrt{\rho}\,V^3$ (Stagnation Heat Flux). Insert the Allen–Eggers velocity profile $V = V_0\exp[-B\,e^{-h/H}]$ with $B = \rho_0 H/(2\beta\sin\gamma)$ (Reentry Velocity) and the exponential atmosphere $\rho = \rho_0 e^{-h/H}$:

$$q \propto \rho^{1/2}V^3 = \rho_0^{1/2}e^{-h/2H}\,V_0^3\exp\!\big[-3B\,e^{-h/H}\big].$$

Take the logarithm and differentiate with respect to $h$ to find the maximum:

$$\frac{d}{dh}\Big[-\frac{h}{2H} - 3B\,e^{-h/H}\Big] = -\frac{1}{2H} + \frac{3B}{H}e^{-h/H} = 0 \;\Longrightarrow\; e^{-h/H} = \frac{1}{6B}.$$

The peak-heating density is therefore $\rho^* = \rho_0 e^{-h/H} = \rho_0/(6B) = \dfrac{\beta\sin\gamma}{3H}$, and solving $e^{-h/H} = 1/(6B)$ for the altitude:

$$h_q = H\ln(6B) = H\ln\!\frac{6\rho_0 H}{2\beta\sin\gamma} = H\ln\!\frac{3\rho_0 H}{\beta\sin\gamma}. \qquad\blacksquare$$

Compare with peak deceleration, where $q \propto \rho V^2$ gives $e^{-h/H} = 1/(2B)$ and a lower altitude $h_a = H\ln(\rho_0 H/(\beta\sin\gamma))$. Because heating carries a higher power of velocity, its peak comes earlier (higher, faster) than the deceleration peak — a key fact in sequencing heat-shield and structural loads.

Dimensional check. The argument $3\rho_0 H/(\beta\sin\gamma) = (\text{kg}/\text{m}^3)(\text{m})/(\text{kg}/\text{m}^2) = \text{(–)}$ is dimensionless (as a log argument must be), and $H\ln(\cdots)$ has units of metres. $\checkmark$

History and Development

Related Concepts: Reentry Velocity, Ballistic Coefficient, Reentry Deceleration, Stagnation Heat Flux, Reentry Heat Load, Thermal Protection Thickness, Reentry Heat Flux

Notes: Registry calculator peak-heating-altitude (unverified). Registry form $H\log(2\beta H)$ is non-standard — log argument not dimensionless, $\beta$ inverted, omits $\rho_0$/$\sin\gamma$; correct Allen–Eggers result is $h_q = H\ln[3\rho_0 H/(\beta\sin\gamma)]$. Peak heating is higher than peak deceleration (heating $\propto\sqrt\rho V^3$ vs $\rho V^2$). Flagged in Known Issues.

← Back to the workspace  ·  All help pages  ·  Getting started