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

Drag — where it appears

The force on a body along the direction of the flow, made up of a pressure contribution and a friction one. Exact inviscid theory says it is zero, every measurement says otherwise, and the gap between the two is where most of this subject lives.

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

d'Alembert's paradox, measured. Surface pressure round a cylinder in ideal flow, plotted against angle. The distribution is symmetric front to back, so every push on the front is matched by an equal push on the back, and the total force along the stream is exactly zero.

The exact theory says nothing has any drag

Solve the flow past a body in a fluid with no viscosity and the answer is beautiful, closed-form, and predicts that a cyclist needs no legs and an airliner no engines. This is not a small error, and it is the most useful failure in the subject.

inviscid · Dalembert Creeping flow, and the same body with inertia. The exact creeping-flow solution beside a solved field at a Reynolds number where inertia matters. The creeping flow is a mirror image of itself front to back — a photograph of it run backwards is a photograph of it — and the field with inertia has a wake, which is what a direction of time looks like.

The world with no inertia

Drop the viscosity and the equations become exactly solvable and wrong about drag. Drop the inertia instead and they become exactly solvable again — and for a cylinder in an unbounded fluid there is no solution at all, which took fifty years to notice and longer to fix.

regimes · Reynolds The speed for least drag is not the speed for least power. Drag and power against airspeed, in units of the speed at which drag is least. Power is drag times speed, so its minimum sits slower — at the fourth root of a third of the least-drag speed, which is 0.76 of it. Flying for range and flying for endurance are therefore different speeds, and the difference is not a rule of thumb.

The cheapest way to stay up

There is a speed at which an aircraft's drag is least, and a different, slower speed at which its power is least. The ratio between them is the fourth root of a third — a number that does not depend on the aircraft, the altitude, or anything else about the flight.

viscous · Drag budget Three times the wind, on a reach. The polar diagram: boat speed in every direction, as a multiple of the true wind speed, for drag angles of 14 and 6 degrees. The shaded wedge at the top is the no-go zone, whose half-angle is exactly the sum of the two drag angles. Everywhere outside about twice that angle the boat is faster than the wind, and the maximum is 2.92 times the wind at 110 degrees — which is 1/sin λ at 90° + λ, both checked.

Faster than the wind that drives it

An ice yacht in a fifteen-knot breeze does forty. That is not a trick and it does not need a special sail — it follows from two drag angles and a triangle, and the best speed a boat can reach is one over the sine of their sum.

applied · Sailing How far downstream the heat is still being made. The dissipation accumulated from a station ahead of a cylinder to a station behind it, as a fraction of the whole of what is made inside the frame, at five Reynolds numbers. At Reynolds number 1 the fluid has finished paying by about a diameter behind the body. At 100 it has not finished at five, and the curve is still climbing at the edge of the picture — the drag is a force on the body, and the heat it stands for is somewhere else.

Where the heat of a drag is made

The power it takes to tow a body through a fluid becomes heat, all of it, eventually. None of the interesting words in that sentence are the first four. It is the "eventually" that decides how a wake behaves, how far a disturbance reaches, and why no box drawn round a body contains its own bill.

viscous · Dissipation Two things moving by a million, and their product standing still. Viscosity, the squared velocity gradient at the dissipation scale, and their product, across six decades of Reynolds number at a fixed large-scale flow. The viscosity falls by a factor of 1e+6; the squared gradient rises by exactly the same factor, because η falls as Re^(−3/4) and u_η as Re^(−1/4); and the dissipation ν(u_η/η)² does not move at all. That is the dissipation anomaly stated as arithmetic: the limit of the dissipation as viscosity vanishes is not the value it takes when viscosity is zero.

The limit that is not the value

Dissipation is viscosity times the square of a velocity gradient, so it ought to vanish as the viscosity does. It does not. The gradient rises by exactly the factor the viscosity falls by, the product stands still, and a fluid with no viscosity at all would dissipate nothing — which is why the limit and the value are different numbers.

turbulence · Dissipation A sphere does not go at the speed of the flow it is in. How far a force-free sphere on the axis of a round tube lags the fluid at its own centre, against its size. The lag is two-thirds of the square of the size ratio and comes from the reciprocal theorem in one line: the sphere moves at the average of the ambient flow over its own surface, and a parabolic profile is slower everywhere on that surface than at the middle. At a fifth of the tube's radius the lag is 2.7 per cent, which is the difference between a tracer and a thing being measured.

A force without the flow that makes it

A sphere carried along by a flow does not travel at the speed of the fluid at its centre. It travels at the average of the flow over its own surface — and getting that result needs no solution of the flow around the sphere at all, only the answer to a completely different problem that everybody already knows.

viscous · Reciprocity 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 Two-thirds to one, and nothing outside it. The drag on a spherical drop as a fraction of Stokes' drag on a rigid sphere of the same size, against the ratio of the viscosities. It runs from exactly two-thirds for a clean gas bubble to exactly one for a rigid particle and takes no value outside that range, however light or heavy the drop is — the formula has no density in it anywhere. A drop of water in air and a drop of mercury in water sit at opposite ends of it.

The surface that moves with the flow

A clean gas bubble feels two-thirds of the drag a rigid sphere of the same size would, and the formula has no density in it anywhere. What buys the third is that the bubble's surface is free to move — and real bubbles in ordinary water do not get it, for a reason that is a millionth of a per cent of the water by mass.

viscous · Mobile interface A fifth of the drag power is not heat yet. The energy account of towing a flat plate a metre long through air at thirty metres a second. The power it takes is the drag times the speed. The heat made inside the boundary layer is half the free-stream energy times the energy thickness, and it is less — 78.6 per cent of what was paid. The rest has not been destroyed: it is kinetic energy still in the wake, which will become heat somewhere downstream, in fluid that is no longer touching the plate.

The third thickness

A boundary layer has no edge, so every thickness quoted for it is an integral of the profile against some weight. Two of them are famous. The third answers a question the other two cannot — how much of the power spent towing a plate has actually become heat by the time the fluid leaves it — and the answer is 78.6 per cent.

viscous · Energy thickness 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 integrands, and the wrong one claims 9.58 per cent more drag. The two things that get integrated across a wake, each scaled to its own peak so the shapes can be compared. The momentum integrand u(U − u)/U² is the drag; the mass integrand (U − u)/U is the displacement thickness, and it is not a drag at all. They differ by a factor of u/U inside them, so the mass one is fatter wherever the deficit is deep — and its integral here is 1.1 times the momentum one's. That ratio is decided by how deep the wake is rather than by how wide: at a twentieth of this momentum thickness it falls to 1, and at twice it rises to 1.25. The error is smallest exactly where a survey is properly done, far downstream where the wake has spread and shallowed.

Weighing what is missing

A control volume drawn round a wing gets the lift out of it, on a flow that was solved exactly. The wake survey asks the same box for the drag on a flow nobody has solved, and it is how a real aerofoil's drag is known — with three assumptions, all of which are checkable, and one integral standing next to it that is wrong.

misconceptions · Momentum lift The only candidate the far field allows, and the wall it slips past. The general Stokes solution has four constants; the condition at infinity kills two of them and fixes a third, leaving one to satisfy two conditions at the wall. Setting the stream function to zero there uses it up, and the tangential velocity that remains is exactly twice the free stream — for every radius, every speed, and every fluid.

The flow with no solution

Creeping flow past a sphere has a solution and everybody knows it. Creeping flow past a cylinder has none — not a difficult one, not one needing a clever method. The equations, the no-slip condition and the uniform stream at infinity are inconsistent, and the residual is exactly twice the free stream.

viscous · Stokes' paradox Spin against distance flown, not against time. The spin of a struck golf ball as a fraction of its launch spin, against how far it has flown. The integrated flight sits exactly on an exponential in distance, because the spin-down torque is proportional to speed times spin and the speed cancels. The constant is 1,202 metres and the drive is 200, so the ball arrives with five sixths of the spin it left with.

The ball that never forgets its spin

Commentary explains a late-swerving ball by saying the spin is dying away. The spin-down torque goes as speed times spin, so spin decays over a distance rather than a time — 1,202 metres for a golf ball against a 200-metre drive — and what dies away is the speed, which makes the swerve stronger.

applied · Sport ball The free surface a submerged body leaves behind it, and the flat water in front. The linearised free-surface problem solved as a Fourier integral with a radiation condition. Behind the body a wave train of the wavelength that stands still relative to it, 2 pi U²/g; ahead of it, an amplitude a hundred and twenty times smaller. The asymmetry is the drag: an ideal fluid with a free surface can carry energy away.

The drag that is made of waves

D'Alembert's paradox says a body in a steady, irrotational, incompressible, inviscid flow feels no drag. Put a free surface above it and every one of those words still holds — and the drag is not zero. It is the energy walking away in the wave train behind.

inviscid · Dalembert Four wakes carrying exactly the same drag. Four velocity-deficit profiles behind a body, each normalised so that the integral of the deficit across the wake is exactly the same. That integral is the drag: the far-wake momentum balance says so with no assumption about the shape of anything. The four are a narrow Gaussian, a wide top hat, the two-lobed wake a body with a splitter plate leaves, and a profile with heavy tails.

Exact in the total, free in the profile

A constraint is one number imposed on a function. Four wakes built to carry exactly the same drag differ by a factor of four and a half in their peak deficit, and the general statement behind that is a question about angles: how much of the wanted answer survives being projected off the constraints, and how much does not.

regimes · Dimensional Two forces on one cylinder, a quarter of a cycle apart. The inertia and drag terms of Morison's equation over one wave period, at a Keulegan-Carpenter number of ten. The inertia term follows the acceleration and peaks where the velocity is zero; the drag term follows the velocity and peaks where the acceleration is. They are a quarter of a cycle apart and they are different kinds of quantity.

Two forces, and only one of them remembers

Morison's equation adds an inertia term to a drag term and is usually presented as an empirical patch. It is not: the two terms are the two kinds of memory this collection has been separating, one a function of the present acceleration and one a function of the wake left by the previous half cycle.

regimes · Keulegan–Carpenter number Two wakes with one drag. The two initial velocity deficits: a slab, which is roughly what a bluff body leaves, and a pair of separated lobes, which is roughly what a body with a jet through the middle of it leaves. Their integrals are identical, so the two bodies have exactly the same drag.

A wake that keeps the drag and forgets the body

Two very different wakes with the same momentum deficit converge to the same profile, because the deficit is conserved and everything else diffuses away. The convergence is a power law rather than an exponential, so it takes two hundred widths for twenty per cent agreement and nine hundred for five.

turbulence · Wake Two external flows that agree where it matters and nowhere else. The velocity just outside the boundary layer, for two pressure distributions, against distance along the surface. They cross at the half-way station with the same value and the same gradient, and they have nothing else in common: one accelerates steadily and the other does most of its accelerating at once.

A layer that is an integral of everything upstream

Two surfaces are given external velocity distributions that agree exactly at one station — the same speed and the same gradient. The boundary layers there differ by 38 per cent in momentum thickness, and the two surfaces separate five per cent of their length apart.

viscous · Boundary layer 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 The trailing sheet rolls up into two vortices, and nothing it carries is lost. The trailing vortex sheet behind an elliptically loaded wing, seen in a plane across the wake, at times 0, 0.05, 0.2, 0.6 in units of b²/Γ₀, represented by 160 point vortices with a smoothing length of 0.03 of the span. The tips curl up first and the sheet winds into two concentrated vortices while the whole system sinks under its own induced velocity; by t = 0.6 the pair's centroid has descended 0.122 of the span. Through all of it the crossflow energy — the induced drag — and the separation of the two halves' centroids, 0.7854 of the span, stay exactly what they were.

The drag a wake keeps however it rolls up

A plane drawn across the wake of a finite wing contains its induced drag as the kinetic energy of the swirling crossflow. The trailing sheet then rolls up into two vortices, and the energy does not change at all — roll-up moves the drag around the plane without spending any of it. What does spend it is viscosity, which turns crossflow energy into a total-pressure defect, so a plane farther back reads less induced drag, more profile drag, and the same total.

misconceptions · Momentum lift

Named alongside it

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

MeasurementModel limitWakeControl volumeCreeping flowModel validityDissipationMemory kernelReynolds numberBoundary conditionBoundary layerd'Alembert's paradox

All concepts