Fluid Flow & Aerodynamics
why things fly, and why the usual explanation is wrong
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 18 essays across 6 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 this site cannot draw

The alternating wake behind a cylinder is the most photographed structure in fluid mechanics, and this site's solver 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

Named alongside it

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

EigenvalueModel limitInflection pointReynolds numberRayleigh–Bénard convectionShear layerThresholdTransitionBoundary layerBuoyancyConvergenceMeasurement

All concepts