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

Threshold — where it appears

A value of a parameter at which a flow's behaviour changes qualitatively rather than gradually. Where one exists it is usually the sharpest thing in the problem, and computing it exactly is often easier than computing anything on either side of it.

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

Particles at St = 1, against the flow that carries them. Particle paths and the streamlines they were released on, in this site's exact cylinder solution. At small Stokes number the two are indistinguishable and the body catches nothing; as the particles get heavier their paths straighten, cross the streamlines, and begin to strike. The paths are integrated with Stokes drag and nothing else — no gravity, no lift, no effect of the particles on the flow.

Whether the droplet turns

The air goes round the wing. Whether what is carried in it goes round too is decided by one number — and below a critical value of that number the body collects nothing at all, however many droplets are thrown at it, because the flow turns every one of them in time.

regimes · Particle Where the number says, and where it happens. Fourteen dimensionless groups on one logarithmic axis. The open circle on each row is the value at which the two terms the group compares are equal, which is one by the way the group is formed; the filled mark is the value at which the thing a reader cares about first changes by 1%. The bar between them is the distance the folklore phrase "of order one" hides, and it runs from nothing at all to a factor of 594.

What "of order one" is worth

A dimensionless group is built by comparing two terms, so it is one when the terms are equal — and that is the only thing it says. Where the behaviour actually changes is a separate question with a separate answer, and across fourteen groups on this site the two numbers differ by factors from one to five hundred and ninety-four.

regimes · Crossover How wrong the ordinary answer already is. The error in a no-slip continuum calculation of the flow through a channel, against the Knudsen number, both logarithmic. The rule at Kn = 1 is where a molecule crosses the whole channel between collisions — the value the number is named for. The error is one per cent at Kn = 1/594, five per cent at 1/114 and ten at 1/54, so the continuum regime of the usual classification, which runs to Kn = 0.01, is a region in which the continuum answer is already six per cent out at its far end.

Where a fluid stops being one

The Knudsen number is the mean free path over the size of the thing, and at one a molecule crosses the whole channel between collisions. The continuum equations with a no-slip wall are already one per cent wrong at one part in five hundred and ninety-four, which is a factor nothing about the definition would suggest.

regimes · Knudsen Quasi-steady stops being true a long way before one. The magnitude of Theodorsen's function, which is the factor a quasi-steady lift calculation is wrong by, and the phase the lift lags the motion. Quasi-steady means C = 1, and the amplitude is already one per cent low at k = 0.0061 and fifteen per cent low at k = 0.1 — a reduced frequency at which nobody hesitates to call a flow quasi-steady. The lag is worse: it reaches a degree at k = 0.003, and a flutter calculation is decided by phase rather than by amplitude.

Slow enough to be steady

A wing moving slowly enough is assumed to carry the lift its instantaneous angle asks for. The reduced frequency has two thresholds — one where the apparent-mass and circulatory lifts are equal, and one where the quasi-steady answer stops being right — and they are a hundred and seventy-eight apart.

regimes · Reduced frequency How small is small enough. The error in Stokes' law for the drag on a sphere, against the Reynolds number, both logarithmic. Oseen's correction is the first term the neglected inertia puts back, and it says the error is 3Re/16: one per cent at Re = 16/297 = 0.054, five per cent at 0.28, and already sixteen per cent at Re = 1 — which is the value at which the two terms the Reynolds number compares are equal, and is where every textbook draws the boundary of creeping flow.

How small is small enough

Creeping flow drops the inertia terms and gets an exact answer for the drag on a sphere. The Reynolds number at which those terms are equal is one — and at one, Stokes' law is already sixteen per cent low. The honest boundary is 0.054, and the reason it is so far down is the same reason the theory needed repairing in the first place.

regimes · Reynolds How big before gravity shows. A drop's height over its width against the Bond number, which is the ratio of its weight to the force its own skin can supply. The number is one where those two are equal, and by then the drop is a bun: it is one per cent from a ball at Bo = 0.0079, five per cent at 0.054 and ten at 0.13. Every one of those is below one, and the first is below it by a factor of a hundred and twenty-six.

The size a drop is allowed

The Bond number sets a drop's weight against the force its own skin can supply, and it is one when they are equal. By then the drop is a bun — it is a per cent from being a ball at Bond number 0.0079, which is a water drop half a millimetre across.

regimes · Bond A falling drop is a bun, not a tear. The shape of a drop pressed on by the air it is falling through, computed by matching the Legendre component of a sphere's own potential-flow pressure distribution against the change in curvature it produces. The result is oblate — flattened along the direction of travel — because the pressure is high at the poles and low round the equator. The teardrop of every illustration has the wrong symmetry entirely: a drop in free fall has no up, and one at terminal speed is being pressed on from in front.

The drop that is not a tear

A falling raindrop is flattened along the direction it is going, by the pressure of the air passing it rather than by its own weight, and the group that decides is the Weber number. The teardrop of every illustration has the wrong symmetry entirely — there is no up in the problem it is drawn for.

regimes · Drop shape Which speed the number is formed on. The fractional change in air density at three places on a body, against the free-stream Mach number. At a stagnation point the density rises, and it reaches five per cent at M = 0.314 — which is where the familiar 0.3 comes from, and it is a five per cent tolerance rather than a physical boundary. At the suction peak the density falls instead, and how fast depends on the body: a lightly loaded section is milder than its own nose, and one working at cp₀ = −2 reaches five per cent at M = 0.22 and is at Mach 0.55 over its shoulder while the free stream is at 0.3.

Which speed goes in the number

The most quoted threshold in the subject — air is incompressible below Mach 0.3 — is a five per cent tolerance on the density at a stagnation point wearing a physical boundary's clothes. A wing working for its living is at Mach 0.55 over its shoulder while the free stream is still at 0.3.

regimes · Mach The error in the balance is the number itself. The fractional error in the geostrophic wind, against the Rossby number, on logarithmic axes. It is a straight line of slope one through the origin, and that is not an approximation: keeping the centripetal term gives V_g/V = 1 ± Ro exactly, so the error and the number are the same quantity. The geostrophic wind is one per cent right at Ro = 0.01 and a hundred per cent wrong at Ro = 1, which is the value the number is named for and is quoted as the boundary of the approximation.

The balance that is its own error

Geostrophic balance is licensed by the Rossby number being small, and the fractional error in the geostrophic wind is the Rossby number — exactly, not approximately. So the balance everybody uses at Ro of order one is a hundred per cent wrong, and the same quadratic has a hard limit at a quarter that no anticyclone can pass.

regimes · Geostrophy The relaxation time a particle actually has. Two quantities against the particle-to-fluid density ratio. β = 3ρ_f/(2ρ_p + ρ_f) is three for a bubble, one for a neutrally buoyant particle and nearly zero for anything heavy; it is the factor by which the fluid's own acceleration is felt. The other curve is the true relaxation time over the usual formula's, which is one for a heavy droplet, exactly three halves for a neutrally buoyant tracer, and unbounded for a bubble — the usual formula gives a bubble a relaxation time of zero, and therefore no dynamics at all.

The tracer that is not one

Every Stokes number is built on a relaxation time that counts the particle's own inertia and nothing else. Adding the two terms it leaves out gives a bubble a relaxation time where the usual formula gives zero, makes a neutrally buoyant tracer half as slow again as advertised, and sends bubbles into vortex cores that droplets are flung out of.

regimes · Particle 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 Where a profile stops being a parabola. Two measures of how far Womersley's solution has left the quasi-steady parabola, against the Womersley number, both logarithmic. The amplitude deficit reaches a hundredth at α = 0.28 and the phase lag reaches a hundredth of a right angle at α = 0.22 — both well below α = 1, which is where the unsteady and viscous terms are equal and is the value the number is named for. By α = 1 itself the flow is already a tenth short and eight degrees late.

Where the parabola goes

A pipe carrying a steady flow has a parabolic profile, and the Womersley number is supposed to say when an oscillating one still does. For the flow rate it is very nearly honest — one per cent at α = 0.91 — and for the phase it is out by a factor of three, because a lag is second order in the number and an amplitude deficit is fourth.

regimes · Womersley Five numbers, one flow, one name. A hyperbolic-tangent shear layer with a matching density profile, measured five ways, each of which appears in the literature as "the Richardson number". The minimum gradient value is what Miles' and Howard's theorem is about and is the only one the quarter belongs to. The bulk numbers depend on which thickness and which velocity difference were used; the depth-averaged one depends on how far from the layer the measurement extended, and grows without limit as it extends further.

Five numbers, one name

The Richardson number has a threshold at a quarter, and the quarter belongs to one of the five quantities that go by the name. On a single tanh shear layer they run from J to seventeen J, and the largest of them grows without limit as the measurement is extended further from the layer.

turbulence · Stratification The number a duct settles at is an eigenvalue. The local Nusselt number against x⁺ = x/(D·Re·Pr), from a Crank–Nicolson march that contains no eigenvalue anywhere. It settles at 3.6568, which is λ₀²/2 for the Graetz eigenvalue problem — a completely separate calculation. The mark at x⁺ = 0.05 is the entry length every textbook quotes: it delivers a Nusselt number 1.45 per cent above the developed value, which is a perfectly reasonable tolerance and is never the one stated.

How far before the heat arrives

A duct's thermal entry length is quoted everywhere as x/(D·Re·Pr) = 0.05, with no tolerance attached. Working out what it delivers gives a Nusselt number 1.45 per cent above the developed value — and the developed value itself is not a term ratio at all but an eigenvalue, 3.6568, which is also the rate at which the duct forgets its inlet.

regimes · Peclet A limit that exists and is never reached. The exponent of the best power law fitted across the overlap layer, against the friction Reynolds number. A logarithm is the zero-exponent member of that family, so the log law is what this sequence is heading for — and it heads there as 1/ln Re_τ, which is the slowest useful way of approaching anything. The exponent is still 0.102 at Re_τ = 10⁶, and driving it to a hundredth needs a Reynolds number with a hundred and fourteen in its logarithm.

A limit nothing reaches

A dimensional argument that succeeds says a variable has dropped out of the answer. The Blasius profile has no Reynolds number in its shape at any Reynolds number; the overlap layer's power-law exponent is still 0.102 at Re_τ of a million and falls as a logarithm, so the limit exists and nothing ever gets there.

regimes · Dimensional The pressure drop stops rising at 0.213 m/s. The pressure drop across a bed of 500 µm sand, against the velocity through it, in units of the fluidisation velocity. The rising branch is Ergun's resistance and the flat one is the bed's buoyant weight, which the flow cannot exceed however hard it is pushed: past the corner the bed expands rather than resisting more. That flat line is the reason fluidisation is unmistakable in practice — the corner is a crossing of two curves rather than a gradual departure, and it can be read off a gauge.

The bed that weighs itself

Blow hard enough through a pile of sand and the pressure drop stops rising. It cannot rise — a control volume round the bed says the drop can never exceed the buoyant weight of the solid in it, and at the velocity where the two meet the bed stops being a structure and starts being a fluid.

applied · Porous The one number that really is one. Three quantities against the Froude number. The upper line is the speed of a surface wave travelling downstream and the lower one the speed of the same wave travelling upstream, both in units of the wave speed itself; the second changes sign at Fr = 1 and not near it. That sign change is not a comparison of two term sizes going through unity — it is the moment a signal stops being able to reach upstream at all, so the equations change from elliptic to hyperbolic and the flow stops knowing what is ahead of it. The specific energy, drawn beneath, has its minimum at the same place, and for the same reason.

The number that really is one

Almost every threshold in this subject sits somewhere other than where its dimensionless group is one. The Froude number does not. At Fr = 1 a disturbance stops being able to travel upstream, the specific energy is least and the equations change type — three statements, one number, and no tolerance anywhere in it.

regimes · Froude Storage and loss moduli, for one relaxation time and for a spectrum. The two moduli of a Maxwell fluid and of a Rouse chain, against frequency in units of the longest relaxation time. One relaxation time makes the storage modulus overtake the loss modulus at λω = 1 and then leave it behind without limit. A spectrum of 1000 modes makes them rise together as the square root of frequency and stay a fixed ratio apart, so the material never becomes the solid the single time predicts.

A solid, if it is not given time

The Deborah number is the only group on this site with no fluid in it — two times and nothing else — and it says a material is a solid or a liquid depending on how long anybody watches. What it throws away is that no real material has one time, and the spectrum it replaces changes the answer in kind rather than in degree.

regimes · Deborah 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 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 A threshold that is a curve, not a pressure. The tension a rectangular pulse has to reach to make a five-micron bubble run away, against how long the pulse lasts. At a tenth of a microsecond it takes forty bar; at a hundred microseconds it takes 1.05, which is within about one per cent of the static threshold. There is no such thing as the cavitation pressure of this bubble on its own.

A threshold that is also a duration

The cavitation number treats inception as a pressure: below it the liquid tears, above it does not. A five-micron bubble asked to grow in 0.3 microseconds needs 43 bar of tension and the same bubble given a hundred needs 1.05, because it has to make a journey and not merely respond.

applied · Cavitation The exact solution and its three approximations, at ε = 0.02. The outer solution is excellent everywhere except in a layer of width ε at the left, where it is wrong by a whole unit. The inner solution is excellent inside that layer and wrong everywhere else. The composite is their sum less the part they agree about, and it is within order ε of the exact solution across the whole interval — which is the entire content of matched asymptotics, drawn.

One formula for both ends

Two limits, each with its own description, neither valid everywhere. The composite is the sum less the part they agree about, and it is uniformly good — but the overlap region that justifies the construction does not exist at ε = 0.01, and the composite is still accurate to two per cent there.

regimes · Crossover The core that does not move, and where its edge is. Four velocity profiles at the same pressure gradient and four yield stresses. The shear stress in a pipe is G r/2 whatever the fluid is — that is a force balance and not a constitutive law — so the fluid is unyielded exactly inside r = 2 tau_y/G, and the plug radius is known before anything is solved.

The core that does not move

Toothpaste in a tube has a region in the middle that is not being sheared at all, and its edge is at a radius you can write down before solving anything. Four things about that core are exact, and one of them is that the fluid does not flow below a threshold — not slowly, at all.

viscous · Non-newtonian 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 profile a diverging channel flattens into, and then cannot hold. Five purely outward profiles in a wedge of 0.2 radians, at rising flux, each normalised to its own centreline value. As the flux rises the profile flattens in the middle and steepens at the walls — and then it stops. The last one has zero slope at the wall, which is separation, and beyond it no purely outward profile of this form exists at all. Nothing was added to the equation to make that happen: the wall shear is the square root of a cubic and the cubic runs out.

One channel, one flux, two flows

Flow between two plane walls meeting at a line has an exact solution. Past a threshold that turns out to be a ratio of gamma functions, it has two — the same wedge carrying the same flux, once outward everywhere and once with the fluid running backwards along both walls, and nothing in the equations chooses.

viscous · Exact layer The arrival map, and the place where it goes backwards. Where each boom ray lands on the ground, against the Mach number the aeroplane was doing when it launched it, for four accelerations from 15 km. In level flight this would be a straight line rising at the aeroplane's own speed. Here it falls before it rises: a ray launched at a higher Mach number is shorter and steeper, and near the cut-off it shortens faster than the aeroplane advances. Every minimum in these curves is a fold — two emission times delivering to one place — and on the fold itself neighbouring rays converge onto a single line.

The carpet an accelerating aeroplane folds

Level flight launches every boom ray with the same invariant, so they run parallel and each place hears one boom. Accelerate, and each successive ray is shorter than the last — shorter, near the cut-off, than the aeroplane's own advance — so later rays overtake earlier ones and the arrival map folds onto a line.

compressible · Sonic boom The characteristic gets a wall, and the wall does not care about the discharge. One jet pump's characteristic, at an area ratio of 0.275, with the flow ratio at which its throat entry reaches vapour pressure drawn for three values of the cavitation parameter σ = (Pₛ − pᵥ)/(Pₘ − Pₛ). Left of a wall the machine runs on its curve. At the wall no lower discharge pressure raises the flow: the head ratio can fall to zero along the vertical and the flow ratio stays where it is. The wall's position contains the nozzle and suction losses and nothing downstream of the throat entry.

The wall the suction puts in the curve

A liquid jet pump's lowest pressure is where the entrained stream enters the throat, and when that reaches vapour pressure the pump curve stops being a curve. The flow ratio freezes at a value no lower discharge pressure can move — and raising the motive pressure, the obvious cure, brings the wall closer.

applied · Ejector Where the eddies outconduct the molecules. The ratio of turbulent to molecular heat diffusivity across the pipe at Reτ = 2000, for five fluids, on logarithmic axes. Wherever it is above one the eddies carry more heat than conduction does. For air it passes one inside the buffer layer and reaches 135; for water, earlier and higher. For liquid sodium it never reaches one anywhere: at its peak, halfway to the axis, the eddies carry just over half what conduction carries, and the temperature profile is set by conduction across the whole pipe.

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 A boom gathers most of its age in the thin air near the aeroplane. The share of the total age gathered above each height, for the ray under the track and the last ray computed near the carpet's edge, with the share of the path length travelled above each height for comparison. Under the track 66 per cent of the age is gathered above the tropopause in 26 per cent of the path. The edge ray spends most of its path in the lowest few kilometres and gathers only 11 per cent of its age below 3 km, because the same pressure distorts dense air far more slowly than thin air.

A boom is aged in the thin air it starts in

The rays that reach the edge of a sonic-boom carpet travel two and a half times as far as the one under the track, and it is natural to expect their signatures to have aged accordingly. They have not. A pressure wave distorts thin air far faster than dense air, so two-thirds of a boom's ageing is done in the stratosphere near the aeroplane, and the extra kilometres near the ground add little — which decides how far out a boom shaped to be quiet stays quiet.

compressible · Sonic boom

Named alongside it

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

DimensionlessMeasurementToleranceEigenvalueModel limitRegimeAsymptoticsConvergenceDimensionless numberLinear stabilityModel validityBoundary condition

All concepts