Concept

Linear stability — where it appears

The question of whether an infinitesimal disturbance to a flow grows or decays, answered by an eigenvalue problem. It says what happens at the very start of a departure from the laminar state, and it is silent about everything after that.

Named by 30 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
Two drag laws, one derived and one fitted. Flat-plate drag coefficient against Reynolds number, laminar and turbulent, on log axes. The laminar curve is Blasius' similarity solution, solved by shooting; the turbulent one is the 1/7-power correlation, which was fitted to experiment and is drawn dashed-in-kind to keep the difference visible. Their slopes differ — −1/2 against −1/5 — so the gap between them widens rather than staying put.

The number that is not a number

Transition Reynolds numbers are quoted to three figures and vary by two decades. That is not sloppiness in the measurement — it is the honest report of a quantity that depends on the laboratory as much as on the fluid, and knowing which part is which decides what may be designed on it.

turbulence · Transition
A shear layer, and the point of inflection in it. The velocity profile U = tanh y across a layer of finite thickness, with the inflection point located by searching for a sign change in the second derivative rather than by reading it off the algebra. Rayleigh's theorem says an inviscid parallel flow can only be unstable if such a point exists — a necessary condition, not a sufficient one.

A layer with a kink in it

Rayleigh proved in 1880 that an inviscid shear flow cannot be unstable unless its velocity profile bends the other way somewhere. It is one line of algebra, it is necessary and not sufficient, and it ties instability to separation through the sign of a single derivative at the wall.

turbulence · Instability
Every wavelength grows, and the shortest grows fastest. Growth rate against wavenumber for a vortex sheet of zero thickness and for a layer of thickness δ. The sheet's curve rises without limit, which is the model announcing that it has no shortest scale in it; the layer's turns over and dies at kδ = 1, the neutral wavenumber computed from Rayleigh's equation rather than fitted.

Every wavelength at once

A vortex sheet of zero thickness is unstable at every wavelength, and the shorter the wavelength the faster it grows. The answer has no smallest scale in it, which is not a fact about fluids — it is the model reporting that it left something out.

turbulence · Instability
Every disturbance has its own threshold; one of them is lowest. The Rayleigh number at which a disturbance of horizontal wavenumber a becomes neutral, Ra = (π² + a²)³/a². Every wavenumber has a threshold and the layer goes unstable at the lowest of them, which a golden-section search on this curve puts at a = 2.221441 and Ra = 657.5114 — the exact π/√2 and 27π⁴/4 to fourteen digits. Below the curve the layer conducts and nothing moves.

A threshold with a closed form

A layer of fluid heated from below sits still until buoyancy overcomes both diffusions at once, and then it convects. Unlike every other threshold in this field, that one is an eigenvalue with an exact answer, and the answer is 27π⁴/4.

turbulence · Convection
The street, as two rows of point vortices. The exact velocity field of a staggered double row of point vortices at the stable spacing ratio, with the row spacing and the circulation of one core printed from a line integral of the field rather than from the number that built it. This is a model of a wake. It contains no body, no viscosity and no mechanism that would shed anything, and the viscous stepper used here does not produce a street at any Reynolds number.

The street a steady calculation cannot draw

The alternating wake behind a cylinder is the most photographed structure in fluid mechanics, and a steady solution does not produce one. What can honestly be drawn instead is a model of it — and the model settles one thing exactly, which is the spacing.

turbulence · Wake
Stable in every mode, and 100 times larger first. The energy of the worst-case disturbance against time, on a logarithmic scale, at four Reynolds numbers and at the one the slider selects. Time is in units of Re, which is what makes the four curves the same shape; what changes with Reynolds number is the height, and it changes as the square. Every eigenvalue of this operator is negative throughout, so nothing that grows here is an instability in the sense a stability analysis reports.

Every mode decays and it grows anyway

A stability analysis asks whether any mode of a flow grows, and for pipe flow the answer is no, at every Reynolds number, which the pipe disagrees with. The missing ingredient is that the modes are not perpendicular — a disturbance made of two nearly parallel decaying pieces can grow by a factor of Re²/16 before it dies.

turbulence · Transition
The jet has to be fed from the sides. The velocity field of the plane jet with streamlines integrated through it. The seven central streamlines run down the jet and spread; the ten started at the top and bottom edges bend inwards and join it, which is entrainment and is a consequence of the solution rather than an addition to it. The dashed lines are the half-speed edges, widening as x^{2/3}. The transverse velocity far from the axis is 5.70e-3 m/s at this station, inward on both sides — a jet is a sink as seen from a distance, which is why two parallel jets pull together.

What a jet keeps, and what it collects

A jet leaving a nozzle into still fluid has no boundary anywhere and one conserved quantity. Its momentum flux is exactly the same at every station downstream; its mass flux is not conserved at all and grows without limit, because a jet is a machine for acquiring fluid it did not start with.

viscous · Free shear
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
Three thresholds, and none of them is one. The neutral curves for a layer heated from below, computed by taking the smallest eigenvalue of the marginal-stability operator at each horizontal wavenumber. Every curve diverges at both ends — a cell wider than the layer has to carry heat sideways for ever, a narrower one loses it to conduction — so each has a minimum, and that minimum is the critical Rayleigh number. Two free surfaces give 657.5, one rigid and one free 1100.7, two rigid walls 1707.8. Nothing but the boundary condition differs, and it carries a factor of 2.6.

The threshold the walls decide

A layer heated from below convects at a Rayleigh number of 1707.762, and nothing whatever happens at one. The free–free case has a closed form and the two that do not differ from it by a factor of 2.6 — produced by nothing but what the top and bottom surfaces are permitted to do.

turbulence · Convection
The marginal Taylor number against the axial wavenumber. The smallest Taylor number at which a disturbance of a given axial wavenumber is neutral. Its minimum is 1707.757 at a wavenumber of 3.1158, which are the critical Rayleigh number and critical wavenumber of a layer of fluid heated between two rigid walls — the same numbers, because in the narrow-gap limit the two problems are the same sixth-order eigenvalue problem. This curve is computed by that essay's own solver.

A transition that needs a second number

Fluid between rotating cylinders goes unstable at a Taylor number of 1707.762 — which is the same number, to every digit, as a layer of fluid heated between two rigid walls. It is not an analogy. And two things the number cannot carry decide whether the transition happens at all and what it looks like when it does.

regimes · Taylor
Every unstable mode of every profile, inside one circle. The complex phase speeds of the unstable modes, scaled so that each profile's own semicircle is the unit one. Howard's theorem says every one of them must lie inside — the centre and the radius are the mean and half-range of the velocity profile and nothing else — and every one of them does, with the closest approach at 0.915 of the radius.

Every unstable wave is inside one circle

Before solving anything, you know where the answer is. Howard's theorem says the complex phase speed of any growing disturbance in a shear flow lies inside a circle fixed by the fastest and slowest parts of the profile — and by nothing else about it at all.

turbulence · Instability
Instability up to a quarter, and none past it. The fastest growth rate of a stratified shear layer against its Richardson number, on a profile whose gradient Richardson number is the same at every height. It falls smoothly towards zero and reaches it at a quarter: at Ri = 0.2499 the fastest mode still grows at 0.00106, and at 0.26 the solver finds no unstable mode at all.

Sufficient, and not necessary

A stratified shear layer whose Richardson number exceeds a quarter everywhere cannot go unstable. That is a theorem with an exact number in it. What it does not say — and what it is constantly read as saying — is that a layer below a quarter will.

turbulence · Stratification
Below Thoma's 6.07 m² the governed tank's swing grows; above it, it dies. The tank level after the turbine's power demand drops by two per cent, with a governor holding the power constant, for tanks of 0.7 and 1.3 times Thoma's area of 6.07 m² — a tank 2.78 m across. The smaller tank's swing grows by a factor of 1.47 every 74 s cycle and has reached −12.20 m by 427 s; the larger one's keeps 0.75 of itself every 100 s and is barely visible. Carried on, the smaller tank's run is refused at 794 s, where the head at the turbine has fallen below a quarter of its design value and the governor would be asking for a flow no turbine passes. The instability has nothing to do with the tank's height: it is the governor drawing more water as the level falls, which feeds the swing, against the tunnel's friction, which is the only thing damping it.

The better tunnel needs the bigger tank

A turbine governed to hold its power opens further when the head at it falls, and draws the tank down harder. That makes it a negative resistance, the tunnel's friction is the only thing damping the swing against it, and so the smallest stable tank grows as the friction shrinks — 2.78 metres across for five metres of friction, 6.09 for one.

applied · Water hammer
The damping crosses zero once, and the crossing is the boundary. The least damping ratio of the four aeroelastic modes against speed. It falls through zero at 80.843 metres a second, and at that speed the crossing root's real part is four parts in 10¹⁸ — which is what an algebraic condition on a quartic with real coefficients looks like when it is solved numerically.

The speed where the damping is exactly zero

Flutter is an algebraic condition on a quartic: one root crosses the imaginary axis, at one speed, exactly. And the quantity that locates it is so nearly flat there that the standard way of finding it from flight test overshoots by a quarter.

circulation · Flutter
The same wavelength drawn as rolls and as hexagons. Plan views of a convecting layer with one critical wavelength, drawn from the amplitude equations' two stable states. On the left, rolls: a single set of parallel bands, rising fluid along one set of lines and sinking along the next. On the right, hexagons: three sets of rolls at 120° to each other with equal amplitudes, whose sum has its maxima on a triangular lattice with spacing 2/√3 of the wavelength, each maximum at the centre of a hexagonal cell. At ε = 0.0250, inside the window, both are stable: rolls with amplitude 0.158 and hexagons with 0.081 in each of their three rolls. A hexagon is not a different kind of cell; it is three roll patterns that the quadratic term lets reinforce one another.

Hexagons remember how the heat was turned up

A layer heated from below convects in rolls, unless its top and bottom are not mirror images of each other. Then three sets of rolls at 120° can feed one another through a term the symmetry used to forbid, hexagonal cells appear before the layer is formally unstable, and there is a range of heating in which rolls and hexagons are both stable — so the pattern a layer shows depends on whether the heat was turned up or down to get there.

turbulence · Convection
A band of growing waves that opens at 5772 and narrows as the viscosity goes. The wavenumbers at which a two-dimensional wave on plane Poiseuille flow neither grows nor decays, against the Reynolds number on a logarithmic axis. Inside the tongue waves grow; outside they decay. The tongue's tip is the critical point. Both edges slope downward and towards each other in wavenumber as the Reynolds number rises, so the band of unstable waves shrinks towards long waves — the direction in which the inviscid problem, which has no growing wave at all, is reached.

The profile Rayleigh cleared and viscosity did not

Flow between two plates has no inflection point, so without viscosity no wave on it can grow. With viscosity one does, above a Reynolds number of 5772. Taking the viscosity away again slows that wave and narrows the band it grows in, because the stress that feeds it is made by viscosity in the first place.

turbulence · Instability
Rolls turning side by side, with the fastest downwind water where they sink. The fastest-growing mode at a Langmuir number of 0.13, looking downwind, over two roll spacings of 2.89 decay depths and 4 decay depths down. The closed curves are streamlines of the overturning; the dashed curves are contours of the downwind velocity the rolls carry, positive under the lines where the water sinks. At the surface the cross-wind flow converges onto those lines, which is where floating foam and weed collect as windrows. The amplitude is arbitrary, as in any linear mode.

The drift that turns a current into rolls

A current carrying a Stokes drift feels a force the drift makes out of the current's own vorticity, and under a wind that force is unstable. It turns the surface layer into rolls lined up downwind, with windrows where they sink. The rolls need both the current's shear and the drift's; their growth rate sees only the product; and the split between the two decides which motion gets the energy.

kinematics · Stokes drift
Three speeds, and the Froude number at which two of them cross. The kinematic wave speed and the two dynamic wave speeds, all divided by the speed of a shallow-water wave on still water, against the Froude number. The kinematic wave is three halves of the water speed and the downstream dynamic wave is the water speed plus one, so they cross at Froude two — and beyond that crossing the news of a change in flux arrives before the news of a change in depth, which uniform flow cannot survive.

When the flux outruns the pressure

A kinematic wave is derived by throwing the momentum equation away, and the derivation never says when that is allowed. It is allowed until the kinematic wave, which travels faster than the water, overtakes the downstream dynamic wave as well — at which point uniform flow ceases to exist and a concrete chute carries a train of surges instead of a sheet.

kinematics · Kinematic waves
Which bends of a vortex pair grow, and how fast. The square of the growth rate of a sinusoidal bend of a trailing vortex pair, in units of one e-fold per descent time, against the bend's wavenumber times the spacing, for a core of elliptic loading's size. Where the square is positive the bend grows. The symmetric mode, in which the two vortices bend as mirror images, grows over a long-wave band and fastest at a wavelength of 8.54 spacings, by 0.827 e-folds in the time the pair takes to sink one spacing. The antisymmetric mode, bending in step, is stable across the same band.

A wake ends by bending, not by fading

The two vortices behind an aeroplane do not simply weaken until they are gone. In still air each one bends in the other's strain, the bend grows by a factor of e every thirty-three seconds behind an airliner, and after a couple of minutes the pair has pinched itself into a chain of rings. The wavelength it chooses is eight and a half times the spacing, and the vortex cores hardly enter.

circulation · Tip vortex
The drip spacing is set by the depth of the layer. The spacing of the fastest-growing wave, in units of 2π capillary lengths, against the depth of the hanging layer in capillary lengths, for water, glycerol and honey. A thin film of any of them drips at √2. Deep water drips at √3. Deep glycerol and honey keep growing past √3, their fastest wave as long as the layer is deep, because viscosity slows short waves more than long ones.

How far apart a ceiling drips

A layer of liquid hanging from a ceiling is heavy fluid over light, and every ripple on it longer than about seventeen millimetres grows. Which ripple grows fastest, and so how far apart the drips form, is usually given as one number. It is at least three, and what chooses between them is not the liquid's surface tension but the depth of the layer.

turbulence · Instability
The unburnt gas bends round a wrinkle, and the bending makes it grow. Streamlines of the unburnt gas flowing up towards a flame with a small sinusoidal wrinkle, one wavelength across, for a density ratio of seven, with the wrinkle exaggerated. The expansion behind the flame pushes back on the gas ahead of it where the flame bulges forward, so the streamlines spread there and the gas arrives slower, and they crowd together where the flame lags, so the gas arrives faster. A flame that burns into the gas at a fixed speed therefore advances further where it was already ahead.

A flat flame is unstable at every size

A flame expands the gas it burns, and the expansion pushes back on the fresh gas ahead. Where the flame bulges forward the fresh gas slows, so the bulge burns further forward; where it lags, the gas speeds up and it falls further back. Every wrinkle grows, the shorter ones faster, and what finally gives a real flame a size is how its burning speed responds to its own curvature.

kinematics · Dilatation
Two eddies that stand still behind a cylinder in ideal flow. A uniform stream past a circular cylinder with a pair of opposite vortices standing 2 radii behind its centre and 0.859 either side of the axis, each of strength 10.31 in units of the stream speed times the radius — Föppl's equilibrium. The shading is the stream function: the fluid between the vortices and the cylinder circulates in two closed eddies and never leaves, while the stream divides round the whole assembly as if it were one longer body.

Two eddies can stand behind a cylinder, but not for long

Ideal flow, which has no viscosity and no wake, can still hold a pair of eddies standing behind a cylinder — at any distance behind it, with a strength fixed by the distance. They cost the cylinder nothing. And they cannot stay: nudged sideways by a thousandth of a radius, the pair grows its displacement exponentially and leaves, which is the first step of shedding a wake.

inviscid · Images
Whether a sloping ceiling drips in place is one ray's growth. How fast a disturbance grows as seen by an observer moving along the ceiling at speed v, for four speeds at which the film carries its disturbances downhill, in the film's own units. On a flat ceiling (V = 0) the growth peaks at a quarter for the observer standing still and falls to zero for one running at ±1.622. Tilting the ceiling slides the whole curve downhill. While the observer at the point disturbed, v = 0, still sees growth, the ceiling drips where the disturbance began; at V = 1.622 that observer sees none, and above it the disturbance grows only while it is carried away.

A sloping ceiling drips downhill, or not at all

A film hanging from a level ceiling drips where its ripples form. Tilt the ceiling and the film flows downhill, carrying its ripples with it, and at some slope they are carried away faster than they spread and the ceiling stops dripping in place. That slope is set by one number, 1.622 — the speed at which a level ceiling's disturbance spreads — and for water it is tiny: a degree for a tenth of a millimetre of film, six for half a millimetre. Past it the ripples still grow, and still drip, but downhill, at a distance that grows with the slope: a ceiling shorter than that delivers its water to the edge.

turbulence · Instability

A strained vortex holds until it has no shape to hold

A patch of vorticity in a strain has two steady shapes for every strain below 0.150 of its vorticity, a rounder one and a longer one, and none above. The longer one comes apart at the slightest nudge. The rounder one, computed with disturbances of two, three, four and five lobes, never does: it nods and holds right up to the strain at which it ceases to exist. So the existence limit is the real limit, and past it a vortex is not shattered but stretched — lingering first near the shape it has lost, for a time that grows as the fourth root of how close the strain is to the limit.

inviscid · Vortex patch

Four vortices bend faster, and bend the wrong one

With its flaps down a wing sheds four vortices, not two: a strong one at each tip and a weaker one at each flap's edge, and on each side the two circle one another as the wake sinks. The orbit makes the wake unstable at waves a fifth of Crow's length, growing twice as fast. But the wave that grows is the weak flap vortex being bent round the strong tip vortex, which hardly moves, while the long bend that pinches the tips together grows no faster than before. Only a counter-rotating inner vortex makes the whole wake unstable, and it does so at every wavelength.

circulation · Tip vortex

A thin interface keeps its waves past a quarter

Miles' quarter rules a shear layer whose density changes over the same depth as its velocity. Make the density interface three times thinner and the stationary billow dies early, but a pair of travelling waves takes its place and is still growing at five times the quarter. No theorem is broken: at the edges of the shear, where the waves draw their energy, the local Richardson number has fallen to nothing.

turbulence · Stratification

The fastest finger is set by the gap and one number

Push air into oil between two glass plates and the interface breaks into fingers. The averaged law of the cell makes the instability exact in its linear stage: every wavelength is driven in proportion to its wavenumber, surface tension holds back the short ones as the cube, and the fastest finger is π gap widths divided by the square root of a capillary number. A narrow channel stays flat; a wide one chooses how many fingers to make; and a growing bubble of air makes more of them the larger it gets.

inviscid · Hele-shaw

What holds the shear decides whether the waves last

Holmboe's travelling waves need a density interface about half as thick as the shear around it, and diffusion thickens both. If the shear layer is held at its thickness, the interface catches up — in under a minute for heat in a laboratory tank, over an hour for salt — and the waves get two e-folds or a hundred and seventy. If the shear layer is free it thickens faster than the interface, the ratio heads to the root of the Prandtl number, and the waves last until the Richardson number ends them instead.

turbulence · Stratification

The meniscus slows the finger and hardly moves it

Air pushing oil through a narrow gap fingers at a wavelength the gap and one number set, on the assumption that surface tension acts only on the interface's curve in the plane. The interface is really a meniscus across the gap, leaving a film behind, and Park and Homsy's correction for it does two things that nearly cancel in the wavelength: it cuts the tension to π/4 and adds a drag that grows with speed. At bench speeds the fastest finger barely moves; it grows a quarter more slowly, and a channel goes unstable at four-fifths of the speed.

inviscid · Hele-shaw

Named alongside it

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

Model limitEigenvalueBuoyancyShear layerInflection pointReynolds numberSurface tensionThresholdWakeDispersion relationInstabilityRayleigh–Bénard convection

All concepts