Concept

Turbulence — where it appears

The state a flow takes when the laminar solutions stop being the ones it chooses, containing motions at every scale from the flow's own size down to the millimetre. Its equations are known, deterministic and unsolvable, so it is described statistically.

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

The laminar line does not end; the flow leaves it. Friction factor against Reynolds number in a pipe. The laminar law f = 64/Re is exact and is drawn continuing past the transitional Reynolds number, faintly, because it remains a solution there — the flow simply stops taking it. The turbulent branch is Blasius' correlation and begins where experiments find transition, not where any calculation puts it.

The solutions stop being chosen

Hagen and Poiseuille's pipe profile is an exact solution of the Navier–Stokes equations at every Reynolds number, and it is linearly stable at every Reynolds number. Something else happens at 2300 anyway, and it is not that the solution stopped being one.

turbulence · Transition

The chart, with one exact line on it. The friction factor of a pipe against Reynolds number, for five relative roughnesses. Every curve here except one is Colebrook's correlation, solved by iteration rather than read off a chart. The exception is the short straight line at the left: f = 64/Re is the laminar solution and it is exact. The curves flatten to the right because once the roughness pokes out of the viscous layer the Reynolds number has nothing left to change.

The roughness a wall cannot feel

A rough pipe and a polished one carry the same flow for the same pressure over three decades of Reynolds number, and then suddenly they do not. What changed is not the pipe. It is the thickness of the film of fluid at the wall, which is the only part of the flow that can see the roughness at all.

applied · Internal flow

Ten unknowns, four equations. What is left after the Navier–Stokes equations are averaged. The mean velocities and mean pressure were there before; the six Reynolds stresses are new, and they arrived from the one term that does not average away. Nothing in the count is an approximation — the averaged equations are exact — and that is what makes the gap uncomfortable.

What averaging costs

Split the velocity into a mean and a fluctuation, average the equations, and the result is exact. It is also short of six equations, because the one nonlinear term does not average away and leaves six new unknowns behind that nothing determines.

turbulence · Closure

The ladder that never closes

Being six equations short is a problem with an obvious remedy — derive six more. The remedy works, produces an exact equation for the Reynolds stress, and leaves ten new unknowns behind. The gap does not narrow at any level, and the counting says why.

turbulence · Closure

The inertial range, and the slope read back off it. The model energy spectrum at Re = 1e+6, with production rolling off below the integral scale and dissipation cutting it off above the Kolmogorov scale. The straight middle is the inertial range, and the number printed beside it is the slope least-squares fitted to the drawn points over the middle of that range — not the −5/3 that went in.

Where the energy goes

Energy enters a turbulent flow at the largest scale and leaves it at the smallest, and in between there is nothing for it to depend on but the rate at which it is passing through. Two quantities and one dimensional argument fix the shape of the spectrum, and the exponent is −5/3.

turbulence · Cascade

Grid points against Reynolds number, and where a wing sits. The number of grid points needed to resolve every scale of a turbulent flow, which is Re^(9/4) — the cube of the ratio between the largest scale and the Kolmogorov scale. The line is the arithmetic and the marks are flows a reader can picture. An airliner's wing needs about 10¹⁷ points, and the largest calculations ever run are around 10¹².

The grid nobody can build

Resolving every scale of a turbulent flow needs Re to the nine-quarters grid points and Re cubed point-updates. An airliner's wing comes to 2·10¹⁷ points against the 10¹² of the largest calculation ever run, and no amount of patience closes a gap of five orders of magnitude.

turbulence · Cascade

80 per cent of the energy goes to the larger scale. One unit of energy is taken out of the middle wavenumber and shared between its two neighbours. Two conservation laws decide the split completely: the energy must add up, and so must the enstrophy, which weights each wavenumber by k². For (1, 2, 4) the answer is that 80.0 per cent of the energy goes up in scale and 80.0 per cent of the enstrophy goes down. There is no model of turbulence anywhere in that: it is two linear equations in two unknowns, and it is why a two-dimensional flow organises itself into large vortices while its gradients get finer.

The cascade that runs backwards

Three-dimensional turbulence carries energy from large scales to small ones and dissipates it. Take away one dimension and the term that does it vanishes identically, a second quantity becomes conserved, and two conservation laws between them force the energy to go the other way — up in scale, into ever larger vortices.

turbulence · Two-dimensional

A cross at 30.0° for every wavelength it makes. A body oscillating at ω = 0.5N in a stratified fluid, and the four beams along which its energy leaves. The angle is arccos(ω/N) from the vertical — 30.0 degrees from the horizontal here — and it is the same for every wavelength the body excites, because the dispersion relation has no length in it. The short strokes are the crests, which lie along the beams rather than across them: the phase advances perpendicular to the energy, and the two are exactly at right angles. Raise the frequency and the cross closes towards the vertical; reach ω = N and it shuts entirely, because nothing above the buoyancy frequency propagates.

The number that stops the mixing

A fluid whose density falls with height resists being stirred, and the resistance has a threshold at exactly one quarter, from an energy balance with no fluid mechanics in it. The waves such a fluid carries are stranger still — their frequency decides the direction they travel in and says nothing about their wavelength.

turbulence · Stratification

A millionth, doubling every three-quarters of a second. The separation of two trajectories started a millionth apart, on log axes against time. It grows as a straight line until it saturates at the size of the attractor, and the slope of that line is the largest Lyapunov exponent — measured here on the system by renormalising a nearby pair, not quoted. Sensitive dependence is what the straightness of the line means.

The randomness that is not in the equations

Turbulence is described in the language of statistics — means, variances, spectra, probability distributions — and none of that language appears in the equations it is a description of. The Navier–Stokes equations have no random term anywhere in them. What is random is the observer's ignorance of the initial data, and the flow's habit of amplifying it.

misconceptions · Randomness

The plateau that is the log law. y⁺ du⁺/dy⁺ across a channel, at five Reynolds numbers. Millikan's argument says this quantity must be constant wherever neither the viscous length nor the channel width may appear, and its value there is 1/κ. At Re_τ = 180 there is no flat part at all; at Re_τ = 100,000 it is flat over 2.16 decades and gives κ = 0.4120. The log law is a statement about a limit, and this is the picture of the flow approaching it.

The layer with no length in it

The logarithm in a turbulent wall profile does not come from any model of turbulence. It comes from a region where neither of the flow's two lengths is allowed to appear, and where a velocity gradient therefore has nothing to depend on but the distance to the wall. The constant in it has never been derived from anything.

turbulence · Wall law

Each cancellation costs two powers of the Mach number. Radiated power against compactness for three source clusters: a single monopole, two of opposite sign, and four on a square with alternating signs. The fitted slopes are 0.00, 2.00, 4.00 — zero, two and four in (kd), measured by integrating the far field over a sphere rather than assumed. A turbulent eddy turns over in about the time sound crosses it, so kd is of order the Mach number, and those exponents become the fourth, sixth and eighth powers of speed. A flow with no moving surfaces has no monopole and no dipole available to it, which is Lighthill's whole argument, and the eighth power is what is left.

The sound that only leaves

A flow is a catastrophically bad radiator, and the reason is that it has no monopole and no dipole available to it. What is left is the eighth power of speed — and the equation is equally happy with sound converging on a jet, which is ruled out by a condition imposed at infinity.

compressible · Aeroacoustics

Where the first grid point may go. The error in the friction a wall treatment infers, against the height of the first grid point, at Re_τ = 20,000. The velocity fed to each treatment is the closure's own, so what is plotted is the modelling of the boundary condition with nothing else in it. The log-law function is exact between y⁺ 30 and 10261 and is 64 per cent wrong at y⁺ = 1; the sublayer treatment is exact below y⁺ = 5 and hopeless above it; and the blend that most codes ship is within a few per cent everywhere and exact nowhere.

What a code says to a wall

A calculation that cannot afford to resolve the viscous sublayer has to tell the wall something else instead, and what it tells it is the law of the wall — an asymptotic result, applied at one grid point, on the assumption that the point lies in a region the calculation has not checked exists. Where it does, the answer is exact. Where it does not, the friction is out by tens of per cent, and refining the grid makes it worse.

turbulence · Wall law

One spectrum, two fields. Two one-dimensional fields built from the same amplitudes and different phases. Their energy spectra are identical to the last bit, their two-point correlations agree to 2e-15, their second-order structure functions are the same function — and no measurement of either kind can tell them apart. Everything that distinguishes them is in the phases, which is where the cascade lives.

The moment a spectrum cannot hold

A spectrum and a two-point correlation are the same object, transformed, so everything one of them records the other records too. Give a field the exact Kolmogorov amplitudes and independent phases and every second-order measurement comes out right — while the cascade, which lives in the phases, is not there at all.

turbulence · Spectrum

The one exact law, and where viscosity takes it back. The two terms of the Kármán–Howarth relation against separation, measured in Kolmogorov lengths. The four-fifths term rises linearly with r and is the whole of the law at large separation; the viscous term, 6ν dS₂/dr, falls as r^(−1/3) and takes over below 5.64η. Nothing in the calculation was told what η is: the crossing is at the same multiple of it at every viscosity tried, which is what makes the dissipation scale the lower end of the inertial range rather than a separate assumption.

The one exact result

Almost nothing in turbulence follows from the equations without a model in it. One thing does — Kolmogorov's four-fifths law, which fixes the third moment of the velocity differences at −(4/5)εr with no adjustable constant anywhere in it. The companion two-thirds law is not exact, and the 4.02 everybody quotes in it turns out to be a gamma function.

turbulence · Structure function

Two things moving by a million, and their product standing still. Viscosity, the squared velocity gradient at the dissipation scale, and their product, across six decades of Reynolds number at a fixed large-scale flow. The viscosity falls by a factor of 1e+6; the squared gradient rises by exactly the same factor, because η falls as Re^(−3/4) and u_η as Re^(−1/4); and the dissipation ν(u_η/η)² does not move at all. That is the dissipation anomaly stated as arithmetic: the limit of the dissipation as viscosity vanishes is not the value it takes when viscosity is zero.

The limit that is not the value

Dissipation is viscosity times the square of a velocity gradient, so it ought to vanish as the viscosity does. It does not. The gradient rises by exactly the factor the viscosity falls by, the product stands still, and a fluid with no viscosity at all would dissipate nothing — which is why the limit and the value are different numbers.

turbulence · Dissipation

Ballistic first, diffusive later. The mean square displacement of a parcel carried by a fluctuating velocity, from Taylor's 1921 integral with an exponential correlation of time scale 1. Below the correlation time it follows the straight-line law u²t² exactly — the parcel has not yet changed its mind — and above it the curve joins the diffusion law 2u²T·t, offset by the head start the ballistic phase gave it. The eddy diffusivity is u²T = 1.000, and it is a property of the correlation rather than of any equation of motion.

How far a parcel gets

Turbulent transport is one integral. A parcel carried by a fluctuating velocity goes in a straight line while it still remembers its own motion and performs a random walk once it has forgotten, and the crossover is the correlation time — so an eddy diffusivity is not a property a flow has, it is the limit of a measurement that has run for long enough.

turbulence · Mixing

The cheapest shape the walls allow

The parabola in a pipe is usually presented as what the equations give. It is better understood the other way round — of every profile that could carry that flow between those walls, it is the one that destroys the least energy, and the equations give it for that reason.

viscous · Minimum-dissipation

Most of the fluctuation is the switching. What a hot wire at the edge of a shear layer records, taken apart. The conditional intensity — the fluctuation inside the turbulent fluid — is flat at 0.16 by construction. The measured intensity peaks 20 per cent above it, and the excess is the third curve: the variance of a signal that keeps switching between two mean velocities, which peaks where the switching is most even and is not turbulence at all. A model calibrated against the measured curve is being fitted to an artefact of averaging.

Turbulent some of the time

At the edge of a jet or a wake a probe is inside turbulent fluid for part of the time and in perfectly smooth flow for the rest, and an ordinary time average mixes the two. Most of the fluctuation it records there is not turbulence at all — it is the switching between two states, and it peaks where the switching is most even rather than where the turbulence is strongest.

turbulence · Intermittency

What a probe in a flame sees. Half the time hot light fluid at a hundred metres a second, half the time cold heavy fluid at twenty. That is what intermittency in a jet flame looks like at a point, and it is the simplest field in which the two averages of the velocity are different numbers.

Two averages of one flow

In a flow whose density varies there are two mean velocities, they are both correct, and across a flame they differ by a factor of two. One of them is what a hot wire returns; the other is what every compressible turbulence model is written in; and the mass flux is the single product they agree on.

kinematics · Averaging

The scalar spectrum, with its two ranges. A model scalar spectrum at a Schmidt number of two thousand — dye in water. Below the Kolmogorov wavenumber it is Obukhov and Corrsin's five-thirds, inherited from the velocity; above it there is no turbulence left and the spectrum is Batchelor's minus one, which contains no velocity spectrum at all.

The scalar has its own cascade

Below the Kolmogorov scale there is no turbulence left, and a dye stirred into the flow goes on cascading anyway — on a spectrum whose exponent is minus one and whose amplitude contains no velocity spectrum at all. Resolving it costs the three-halves power of the Schmidt number, which for dye in water is a factor of ninety thousand.

turbulence · Mixing

What a mean profile cannot tell anybody

Two flows are built here with mean velocity profiles that agree to four parts in 10¹⁷. One of them carries momentum across the shear and dissipates forty per cent more energy; the other carries nothing. Everything that distinguishes them is second order in a disturbance the mean cannot see at all.

kinematics · Averaging

The production jumps and the dissipation does not. Production and dissipation against time, through a step change in the strain rate. The production follows the strain immediately — it is the strain squared times an eddy viscosity — and the dissipation moves by less than one per cent at the instant of the step, because it is set by a cascade that has not been told yet.

A dissipation that lags its production

Change the strain rate on a patch of turbulence and the production of energy follows instantly — it is the strain squared. The dissipation moves by less than one per cent, because it is the far end of a cascade that has not been told yet, and the two are out of balance by a factor of six for the next turnover.

turbulence · Dissipation

The stress a closure predicts, and the stress there is. The Reynolds-stress anisotropy through a step change in the strain rate, against what an eddy viscosity gives — which is the equilibrium value at every instant. The real stress takes about a turnover to get there, and during that turnover the closure is wrong by up to sixty per cent.

A closure with no memory at all

An eddy viscosity says the Reynolds stress is the mean strain rate times a number, now. The stress it is standing in for takes a turnover to arrive, so the closure is the zero-frequency limit of a response that has a lag in it — and the curve it is the limit of is the same shape as an aerofoil's lift deficiency.

turbulence · Closure $A particle is a low-pass filter. The fraction of a fluctuating flow's velocity that a particle follows, and the phase by which it lags, against the Stokes number. At one the particle follows 71 per cent of the motion and lags by 45 degrees — which is the corner frequency of a first-order filter, arrived at from mechanics rather than from electronics.$

A particle is a low-pass filter

A seeding particle does not report the flow; it reports the flow through its own transfer function. At a Stokes number of one it follows 71 per cent of the motion and lags it by 45 degrees, and a fifty-micron droplet at a kilohertz is following two per cent of what it is supposed to be measuring.

regimes · Particle

The heat the eddies do not carry

In a laminar layer the temperature and the velocity have different thicknesses in every fluid but one. In a turbulent pipe the eddies carry both, and the difference nearly vanishes — for air, water and oil alike. It does not vanish for a liquid metal, whose molecules conduct heat faster than the eddies can, and the boundary between the two behaviours is a Péclet number of about four hundred at every Reynolds number.

regimes · Peclet

Named alongside it

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

Dissipation Measurement Closure Correlation Boundary layer Cascade Inertial range Mixing Model validity Regime Reynolds stress Averaging

All concepts