Compressible flow

What a shock costs

All the heat survives a shock and none of it is lost. What is lost is the ability to turn that heat back into speed — 27.9 per cent of it at Mach 2 and 93.8 per cent at Mach 5 — and every supersonic intake ever built is a scheme for paying less.

Worth reading first: The only law that forbids it.

The first thing to get straight about a shock is what it does not do. It does not destroy energy. Every joule that went in comes out, and the stagnation temperature — which is the energy, expressed as a temperature — is identical on both sides.

What it destroys is the flow’s ability to give that energy back as ordered motion.

Everything a normal shock does, against the Mach number in front of it. Four quantities across a normal shock, each scaled to fit one axis. Pressure and density rise without limit and without bound as the Mach number grows; the Mach number behind falls towards a floor it never passes; and the total pressure — the flow's ability to be turned back into speed — collapses. That last curve is why a supersonic intake is designed around avoiding a single strong shock.
Fig. 1 Four quantities across a normal shock, against the Mach number in front. Pressure and density rise without bound, the Mach number behind falls to a floor at 0.378 and stops, and the total pressure collapses — 72.1 per cent surviving at Mach 2, 32.8 at Mach 3, 6.2 at Mach 5. That falling curve is the cost, and it is the only one of the four that an engine designer loses sleep over.

Two totals, and only one of them is a conserved quantity

The energy equation gives the stagnation enthalpy, and hence the stagnation temperature, as constant along a streamline through anything adiabatic — smooth or discontinuous, reversible or not.

The stagnation pressure is a different animal. It is the pressure the flow would reach if brought to rest isentropically, and that word is doing all the work. A flow that has produced entropy cannot be brought back to its original pressure, however carefully it is decelerated, because the process that would take it there is not available.

assertTotalsAcrossShock checks both halves. It recomputes the stagnation temperature independently on each side of the jump from that side’s own static state and Mach number, and requires them equal to within 101210^{-12}; and it requires the total pressure ratio to lie strictly between zero and one. The rejection test moves the temperature ratio by one part in a thousand, an amount no figure would render visibly, and requires the refusal.

The currency

Total pressure is the working currency of compressible aerodynamics, and it is worth understanding why rather than accepting it as convention.

Everything a gas can be asked to do — accelerate through a nozzle, push a turbine, produce thrust — depends on the pressure ratio available between where it is and where it is going. A stream with a high total pressure and a stream with a low one at the same total temperature carry the same energy and are not equally useful, and the difference is precisely what an engine sells.

The relation to entropy is exact and short:

p02p01=eΔs/R\frac{p_{02}}{p_{01}} = e^{-\Delta s / R}

so total-pressure ratio is the entropy rise, exponentiated. Two ways of saying one thing, and the pressure form is preferred by engineers because a pitot tube measures it.

The numbers, and where they become intolerable

A single normal shock at the flight Mach number gives these total-pressure recoveries:

Mach recovery
1.2 0.9928
1.5 0.9298
2.0 0.7209
2.5 0.4990
3.0 0.3283
5.0 0.0617

Below about Mach 1.5 a normal shock is nearly free and nobody bothers avoiding it. At Mach 2 losing 28 per cent of the total pressure is a serious hit but survivable. At Mach 3 a normal-shock intake throws away two thirds of what the flight condition offered, and no engine downstream can recover it.

At Mach 5 the recovery is six per cent, which means an engine fed by a single normal shock has essentially nothing to work with. That is a hard bound on the architecture and it is why hypersonic propulsion looks the way it does — and, past a point, why it stops decelerating the flow at all.

The cube law, which is the way out

The escape is the fact that entropy production goes as the cube of shock strength.

The second law is the only thing that forbids the other half of this curve. Entropy change across a normal shock, against the Mach number in front of it. The solid branch is the compression shock that exists. The dashed branch below Mach one is the expansion shock, and it satisfies mass, momentum and energy exactly — the residuals are zero to machine precision. It is refused by the second law alone, the one statement in the problem that no conservation residual can show.
Fig. 2 The currency itself. Entropy change across the shock against the Mach number in front of it: the solid branch is what happens and the dashed one below Mach one is the expansion shock, which satisfies mass, momentum and energy exactly and is refused by this curve alone. Everything this essay prices is read off it.

The exponent is fitted rather than quoted. weakShockEntropyExponent computes the entropy rise at four small turning angles at a given Mach number, fits a straight line through log Δs\Delta s against log θ\theta, and returns the slope; assertWeakShockEntropyIsCubic requires it to be 3 within 0.06. At Mach 2 the fitted value is 2.985, the departure from exactly 3 being the higher-order terms that a finite fitting range still contains.

The consequence, in numbers: at Mach 2, a fifteen-degree compression taken as a single oblique shock keeps 95.9 per cent of the total pressure. The same fifteen degrees taken as four shocks of 3.75° each keeps 99.4 per cent. Three quarters of the loss has been avoided by changing nothing but the number of steps.

Every attached shock there is: deflection, shock angle, Mach number. For each Mach number, the shock angle that produces a given deflection. Each curve is double-valued — a weak solution low down and a strong one high up — and each has a maximum. To the right of the dashed ridge through those maxima there is no attached shock at any angle, and the flow answers with a curved bow shock standing off the body.
Fig. 3 The shock angles available at four Mach numbers. Every point on a curve is a shock, and every one of them costs something — but a shock low on a curve turns the flow a little and leans far forward, which makes it weak, while one high up turns it hard and stands nearly normal. Taking a large turn as several small ones means staying near the bottom of these curves, where the cube law makes the entropy negligible.

Which is why intakes look like that

Every supersonic intake is a device for taking the compression in as many small bites as possible.

A pitot intake — a plain hole facing forward — has one normal shock and takes the whole loss. It is used up to about Mach 1.5, where the loss is under seven per cent and the simplicity is worth it.

An external-compression intake puts a wedge or a cone ahead of the opening, so the flow crosses one or more weak oblique shocks before meeting the final normal shock at a much reduced Mach number. Two oblique shocks plus a normal one at Mach 2 recover around 92 per cent against the single normal shock’s 72.

An isentropic intake takes the compression through a continuously curved surface, which is a compression fan of infinitely many infinitesimal waves and produces, in principle, no entropy at all. In practice those waves converge and coalesce into a shock a short distance out, so the ideal is approached and not reached.

The progression is a direct reading of the cube law, and every step of it costs geometry, weight and off-design behaviour. That is engineering: the second law hands over a graded penalty, and the design is a purchase of how much of it to pay.

The efficiency that is not a percentage

A recurring confusion is worth heading off, because it makes the numbers in this essay sound worse than they are and occasionally better.

Total-pressure recovery is not an efficiency in the sense of “what fraction of the energy is usefully employed”. No energy has been lost at all — the stagnation temperature is unchanged — so an intake with 72 per cent recovery has not thrown away 28 per cent of the incoming energy.

What it has thrown away is 28 per cent of the pressure ratio it could have delivered, and the consequence for the engine’s output is neither 28 per cent nor nothing. It depends on the cycle: an engine with a high pressure ratio of its own is relatively insensitive to losing some of the intake’s contribution, while a ramjet, which has no compressor and lives entirely on ram compression, loses almost proportionally.

So the same shock loss is a minor nuisance in one machine and fatal in another, which is why the quantity tabulated is the recovery rather than any efficiency: recovery is a property of the flow alone, and what it is worth is a property of the machine attached to it.

Where the lost availability actually ends up

Saying that the ability to give the energy back has been destroyed is correct and slightly abstract, and the concrete version is worth having because it is measurable.

The entropy rise means the gas can no longer reach its original pressure by any adiabatic deceleration. Read the same statement the other way and it says the flow, brought to any given static pressure downstream, arrives there hotter than it would have without the shock. The energy is all present and some of it is now in a form that is thermally degraded — at a higher temperature and a lower pressure than the same energy occupied before.

That has a direct consequence in a combustion system, and it is the reason intake losses are counted so carefully. A combustor is limited by the temperature its turbine can survive. Air arriving hotter because it has crossed a shock leaves less room between its own temperature and that limit, so less fuel can be burned in it — and the engine’s output falls for a reason that is not a pressure ratio at all. The same shock loss appears twice in the cycle, once as a pressure the nozzle does not have and once as a temperature margin the combustor does not have.

The overall sensitivity is the number a propulsion engineer carries, and it is unforgiving. For a supersonic turbojet, one per cent of intake total-pressure recovery is worth something like one to one and a half per cent of net thrust; for a ramjet, which has no compressor and lives entirely on ram compression, it is worse and approaches proportional. Against the table above — twenty-eight per cent lost to a single normal shock at Mach 2 — that is the whole argument for the geometry of a supersonic intake in one line.

Why total pressure is the currency and not merely a quantity

There is a second reason the loss is tracked in this form rather than as an entropy or an efficiency, and it is the practical one. Losses in total pressure add.

A duct’s shock system, the friction on its walls, the mixing of a boundary layer with the core, the wake of a strut, the turning loss in a bend: each of them produces entropy by its own mechanism, and each appears as a decrement of p0p_0. They can therefore be measured separately, tabulated as loss coefficients, and summed — which is exactly what an internal-aerodynamics calculation does, and why every compressor cascade, diffuser and duct in the literature is characterised by a number of the form

ω=p01p02p01p1,\omega = \frac{p_{01} - p_{02}}{p_{01} - p_1},

a total-pressure loss normalised by an inlet dynamic pressure.

Nothing else in the problem has that property. Static pressures do not add — they are related by the geometry. Efficiencies multiply and are defined against different ideal processes in different components. Entropies do add, and are inconvenient to measure. Total pressure is the one quantity that is both directly measurable with a probe and additive across mechanisms, and that combination is why an entire branch of engineering is written in it.

The other cost, which is drag

A shock takes total pressure, and it also pushes.

The pressure rise across a shock acts on whatever is behind it, and if that is a body, the result is a force along the flow direction: wave drag. It is a drag that exists in a fluid with no viscosity whatever, which makes it a genuine exception to d’Alembert’s paradox and one of the few places on this site where the ideal theory’s famous result fails on its own terms.

Everything a normal shock does, against the Mach number in front of it. Four quantities across a normal shock, each scaled to fit one axis. Pressure and density rise without limit and without bound as the Mach number grows; the Mach number behind falls towards a floor it never passes; and the total pressure — the flow's ability to be turned back into speed — collapses. That last curve is why a supersonic intake is designed around avoiding a single strong shock.
Fig. 4 The same four quantities out to Mach 12, where the bill stops being payable. Total pressure — the flow’s ability to be turned back into speed — collapses, while the static pressure rise it is being bought with grows without limit. That is the whole reason a hypersonic intake is not a normal shock.

The energy accounting behind the two costs is the same accounting. Wave drag times flight speed is the rate at which the vehicle does work on the air, and that work reappears as the energy carried away in the wave system — which is exactly the availability the total-pressure loss measures. The drag and the entropy are two faces of one transaction, and the whole ladder that follows from it is a field of its own.

What it costs a nozzle, which is the same sum backwards

The accounting works in the other direction too, and in a nozzle the loss is easier to see because it shows up directly as thrust that was not produced.

One nozzle, five back pressures, five different flows. The duct above, and the static pressure along it below, computed station by station from the local area. Where a shock stands inside the divergent section its position was solved for rather than placed: the shock spends total pressure, which sets the subsonic Mach number at the exit, which has to match the imposed back pressure.
Fig. 5 One nozzle at five back pressures. The two middle cases carry a shock inside the divergent section, and each of them is throwing away total pressure at exactly the rate this essay’s first figure gives — which then limits the exit velocity, because exit velocity is set by the pressure ratio the flow still has available to it.

A nozzle whose shock stands at an area ratio of 2.43 meets flow at Mach 2.41 and loses 46.6 per cent of its total pressure, and the exit velocity it then achieves is a fraction of the design value. The propellant burned is unchanged; the thrust is not. Everything about how a nozzle’s shock position is set by its back pressure is therefore also a statement about how much of the chamber’s energy reaches the exhaust.

That is why an overexpanded rocket engine at sea level is worse than a smaller nozzle would have been, and why staging exists: it is cheaper to change the nozzle than to fly the wrong one for the whole ascent.

What the shock does to the boundary layer

The pressure rise a shock imposes is, from the point of view of the thin layer on a surface underneath it, an adverse pressure gradient — a violently adverse one, applied over a distance of almost nothing.

A boundary layer can only climb so much of a pressure rise before separating, and a shock delivers more of it faster than any smooth deceleration would. So a shock standing on a surface usually separates the layer beneath it, and the separation feeds back and changes the shock, which changes the separation.

That interaction is where supersonic and transonic aircraft actually get into difficulty: buffet, control reversal, and intake instability are all shock–boundary-layer interaction rather than shock loss. This site’s solver does not compute it. What it can say honestly is the size of the pressure rise being applied — 4.5 times at Mach 2 — against the pressure rise a laminar layer can survive at all, which is a fraction of that.

A normal shock at Mach 1.40, and what crosses it unchanged. The state in front of the shock and the state behind it. Every ratio was computed from the standard jump relations and then substituted back into mass, momentum and energy, which is an independent route — a mistyped exponent in the total-pressure expression cannot survive a momentum balance it never appeared in. The residuals are printed below because a check nobody can see is a check nobody can audit.
Fig. 6 A weak shock, at the sort of Mach number that appears in a supersonic pocket over a transonic wing. Even here the static pressure rises by a factor of 2.12 across a distance of essentially nothing, and the total-pressure loss is only 3.2 per cent. The shock is nearly free thermodynamically and is still an enormous adverse gradient for the layer underneath it — which is why the transonic trouble is a boundary-layer problem rather than a loss problem.

What the solver computes, and how it is checked

Three assertions carry this essay.

assertTotalsAcrossShock proves that the stagnation temperature does not move and the stagnation pressure does, by recomputing each on both sides independently. assertShockRaisesEntropy proves that the loss is in the direction the second law requires. assertWeakShockEntropyIsCubic fits the exponent rather than assuming it, which is the difference between a claim and a measurement.

The third is the one worth dwelling on, because the whole intake argument rests on it and it would be easy to state without checking. A fitted exponent of 2.985 is evidence; “entropy goes as the cube of shock strength, as is well known” is not. The rejection test hands the check an exponent of 2.0, which is what a plausible algebra slip would produce and which would make the multi-shock intake a far weaker idea than it is, and requires the refusal.

Where the model stops

Adiabatic, and no work done. Total temperature is conserved because nothing adds heat and nothing does shaft work. In a real intake there is heat transfer to the structure, and downstream of it there is a compressor doing exactly the work this analysis excludes.

Uniform flow on both sides. The recovery figures above treat the intake shock system as a set of uniform states. Real intakes have non-uniform exit flow, and the distortion it presents to the compressor is a design constraint at least as serious as the mean recovery.

No viscous loss. Everything counted here is shock loss. Skin friction in the duct and separation losses add to it and are frequently comparable, so an intake’s measured recovery is always below the shock-only prediction.

Steady operation only. An intake that unstarts — losing its swallowed shock, which then pops out in front of the inlet — does so violently and unsteadily, and none of this analysis describes the event.

What the picture cannot show

The recovery curve is a clean falling line, and it hides the thing that makes total pressure worth tracking: that it is irrecoverable.

Nothing in a plot of p02/p01p_{02}/p_{01} against M1M_1 says that the loss is permanent. A reader could reasonably take it for a resistance that a sufficiently good diffuser might undo, in the way a pressure drop through a bend can be partly recovered downstream. It cannot. The lost total pressure is not stored anywhere and there is no process that returns it.

That permanence is a thermodynamic statement, and — like everything else on the previous rung — it has no picture at all. The curve shows the size of the loss. Its irreversibility is carried entirely by the caption.

Who found it, and when

The total-pressure accounting is Prandtl’s school, in the 1900s and 1910s, arriving from thermodynamics rather than from hydrodynamics — which is why the vocabulary is a thermodynamicist’s.

The multi-shock intake is Oswatitsch, in 1944, in a report that is one of the cleanest pieces of applied second-law reasoning in aerodynamics: it works out the optimum distribution of turning among a fixed number of oblique shocks, and finds that the shocks should be of equal strength. That result follows directly from the cube law — equal cubes minimise the sum for a fixed total — and it is still how intakes are laid out.

The delay between the theory and the hardware was short here, for the usual reason: the aircraft that needed it were being built at the time.

Where the ladder goes next

Every shock in this ladder has been normal to the flow, and a normal shock is the most expensive kind there is. A shock that leans takes the same flow through the same conservation laws with only its normal component changed, which makes it weaker, cheaper, and able to leave the flow supersonic — and that is the next ladder.

What links here

Computed from the collection rather than written here: the essays that point at this one.

Reads more easily once this is understood

Essays that name this one as worth reading first.

Shares its objects with

Essays naming at least two of the same things, that neither author linked.

Named objects

A dashed tag is an object no other essay names yet.

AvailabilityEfficiencyEntropyIntakeMach numberThe second law of thermodynamicsShock waveStagnation temperatureTotal pressure