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

Optimisation — where it appears

The search for the shape or the distribution that makes a stated quantity least, subject to a stated constraint. In this collection the answers are usually exact — an elliptic loading, a Betz limit, a least-cost radius — and the constraint is where the physics is.

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

The fastest course is never the one towards the mark. Boat speed and speed made good against the course sailed. Speed rises steadily as the boat bears away, but what counts is the component along the direction wanted, and that has a maximum well off the straight line — at 55.0 degrees going upwind and 145.0 going down. Both optima are found here by search and agree with 45° + λ/2 and 135° + λ/2 to six figures.

The fastest way is not the straight one

A boat racing to a mark dead upwind sails seventy per cent further than the distance to it, and arrives first. The best angle is forty-five degrees plus half the apparent wind angle, it comes out of one line of trigonometry, and the same argument says to gybe downwind rather than run.

applied · Sailing Two bills, and the radius that settles them. The cost of a vessel against its radius: the pumping power, which falls as the inverse fourth power, and the price of owning the fluid and the wall, which rises as the square. Their sum has a minimum, found here by golden-section search and agreeing with the closed form to eight figures. At that radius the pumping bill is exactly a third of the total — for every set of constants, because it follows from the two exponents alone.

The radius that costs least

A vessel that carries a flow costs two things to own — the power to push fluid along it and the price of the tissue itself. Minimising the sum gives a best radius, the best radius makes flow proportional to radius cubed, and the rule that follows is a statement about a photograph that came out of a cost function.

applied · Branching 20.45 MPa out of a film 25 µm thick. The pressure along a tapered pad, from the closed-form solution of Reynolds' equation, with the same equation's tridiagonal grid solve drawn over it as points. The peak is 20.45 MPa — enough to yield mild steel — and it sits at 69 per cent of the way along rather than in the middle, because the pressure gradient vanishes where the film equals the harmonic mean of its two ends and the harmonic mean is biased towards the thinner one. The pad's own shape is drawn along the top, to a vertical scale of its own. Nothing pumps this oil: the runner drags it into a narrowing gap and the gap does the rest.

Nothing but the shape of the gap

A machine that holds a steel shaft off its bearing with a film of oil twenty-five microns thick has no pump in it, and the pressure it generates would yield mild steel. The mechanism is not the oil and not the speed; it is that the gap narrows.

viscous · Lubrication Most power at exactly half the jet speed, found by search at 0.5000. The power a bucket takes from a jet, against how fast the bucket runs, for four deflection angles. Every curve is a parabola with roots at zero — where the force is greatest and the bucket is not moving — and at the jet speed, where the bucket is running away and there is no force at all. The peak is halfway between, at U = V/2, and it is there for every angle and every flow rate. A golden-section search that knows none of the algebra puts it at 0.500000.

Half the jet speed takes everything

A bucket standing still feels the largest force and does no work; a bucket running with the jet does no work either. Between them the power peaks at exactly half the jet speed, for every bucket shape and every flow rate — and at that speed a perfect bucket leaves the water motionless.

applied · Jet Five quarters of the span, the same bending moment, 64/75 of the drag. The bell-loaded wing's induced drag and root bending moment against its span, both as fractions of an elliptic wing of unit span carrying the same lift. At equal span the bell is worse: its span efficiency is exactly three quarters. But its bending moment is exactly four fifths, and bending moment grows in proportion to span while drag falls as its square — so there is a span at which the bell has bought back the structure and is ahead on drag. It is at exactly five quarters, where the two curves are at 1 and 0.8533. Both numbers are rational and neither was put in by hand.

The loading nobody used

Elliptic loading is the least-drag answer to a question no aeroplane asks. Constrain the moment the lift makes about the wing root instead of the span, and a different curve comes out — five quarters of the span for sixty-four seventy-fifths of the drag, both exact — with an upwash over the outer wing and a yaw that turns the right way.

circulation · Bell shape The wake, and the velocity it gives itself. The cross-section of the wake far behind a wing with turned-up tips, with the velocity the wake induces on itself drawn as arrows normal to the trace. Induced drag is the integral of the circulation against that velocity and nothing else — this picture contains the entire quantity. It also contains no information whatever about where the surfaces were: two wings a chord apart and two wings ten chords apart produce the same picture and therefore the same drag, which is Munk's stagger theorem stated as a fact about what the arithmetic can see.

A wing that leaves the plane

A winglet is not a fence and it does not block anything escaping round the tip. The induced drag of any system of lifting surfaces depends on one cross-section of its wake and on nothing else whatever, and a wake that reaches upwards is cheaper for the same reason a wake that reaches sideways is.

circulation · Winglet Nothing can deliver more than h/H, and here that is 10.0%. The fraction of the supply a ram can deliver, against the height it is asked to deliver to, from a supply falling 2 m. The upper curve is the exact ceiling h/H, which follows from the energy audit with every loss set to zero and can be reached by no real machine; the lower one is what a ram at 65% efficiency actually sends. Asking for twice the height halves the delivery, exactly, and there is no design that escapes it.

A pump with no engine

A hydraulic ram lifts water uphill using nothing but the water that is already falling. It has one moving part and no power supply, and everything it can and cannot do follows from an energy audit that fits on one line — including a ceiling nothing about its design can move.

applied · Water hammer The parabola is the cheapest shape the walls allow. The dissipation of the profile u = A(1 − |y/h|ⁿ) carrying a fixed flux between fixed walls, against the exponent n, divided by the parabola's. Every member of the family satisfies the same boundary conditions and carries the same fluid; they differ only in shape. The minimum is at n = 2 exactly, which is not a coincidence — it is Helmholtz's theorem, and the parabola is a solution of the equations because it is the least dissipative shape rather than the other way round.

The cheapest shape the walls allow

The parabola in a pipe is usually presented as what the equations give. It is better understood the other way round — of every profile that could carry that flow between those walls, it is the one that destroys the least energy, and the equations give it for that reason.

viscous · Minimum-dissipation A pump curve out of the momentum theorem. The pressure an ejector delivers, against how much it is entraining, at a fixed nozzle. It has the shape of every pump characteristic ever measured — a shut-off pressure with no flow, falling to no pressure at free delivery — and it was obtained from a momentum balance on a tube with nothing in it. The shut-off value here is 36.0 kPa and the machine at its best power runs at 4.56 times its own motive flow.

Mixing is a pump

Two streams at different speeds mixing in a tube destroy energy — the same Borda–Carnot expression a handbook prints beside a sudden enlargement, with two streams in it instead of one. And while they destroy it the pressure rises, which makes the loss the mechanism of a machine with no moving parts.

applied · Ejector Tip to tip, two wings cost exactly half of what they cost apart. The induced-drag saving of a pair of wings against the gap between their tips, as a fraction of what the two pay flying alone. With the tips touching the pair is one wing of twice the span carrying twice the lift, and the arithmetic of that is exact: the drag halves, and the computation returns 0.499257 of the separate figure. Pull them apart and the saving falls away with the square of the distance. Birds fly in a V because the tips are the part worth overlapping, and the spacing that pays is a small fraction of a span rather than any distance a formation could hold by eye.

The lift beside a wing

Fly two aeroplanes with their wingtips touching and the pair costs exactly half what the two cost apart. Not approximately half — the arithmetic is a closed form, because two wings tip to tip are one wing of twice the span, and induced drag goes as the square of it.

circulation · Formation 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 A closed wake, and the loading of least drag on it. The wake of a box wing in the plane that decides its drag, with the circulation of least induced drag drawn as a thickness along it. The horizontal members carry a loading close to elliptic and the vertical ones carry a share that lifts nothing — they contribute no lift, since lift is Γ dy and dy is zero on a vertical, and they change the drag by changing where the wake's vorticity is. At a gap of 20 per cent of span this system costs 67.1 per cent of what a single wing of the same span and lift would.

A wake that closes on itself

Prandtl's best wing system is a rectangle, not a wing. Solving for the loading on a closed wake turns up a circulation that costs nothing and does nothing — a gauge freedom in the middle of an optimisation — and a drag that keeps falling with no floor under it.

circulation · Box wing Four bodies of identical length and volume, and their wave drags. Each body has the same length and the same volume; only the distribution of area along it differs. The Sears–Haack body — the spindle whose area goes as the three-halves power of x(L−x) — has the least wave drag of the four, and every other shape pays between thirty-seven per cent and a hundred and seventy-five per cent more for carrying the same volume the same distance. Nothing about the cross-sections' shape enters: only the area distribution does.

The least drag a volume can have

A body's supersonic wave drag depends on nothing about it except how its cross-sectional area is distributed along its length. Minimising that for a given volume gives one shape — and the answer goes as the volume squared over the fourth power of the length.

compressible · Wave drag So splitting a turn into N ramps costs one over N squared. A twelve-degree compression at Mach 3, done in one ramp and in up to sixty-four. The entropy is N times a cube of one Nth, so it falls as exactly the inverse square of the number of ramps — the measured exponent is −2.00 — and sixty-four ramps cost a two-hundred-and-fifty-sixth of what one costs.

A compression that costs nothing in the end

Turning a supersonic stream away from itself is free and turning it into itself is not. But the price of a compression is the cube of its strength, so splitting one turn into N turns costs one over N squared — and in the limit the compression is free too.

compressible · Expansion fan 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 Same disc, same solidity, different number of blades, different answer. Thrust and power coefficients for six blade counts at a fixed solidity. They span a factor of 1.101 in thrust, and every one of those rotors is the same actuator disc: same area, same blade area, same tip-speed ratio. The disc theory cannot distinguish them because the blade count is not one of its variables.

A disc that knows no blades

Momentum theory replaces a rotor with a surface across which the pressure jumps, and gets the Betz limit, the induced velocity and the whole energy argument out of it. It has no chord, no section and no number of blades — and at one fixed solidity, two blades and twenty give thrust coefficients ten per cent apart.

circulation · Blade-element How the weight is shared between the two surfaces, against the centre of gravity. The wing's load and the tail's, for an aeroplane in level flight, against the static margin. The two must sum to the weight and their moments must cancel, and those two equations decide the split. At a forward centre of gravity the tail carries down and the wing carries more than the weight; the crossing is where the tail carries nothing.

More lift than weight

An aeroplane in level flight is drawn with one arrow up and one arrow down, equal and opposite. That equation collapses the whole configuration onto one number, and it is not true of any aeroplane with a tail behind it: there are two surfaces, two equations, and the second decides the split.

circulation · Trim 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 The excess energy is the energy of the difference, exactly. Add any admissible perturbation to the potential flow and its kinetic energy rises by precisely the energy of the perturbation itself — not approximately, and not to leading order. The measured excess and the perturbation's own energy lie on one another to two parts in 10¹¹ across a sixty-fold range of amplitude.

How much more than the least

Of all the flows that conserve mass and stay inside the walls, the ideal one carries the least energy. That is a theorem, and the useful half of it is the part nobody quotes: it says by exactly how much every other flow misses, and the answer is the square of how wrong it looks.

inviscid · Least energy How flat the optimum is. The same curve near its maximum, with the bands within a tenth and a half of one per cent of the best shaded. Every taper ratio from 0.31 to 0.42 is within a tenth of a per cent of the optimum, and the whole band from 0.25 to 0.51 is within half a per cent. The optimum is exact, and choosing it rather than its neighbour buys nothing a wing can measure.

The optimum that does not matter

Elliptic loading gives the least induced drag, exactly, and the proof is one line: the penalty is a sum of squares. Which is also why the optimum is flat enough that a quarter of the design space sits inside a tenth of a per cent of it.

circulation · Taper 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 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 Six nozzles, six pump curves, and where each one is best. The head ratio a water jet pump delivers against the flow ratio it entrains, for six area ratios from a narrow nozzle to one filling four-fifths of the throat. A wide nozzle makes a tall, steep curve that is finished at a small flow; a narrow one makes a low, long curve. The dots are each curve's best-efficiency point. Nozzle, suction, throat-friction and diffuser losses are included at borrowed representative values, and the mixing loss is computed.

The nozzle that is best at one thing

Put the four losses back into a jet pump and three questions get three answers. The most head comes from a nozzle four-fifths of its throat, in closed form; the best efficiency from one a quarter of it; and the most flow from whichever nozzle is smallest, because flow has no optimum at all.

applied · Ejector Flat until the back pressure reaches a value, then gone. Entrainment ratio against back pressure, as a multiple of the suction pressure, for three mixing-section sizes of one steam ejector driven from a motive supply a hundred times its suction pressure. Each is exactly flat while the entrained stream is choked beside the jet, up to its critical back pressure. The dashed lines join that point to the back pressure at which the entrainment has fallen to nothing, 1.4 per cent higher for the middle machine. The model fixes those two ends and not the path between them. A larger mixing section entrains more and breaks at a lower back pressure.

A choke that belongs to two streams

A steam ejector entrains a fixed amount of gas whatever its back pressure, up to a pressure where it stops. The flat part is a choke, and the entrained gas is not at Mach one when it happens: it is at 0.886, because the supersonic jet beside it is part of the same throat. One area ratio then trades that entrainment for compression, and the trade decides how a vacuum train is built.

applied · Ejector The best efficiency runs from nine-eighths of φ² to one. The best efficiency a peristaltic pump can reach, against the fraction of the channel its wave closes, with its two limits. For a shallow wave it is 9φ²/8, which is small — a wave closing a fifth of the channel is at best 4.5 per cent efficient. As the wave closes the channel the best efficiency tends to one, and its shortfall shrinks in proportion to the remaining gap: about 1.9(1 − φ). Nothing in between is independent of the amplitude.

The pump that is better the more it squeezes

A waving sheet swims at a cost per metre with no amplitude in it. A waving wall pumping fluid is the same mechanism turned round, and its efficiency is nothing like that: it starts at nine-eighths of the amplitude ratio squared, is exactly 2 − √3 at half closure, and rises towards one as the wave closes the tube — where the pump stops being a wave and becomes a piston.

kinematics · Continuity What every wing pays to bend its root less. Least induced drag against the root bending moment of the lift, both as fractions of the elliptic monoplane of the same span and lift, for the monoplane and for box wings with gaps of a tenth, a fifth and two-fifths of the span. Each curve is a parabola with its minimum at that wing's unconstrained optimum. The box wings' minima sit to the right of the monoplane's — they load their lift further out — and their parabolas are shallower: the gap-of-a-fifth box can bend its root 21 per cent less than the elliptic monoplane and still match its drag.

A lighter spar turns a box wing into a biplane

Prandtl's best wing system has two-thirds of a monoplane's induced drag at a gap of a fifth of the span, and part of that saving is carried by circulation turning the corner into its fins. Ask the box to bend its root less and it pays about half what a monoplane pays — but a biplane with no fins pays nearly as little, and by the time the spar is a fifth lighter the fins carry almost nothing and the box has become the biplane it was built from.

circulation · Box wing

Named alongside it

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

Model limitEfficiencySimilarityDimensionlessInduced dragMeasurementCirculationSpan loadingThe Trefftz planeConstraintCorrelationDownwash

All concepts