Fluid Flow & Aerodynamics
why things fly, and why the usual explanation is wrong
Concept

Self-similarity — where it appears

The property of a solution that keeps its shape as it evolves, so that one function of a combined variable describes every time or station. It is what turns a partial differential equation into an ordinary one, and it is an assumption about the solution rather than a consequence of the equation.

Named by 11 essays across 3 fields — each of them below, with the objects they name alongside it.

Two decay laws from two invariants, and nothing in the equations to choose. The energy of a decaying turbulence against time, integrated from dK/dt = −A K^(3/2)/l with the large scales conserving u² l³ in one case and u² l⁵ in the other. The exponents come out at 1.1997 and 1.4282 against the closed forms 6/5 and 10/7. Which invariant holds is decided by the shape of the spectrum at the very largest scales, at the moment the stirring stops.

What decay never forgets

Stir a box of fluid and stop. The turbulence decays, at a rate with no viscosity in it — so the rate cannot come from the fluid. It comes from an invariant of the very largest scales, fixed at the moment the stirring stops, and never revisited.

turbulence · Decay Four sets of scaling exponents, all of them exact at the third moment. zeta_p against p for K41, the beta-model, the log-normal model and She–Leveque. Every one of them passes through zeta_3 = 1 exactly, because the four-fifths law is a consequence of the equations and a model that missed it would be wrong about the one thing that is known. What they disagree about is every other moment.

The exponents that stop being thirds

Kolmogorov's 1941 theory says every moment of the velocity difference scales with the same exponent, p over three, so the distribution keeps its shape at every scale. It does not. The exponents fall below the line, by more the higher the moment, and what the departure measures is a dimension.

turbulence · Intermittency Two wakes with one drag. The two initial velocity deficits: a slab, which is roughly what a bluff body leaves, and a pair of separated lobes, which is roughly what a body with a jet through the middle of it leaves. Their integrals are identical, so the two bodies have exactly the same drag.

A wake that keeps the drag and forgets the body

Two very different wakes with the same momentum deficit converge to the same profile, because the deficit is conserved and everything else diffuses away. The convergence is a power law rather than an exponential, so it takes two hundred widths for twenty per cent agreement and nine hundred for five.

turbulence · Wake Five spectra, five exponents, one linear equation. The energy of a decaying turbulence after the nonlinear term has stopped mattering, computed by integrating the exact modal solution E(k,0)exp(−2 nu k² t) at five different shapes of the spectrum at the origin. Each is a straight line on these axes and no two have the same slope: the exponent is (m+1)/2, where k^m is the spectrum's behaviour at wavenumbers smaller than any eddy. The 5/2 that is quoted as the final period's universal exponent is the m = 4 line and one of five.

Universal, and one of five

A turbulence that has run its Reynolds number down stops being turbulent, the equations go linear, and the decay picks up a new exponent. That exponent is quoted everywhere as 5/2 and as universal. It is neither: it is the same corner of the same spectrum deciding the answer a second time.

turbulence · Decay The coefficient that was a constant, against the number it is said not to depend on. The dissipation coefficient Cε = eps·l/u³ along two decays, plotted against the Taylor-scale Reynolds number they pass through. One is flat because it was put in flat; the other falls as the reciprocal of the Reynolds number, which is what is measured in the near field of a grid. Neither line is a derivation. What is exact is the relation between them, Cε = 15(ℓ/λ)/Reλ, which is a rearrangement of two definitions and holds along both curves to 5·10⁻¹⁶.

The constant that travels

Every decay law in the subject rests on the dissipation being some constant times u³ over a length. The constant is not one. Letting it move the way grid measurements say it moves changes the decay exponent by a quarter — and lands one of the answers five per cent from another that is entirely different physics.

turbulence · Decay The kept transport spirals into nothing as the sea deepens. The net Lagrangian transport as a vector, scaled on the Stokes transport, traced as the water depth increases from a quarter of an Ekman depth to eight, for an 8-second swell with an eddy viscosity of 0.01 m²/s. Shallow water keeps the whole transport pointing with the waves, at the right-hand end. As the sea deepens the vector shortens and swings to the right, crosses the across-wave axis near two Ekman depths, and winds into the origin, which is the open ocean's exact cancellation.

The floor that gives the drift back

In the open ocean the Coriolis force drives a current that cancels a swell's Stokes transport exactly. Over a continental shelf the sea floor holds a stress, and whatever it holds is transport the rotation does not take back. How much survives depends almost only on the depth in Ekman depths; which way it points depends on the wave.

kinematics · Stokes drift Three pairs released at different distances end on one line. The mean square separation of two parcels against time, both in Kolmogorov units, for pairs released one, thirty and a thousand Kolmogorov lengths apart in turbulence whose integral scale is a million of them. Each starts flat — the pair has not yet moved — and each joins the same line, g ε t³ with g = 0.5, after which nothing about where it started survives. Beyond the integral scale all three become ordinary diffusion, growing as t.

The further apart, the faster they part

Release two specks of smoke close together in turbulence and the distance between them does not diffuse the way a single speck wanders. It grows faster the larger it already is, because larger separations are pulled apart by larger eddies — as the cube of the time, forgetting where it started. The law is ninety years old, its constant is still argued about, and it needs a very large flow to be seen at all.

turbulence · Mixing Inside a strong blast: a shell of gas around a hot, nearly empty core. Left, the density, velocity and pressure inside a spherical strong blast wave at γ = 1.4, each as a fraction of its value just behind the shock, against the distance from the centre as a fraction of the shock radius. The density falls to nothing well inside the shock while the pressure levels off at 0.37 of its post-shock value. Right, the temperature, which is the pressure over the density and so rises without limit towards the centre.

The fireball is hollow

Inside a strong blast wave there is almost nothing. Half the air the shock has swept up lies in the outer four per cent of its radius; at half the radius the density is a hundredth of the air's; the centre is empty and, formally, infinitely hot. Integrating Sedov's equations to see this also shows that a published table of his constants was wrong at three of its four entries.

compressible · Blast wave A wrinkled flame settles into arcs meeting at a cusp. The steady front of a flame in a periodic domain 5, 10 and 20 neutral wavelengths wide, from the exact pole solution of the Michelson–Sivashinsky equation, with the burnt gas below and the flame advancing upwards; each is drawn across one period, scaled to the same width. Every one is a single smooth arc bulging into the fresh gas, meeting its neighbour in a sharp cusp pointing back into the burnt gas, and in these units the three arcs nearly coincide: only the cusp sharpens as the domain widens. In physical units the arc's depth grows in proportion to the domain's width, so the three flames are the same shape at three sizes.

A wrinkled flame has one cusp and a speed limit

The linear theory of a flame says every long wrinkle grows and none is favoured. The weakly nonlinear theory — the Michelson–Sivashinsky equation — says where the growth goes: small wrinkles merge, the front settles into smooth arcs bulging into the fresh gas and meeting in sharp cusps, and in a domain of any width it ends with a single arc and a single cusp. That front is an exact solution made of poles in the complex plane, and its speed is a closed form that rises in steps as the domain admits more poles and then stops: beyond about five neutral wavelengths a wider flame is no faster, because it is the same shape at a larger size.

kinematics · Dilatation A bow shock that grows as a line explosion does. The radius of the bow shock around a hemisphere-nosed cylinder, in body diameters, against distance behind the nose: the blast-wave analogy, R/d = 0.795 C^¼ (x/d)^½ with C, the drag coefficient, 0.919, which contains no Mach number, and Billig's correlation of measured bow shocks on spheres at Mach 5, 10 and 20, anchored at the nose and continued as a hyperbola to the Mach cone. Both grow as the square root of the distance; at Mach 20 the analogy's shock is a steady 0.71 of Billig's.

A hypersonic body leaves a line explosion behind it

A blunt body at hypersonic speed does work on the air at a rate equal to its drag, and each slice of air it passes through is struck once and left to expand. Seen from the ground, that is a line explosion, and Sedov's cylindrical blast wave gives the bow shock's width and the pressure on the afterbody without a Mach number in either. The analogy's classical constants come straight out of the blast solution. So do its limits: it holds only while its own shock stays strong, over a length that grows as the square of the flight Mach number, and it puts the body inside a core hotter than anything the flow can reach.

compressible · Blast wave An adverse gradient hollows the layer from the outside in. The velocity profile of an equilibrium layer, as a fraction of the edge velocity against y/δ, at Clauser's β = 0, 2, 10 and 50, with the wake strength from Das's fit. The inner part keeps the wall law; the outer part's wake grows — Π = 0.476, 1.59, 4.67 and 15.2 — and the profile sags until most of the layer is moving slowly over a thin fast wall region.

An adverse gradient grows the wake, not the log law

Push a turbulent boundary layer up a rising pressure and its profile changes from the outside in. The wall region keeps its law; the wake grows, the layer hollows, the shape factor climbs and the friction falls. Clauser found the layers that stay similar while this happens, and their arithmetic has a surprise in it: there is a steepest pressure rise any such layer can climb, close to a free stream falling as the inverse fourth root of the distance, and the friction approaches zero only as the layer becomes all wake.

turbulence · Wall law

Named alongside it

The objects these essays reach for when they reach for this one.

Model limitDecayEddy viscosityIntegral scaleInvariantMeasurementPower lawReynolds numberBlast waveClosureConserved quantityDiffusion

All concepts