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

Scaling — where it appears

The way a quantity depends on the size or speed of the flow, usually as a power. Establishing an exponent takes far less than solving a problem and is often all that is needed, since it says what changes when the aeroplane is built full size.

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

5 quantities, 3 rows, 2 left over. The dimension matrix for the drag on a sphere: one column per quantity, one row per base dimension, and every entry an exponent. Buckingham's theorem is a statement about this matrix and nothing else — the number of independent dimensionless groups is the number of columns minus the rank, computed here by elimination. Nothing about fluids enters until somebody decides which columns to write down.

Counting what matters

Five quantities decide the drag on a sphere, and the experiment that measures it has one curve in it rather than a five-dimensional table. The reason is a rank: the matrix of dimensions has three independent rows, and what is left over is the number of dimensionless groups the answer can possibly depend on.

regimes · Dimensional The bath is ten thousand times too small. The Rossby number — the ratio of the inertial term to the Coriolis term in the momentum equation — for eight flows at 45 degrees, on a logarithmic axis. Above one, rotation is a correction; below one, it is the physics. A draining bath sits at 3.2e+3 and a mid-latitude depression at 1.9e-1. The Coriolis term is not absent from the bath: it is present, computable, and four orders of magnitude smaller than the terms that decide what happens. Saying so is not the same as saying it is zero.

The bath does not know the hemisphere

The Coriolis term in a draining bath is not absent. It is present, computable, and about four hundred times smaller than what a hand left in the water an hour ago is still doing — which is a ratio rather than an opinion, and it is the same kind of argument as a Reynolds number.

misconceptions · Coriolis scale 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 1.139 asks for a Francis. The specific-speed axis, with the four machine types on it and one duty marked: 3 m³/s at 60 m, on a shaft turning at 750 rev/min. The number is 1.1387, and the choice of runner follows from it before any blade has been drawn. What the number contains is a ratio of flow to head; what it does not contain is any size at all, which is why one axis serves a garden pump and a gigawatt turbine.

One number picks the machine

A flow rate, a head and a shaft speed contain exactly one dimensionless combination with no size in it. That combination decides whether a duty wants an impulse wheel, a Francis runner or a propeller — before anything has been drawn, sized, or costed.

applied · Specific speed The shock radius against time, from an equation that was given no exponent. The thin-shell energy balance integrated forward from a small initial radius, on logarithmic axes. The two-fifths power is not put in: the ordinary differential equation is Ṙ = √(E/AρR³), and the straight line is what it does. The fitted slope is 0.39983 and the fitted prefactor is 0.90721 against a closed form of 0.90702 — the small residuals being the integration's memory of where it started, which the similarity solution has no equivalent of.

A radius that gives the energy away

Four quantities, three dimensions, one group. The radius of a strong blast must be a constant times (Et²/ρ)¹⁄⁵, and nothing about the device, the chemistry or the initial size can appear. The exponent is free and the constant is not.

compressible · Blast wave Three bluff bodies, whose Strouhal numbers agree once the wake's width is used. The measured Strouhal number of each body, and Roshko's universal number formed with the wake's width and the speed on the free streamline that bounds it. The raw numbers span a factor of 1.462; the collapsed ones span 1.0011, with a mean of 0.16281. The shedding was never body-dependent — the length in the number was.

The frequency a wake chooses

Bluff bodies shed at Strouhal numbers from 0.145 to 0.212, and the spread is not a fact about shedding — it is a fact about which length went into the number. Change the length to the wake's own width and three bodies agree to a tenth of a per cent. Then let the body move, and the number stops deciding anything at all.

regimes · Strouhal Two different moments of one distribution. The permeability and the specific surface of a log-normal bundle, against the width of its pore-size distribution at a fixed median. The permeability rises by five orders because it is a fourth moment and the widest tubes dominate it; the surface falls because it is a first moment and the narrowest tubes dominate that.

A permeability that is only the geometry

Kozeny–Carman says that a porous medium's permeability follows from its porosity and its specific surface. Both are real, both are exactly measurable, and they do not determine the answer: forty-nine tube bundles built with identical values of each span a factor of eighty-two in permeability.

applied · Porous The junction the minimisation is over. A parent vessel entering from the left and two daughters leaving to fixed points. The radii are settled by Murray's law; what is left free is where the branch point sits, and the cost of the junction depends on it. The point drawn is the one the minimisation finds.

The angle a junction chooses

Murray's law fixes the radii at a branching vessel and is where every account of it stops. The same minimisation fixes the angles completely — 74.93 degrees for a symmetric bifurcation, a right angle for a vanishing side branch — and it does so as a triangle of forces, with tensions proportional to the squares of the radii.

applied · Branching Nu/√Re against the Prandtl number, over eight decades. The whole of the flat plate's heat transfer, as one curve. It is not a power law: it goes as Pr^½ at the bottom, where the thermal layer is far thicker than the viscous one, and as Pr^⅓ at the top, where it is buried inside it. The Pr^⅓ everybody quotes is the upper half. The two asymptotes are drawn beside it, and the low one is √(Pr/π) in closed form.

The number that is an answer

Almost every dimensionless group is a hypothesis: somebody sets the speed, the size and the fluid, and the number licenses a model. The Nusselt number is not. It is what the experiment produces, it sits on the left of the equals sign, and a regime diagram drawn on it is a category error.

regimes · Nusselt Nu/Gr^¼ against the Prandtl number, with the exact solution's points on it. The closed form 0.508 Pr^½(20/21 + Pr)^−¼, over eight decades, with Ostrach's exact similarity values marked. The integral method is two to eight per cent high from Pr = 0.7 upwards and 27 per cent high at Pr = 0.01 — which is where the thermal layer is ten times the momentum layer and giving them one thickness stops being an approximation to anything.

A speed nobody imposed

Every regime number in this collection contains a velocity somebody chose. Natural convection has none: a warm plate makes its own flow, and the Grashof number is what is left when the speed is taken out. The Reynolds number of the result — six thousand, on an ordinary radiator — is an output of the solution rather than a setting on an apparatus.

regimes · Grashof Three exponents for one dimensionless group. The local slope of each error, measured over one decade at a time. The duct's is exactly 1, the long wave's is exactly 2, and the slender body's drifts from 1.900 to 1.733 across the range and never reaches either. The same geometric ratio, in three problems that look alike, and the third one has no exponent at all.

One group, three exponents

Lubrication theory, shallow water and slender-body theory are taught in three places and are one expansion in one group — a ratio of two lengths, with no speed, no viscosity and no fluid in it at all. The error is supposed to be second order. In three problems that look alike it is first order, second order, and an exponent that does not exist.

regimes · Slenderness The convergence exponent depends on the gas, which a dimensional exponent cannot. R ∝ (−t)^α for a converging shock, against the ratio of specific heats, for cylindrical and spherical symmetry. Guderley's exact values are marked and the agreement is to four figures. The Sedov blast's two-fifths is drawn beside them: it is the same for every gas, because it comes from dimensions and a conserved energy, and γ is dimensionless.

An exponent dimensions cannot give

A blast wave's radius goes as the two-fifths power of time, and the two-fifths is arithmetic: count the dimensions and it falls out. A shock converging on a point goes as the 0.717 power, and no amount of counting will produce that number — because it depends on the gas, and γ is dimensionless.

regimes · Similarity 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 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 Two ways to move a duty, along one axis, in opposite directions. A duty at a specific speed of 0.01 and what splitting it does. Dividing the head between stages in series multiplies each stage's specific speed by the number of stages to the three-quarter power, because the group carries (gH)^(−3/4); dividing the flow between units in parallel divides it by the square root of the number of units, because the group carries √Q. Both exponents are read off the group rather than remembered, and they are why the two operations are not interchangeable: three stages buy a factor of 2.28 and three units cost a factor of 0.58. The bands are drawn in the colour this site reserves for a borrowed claim, because where a Francis runner stops is practice rather than a result.

The duty that had no machine

Some duties have a specific speed outside every band, and no runner will do them at any size because the number contains no size. That is true and it is not the end. The same group says how to split the duty until it fits, and the two ways of splitting move it along the same axis in opposite directions, by exponents read straight off it.

applied · Specific speed And they are one piece of arithmetic. All three corrections against one over the logarithm of their own large parameter. Each is a straight line through the origin, with its own slope: one third exactly for the enstrophy range, near one for the wall layer, and near ln 10 for Oseen's per-decade measure. Three disappointments in three different fields, written up separately, are the same function.

The three that never converge

Oseen's drag coefficient, the overlap layer's power-law exponent and the enstrophy range's slope are three separate disappointments in three different fields. They are one piece of arithmetic, and the arithmetic says why none of them will ever be reached.

regimes · Crossover One shape, at every station and for every jet. The axial velocity of Schlichting's round jet at four distances, each scaled on its own centreline speed and its own half-width. The four curves are one curve: the shape function (1 + eta²/4)⁻² has no parameter in it at all, and every station of every jet — laminar or modelled, weak or strong — collapses onto it exactly.

Exactly similar, and one number short

The round jet has an exact solution of the Navier–Stokes equations, and a turbulent jet is modelled by the same formula with the viscosity replaced by a fitted constant. The shape is identical. Nothing a measurement of the profile can do will tell the two apart.

viscous · Free shear A blade root at 412.51 megapascals, and the chord is not in it. The centrifugal stress at a blade root against shaft speed, for one annulus of 0.36 m², in four materials. Integrating the blade's own weight outward gives σ = 2π ρ_b A N² with a taper relief — and the chord has cancelled, the blade count has cancelled, and every property of the gas has cancelled. What is left is an area times the square of a speed, capped by a material. The horizontal lines are each material's allowable stress, and where a curve crosses its own line is the fastest that annulus may be turned in that metal. At 12000 rpm this blade carries 412.51 MPa with a tip speed of 439.82 m/s.

The stress that picks the aerodynamics

The free term in a rotor's pressure rise wants radius and speed, and both are capped by something with no fluid in it. Integrate a blade's own weight outward and the root stress comes out as the annulus area times the square of the shaft speed — with the chord cancelled, the blade count cancelled, and every property of the gas cancelled.

applied · Turbomachine The same data, in a basis the theorem allows just as much. The identical two hundred and forty points, plotted as F/(mu U d) — which is the drag coefficient times the Reynolds number, a perfectly legitimate pi group forming a complete pair with the abscissa. The result is a straight line of slope 0.9985 through five decades with an r-squared of 0.9985, and it contains no physics: the ordinate contains the abscissa.

The groups are not the only groups

Buckingham's theorem fixes how many dimensionless groups an answer can depend on and says nothing about which. Two of the infinitely many legitimate choices are used here on the same data: one manufactures a straight line through five decades out of a constant, and the other erases Stokes' law completely.

regimes · Dimensional How far apart two points can be and still be correlated. The correlation between the logarithm of the dissipation at two points, against how far apart they are in units of the smallest scale, over four decades. It falls as a ratio of logarithms — so it is still a quarter at a thousand smallest scales, and reaches a half only at a hundred.

A dissipation correlated across every scale

The dissipation is supposed to be the most local quantity in turbulence — a thing happening at the smallest eddies, everywhere and independently. A multiplicative cascade makes its logarithm correlated over a distance that is a ratio of logarithms, so two points a thousand smallest scales apart still agree a quarter of the time.

turbulence · Intermittency A pump slowed into a system with a static lift leaves its own specific speed. The specific speed of the operating point, as a share of its value at the best point, against the fraction of the design flow delivered, for a pump drawn for 0.1 m³/s against 40 m at 1450 rpm, specific speed 0.545 at its best point. Under speed control into a system whose static lift is 0 per cent, 30 per cent, 60 per cent, 90 per cent of the design head (lines), and under a throttle at design speed into the 60 per cent system (dashed). With no static lift the speed-controlled pump stays exactly at its best point and its specific speed never moves. At half the design flow it has fallen to 1.000, 0.800, 0.708, 0.652 of the design value as the static share rises, and to 0.598 under the throttle.

The specific speed a pump spends its life at

A pump is chosen by its specific speed at its best point and then run somewhere else. Written in the pump's own coefficients the number is √φ/ψ^¾, a position along its characteristic, and a variable-speed drive keeps it there only when the system it pumps into has no static lift. Every metre of lift moves a slowed pump along its own curve, towards shut-off, and a throttle moves it further.

applied · Specific speed 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

Named alongside it

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

MeasurementModel limitCorrelationDimensionless numberDimensional analysisDimensionlessSimilarityDissipationOptimisationRegimeReynolds numberSimilarity solution

All concepts