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

Similarity — where it appears

The property of two flows that differ only by a change of scale, so that one is the other measured in different units. It holds when the governing dimensionless numbers match, and matching all of them at once is usually impossible.

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

Reynolds number: one number, four different flows. Reynolds number is inertia ÷ viscosity. It is not a property of the fluid or of the shape but of the combination, and crossing a threshold changes the physics rather than the magnitude.

One number decides which physics applies

A bacterium and a whale both swim, and they are not doing the same thing at different sizes. The ratio of inertia to viscosity separates them, and crossing it changes the rules rather than the magnitudes.

regimes · Reynolds Reynolds number: one number, four different flows. Reynolds number is inertia ÷ viscosity. It is not a property of the fluid or of the shape but of the combination, and crossing a threshold changes the physics rather than the magnitude.

The Reynolds number, and the length in it

The most useful number in fluid mechanics has an arbitrary quantity buried in it, and quoting one without saying which length was used makes it meaningless. That detail is where most misuse comes from.

regimes · Reynolds 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 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 correction belongs to the wing, not to the air. Lift-curve slope against Mach number for three aspect ratios, each computed by solving the lifting line for the wing the transformation actually implies — one of aspect ratio βAR — and dividing by β. The dashed curve is the two-dimensional rule, which is what an aerofoil section gets. At Mach 0.7 the aspect ratio of 20 has gained 35 per cent of slope and the aspect ratio of 4 has gained 24, against the 40 per cent the section rule promises both. The β in the finite-wing term cancels the β in front of it, so the shorter the wing the less compressibility does to it.

The wing the equation is really solving

The Prandtl–Glauert rule is usually quoted as a factor on the answer. It is a change of shape — and in three dimensions the shape it changes is the aspect ratio, so a short wing is far less affected by compressibility than a long one.

regimes · Similarity 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 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 coefficient of the equation's second derivative, along a chord. The bracket multiplying the streamwise second derivative in the transonic small-disturbance equation, along a chord at Mach 0.85. Where it is positive the equation is elliptic and the flow is subsonic; where it is negative the equation is hyperbolic and the flow is supersonic. Which it is at a given point depends on the perturbation velocity there, which is the thing being solved for. Forty-two per cent of this chord is hyperbolic, and no amount of inspecting the problem beforehand could have said so.

The equation that changes type inside its own answer

Near Mach one the coefficient of the streamwise second derivative depends on the perturbation velocity, which is what is being solved for. Two solutions of the linear equation no longer add — the leftover is three times the term the linear theory keeps — and the critical Mach number approaches one as the two-thirds power of thickness.

compressible · Transonic 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 roughness function, and the asymptote in which the viscosity has gone. The whole effect of a rough wall on a turbulent boundary layer is one number: the downward shift of the logarithmic profile. It vanishes on a smooth wall, rises through a transitional band, and becomes (1/kappa)ln(k+) + B − 8.5 — at which point substituting it back leaves u+ = (1/kappa)ln(y/k) + 8.5, with the fluid's own length gone from the answer entirely.

A second length at the wall

The logarithm in a turbulent wall profile exists because a region of the flow is not allowed to know about any length except the distance to the wall. Roughen the surface and there is one it does know about, which belongs neither to the fluid nor to the flow — and the slope does not change at all.

turbulence · Roughness A factor of 1.8 between the extremes. The same numbers as a bar. The spread is exactly 96 over 160/3, which is 1.8 — and it is a spread in a quantity that a single correlation using the hydraulic diameter treats as a constant. A rectangle twenty times as wide as it is deep is 78 per cent more resistant than a circle of the same hydraulic diameter.

A number that is only the shape of the hole

The friction factor times the Reynolds number for laminar flow in a duct is a pure number that depends on the cross-section's shape and on nothing else — not the fluid, not the size, not the flow rate. It is exactly 64 for a circle and exactly 96 for a slot, and the correlation that treats them as the same is wrong by a third.

applied · Internal flow 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 Four bodies, four drag coefficients, and no flow was solved. Four bodies with their Newtonian drag coefficients, each computed as a quadrature over its own surface with no flow solution anywhere. The cone's answer is exactly 2sin²δ, checked against the closed form to nine decimal places; the flat disc's is exactly 2, since every element of it faces the stream; and the sphere's is 1 against the classical Newtonian value of 1. Every one of those is an integral of one expression over a shape, and none of them required knowing what the air was doing anywhere.

The only theory simple enough to optimise

Whether Newton's sine-squared law is right has two answers — hopeless at the speeds he argued about, nearly exact behind a strong shock. This asks a different question about the same formula. Its pressure depends only on the local surface angle, so a shape's drag is a quadrature rather than a solution, and the best shape can be found by calculus.

misconceptions · Newtonian 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 A pure number with no fluid in it. The load a plane slider bearing carries, divided by everything dimensional in it. What is left is a function of the convergence ratio alone — no viscosity, no speed, no size, no material — so the tilt that carries the most load is decided before anything about the machine is known. It is the root of one transcendental equation, at 2.1887048.

The gap that carries the most

A bearing's most consequential dimension is decided by a number with no oil in it, no speed and no size — the root of one transcendental equation, at 2.1887048. And the maximum it sits on is flat enough that the value printed in every handbook is not it.

viscous · Lubrication The shaft speed is a window, and cavitation closes the top of it. Two groups against shaft speed for the same duty: the specific speed, which must be inside a band for a runner to exist, and the suction specific speed, which must be below about 3 for the impeller not to cavitate. Both rise with the shaft speed, so raising it to reach a band is also raising it towards the cavitation limit. The window here runs from 274.98 rpm to 962.31 rpm and cavitation sets its top. A duty whose window is empty needs something other than a different machine — a booster, a lower installation, or an inducer.

The group with no head in it

The number that picks a machine says nothing about whether the machine can exist. A second group formed from the same variables, with the delivered head replaced by the margin available at the inlet, decides that — and the delivered head has left the expression entirely, so how far a pump lifts is irrelevant to whether it tears the liquid apart at its own entrance.

applied · Specific speed Fourteen pure numbers, and where each came from. Every one of these is dimensionless, exact and quoted as a fact about fluids. None of them comes from dimensional analysis, which says which numbers an answer may depend on and never what any of them is. They come in four kinds — algebra, an integral, a root and an optimum — and the last two are not equally knowable.

Where a pure number comes from

Sixty-four for a round pipe, sixteen twenty-sevenths for a wind turbine, 0.332 for Blasius. Counting dimensions produces none of them — it produces the list of arguments, and the function has to be solved. Which of the four ways it was solved decides how many digits are worth printing.

regimes · Dimensional 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 A force ceiling ends the similarity: a breeze makes the boat slower as a share of the wind. Speed made good to windward as a fraction of the true wind, on the best course for each wind, for a rig that is flattened once the righting moment binds and for one that is reefed, with the hull's drag angle held at 6°. Below 3.70 m/s the two are the same boat, and the fraction does not depend on the wind — 1.122 at every speed, which is the similarity the two-angle polar rests on. Above it the flattened rig's fraction falls: 0.918 at 6 m/s, 0.572 at 10 m/s, 0.321 at 15 m/s and 0.169 at 20 m/s. The reefed rig's stays at 1.122, because a reefed sail keeps its least drag angle and the hull's angle is held. Once a force is limited, the triangle is no longer the same shape in every wind.

A breeze the boat cannot use

The two-angle polar makes a boat's speed a fixed fraction of the wind, in any wind. A righting moment ends that at 3.7 metres a second on the beat. Past it the crew must spill force, a flattened sail's drag angle climbs, and the best course to windward moves closer to the wind rather than away from it — which 45° + λ/2 cannot say, because λ now depends on the course.

applied · Sailing The shaded face's share of the force is set by K, and it is not small until K is. The fraction of a flat plate's normal force carried by its leeward face against K = M sin α, at Mach 3, 5, 10, 20, with the hypersonic small-disturbance value and the share the leeward face would have at vacuum (dashed). Newtonian theory puts it at zero. The curves collapse on K: the small-disturbance share is 35.7 per cent at K = 0.5, 24.2 at 1, 11.2 at 2 and 3.4 at 4, and the vacuum bound is 28.8, 11.4 and 3.4 per cent at 1, 2 and 4. The zero is a good approximation only where K is large — which is also the only place the Newtonian windward pressure is itself accurate.

The face Newton left in shadow

Newtonian theory gives a surface turned away from the stream a pressure coefficient of exactly zero, and at hypersonic speed the rest of the theory is nearly right. The shaded face is not. Computed exactly on a flat plate, its share of the force depends on the similarity parameter K = M sin α rather than on the Mach number, it is a quarter of the force at K = 1, and it moves a hypersonic plate's best lift-to-drag ratio from 5 to 7 at Mach 10.

misconceptions · Newtonian Eight numbers, one construction. The dimensionless groups this collection has produced, placed on a logarithmic axis at a representative value. Each is a memory time divided by a process time, each was named separately in a different field, and each decides the same question: whether the past is still present.

Every memory number is one time over another

These essays produced eight dimensionless groups in eight different fields, named after eight different people, spanning a factor of four hundred in value. Written out, all eight are a memory time divided by a process time, and all eight govern the same curve.

regimes · Dimensional 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 Three different flows, drawn on top of one another. The surface pressure of three parabolic arcs of different thickness, each at the Mach number that gives them the same transonic similarity parameter, divided by the scale the similarity rule prescribes. They lie on one curve. The three solutions were computed separately from an equation with different coefficients in each case, and nothing in the solver knows what the parameter is — the scaled peaks differ by 0.01 per cent.

One curve for every thickness

Near Mach one there is no linear theory, so there is no superposition, no scaling of a solution by the thing that caused it, and no shortcut of any kind. There is exactly one thing left — a change of variables in which the thickness and the Mach number vanish and only their combination survives — and it is worth more than the theory that was lost.

compressible · Transonic A drag that is nothing until there is something to be discontinuous about. The scaled wave drag of a parabolic arc against the similarity parameter. Above about 1.4 it is not small but zero to five decimal places — d'Alembert's paradox, which the small-disturbance equation inherits. Below it the drag rises by four orders of magnitude over a range of the parameter that corresponds to a few hundredths in Mach number, and the only thing that changed is that the flow now contains a shock.

A drag that needs a discontinuity

A closed body in a steady inviscid flow has no drag, and the transonic equation obeys that as faithfully as any other — until a shock appears inside its own answer. Then the same pressure integral that was zero to five decimal places becomes four decades larger, with no viscosity added and nothing about the body changed.

compressible · Transonic

The aeroplane rises out of its own gust

A gust pushes a wing up, and the wing, being attached to an aeroplane that is free to move, starts rising — which lowers its angle of attack and takes away part of the load the gust was delivering. How much depends on one number, and the formula every airworthiness code used for forty years is a fit to the calculation this essay does again.

circulation · Gust

One fast pump or two slow ones

A pumping station with several identical pumps in parallel has a control that a single pump does not: how many of them are running. With no static lift the answer is all of them, always, slowed together. With lift it is a number that falls with the flow and is almost never whole — the flow times the number of pumps, over the square root of the head the system asks. Staging rounds it, and the right moment to round up comes before the running pumps have run out of speed.

applied · Specific speed

Two small pumps carry two big discs

Split a pumping duty between two pumps of the same head and speed and each one keeps the big pump's impeller diameter, because the head fixes the tip speed. So the pair spins two full-size discs in their casings for the flow one disc used to carry, and the efficiency penalty that a station comparison assumed to be two points is four. That moves the lift at which splitting saves energy from 44 per cent to 65, and below a specific speed of about a half it never does.

applied · Specific speed

Moody's fifth root is true at one size

A model pump is less efficient than the full-size machine, and since 1925 the gap has been bridged with a fifth root of the size ratio. Built from its losses, a pump's inefficiency scales three ways at once — leakage as the inverse size where a clearance cannot shrink, friction as a fifth root until the passages go rough, and a fixed share not at all — so the exponent runs from a third down to a twentieth, and the fifth root is right only near one size. The same arithmetic says a duty should be split by speed, not by width.

applied · Specific speed

Named alongside it

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

DimensionlessModel limitOptimisationMeasurementDimensional analysisReynolds numberScalingSpecific speedEfficiencyAffinity lawsCorrelationMach number

All concepts