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

Skin friction — where it appears

The drag contributed by the shear stress the fluid exerts along a surface, as distinct from the pressure acting on it. On a streamlined body it is most of the drag, and it depends on whether the boundary layer over it is laminar or turbulent.

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

Flow past a cylinder at Re 40. A real fluid past a circular cylinder. At low Reynolds number the flow closes up behind the body much as the ideal theory says; as it rises the flow separates and a region of reversed flow appears behind, which is where drag comes from.

Everything happens in a layer you cannot see

Air has so little viscosity that ignoring it works almost everywhere. Almost everywhere leaves out a film next to the surface, perhaps a millimetre thick, and that film decides drag, stall and whether an aircraft flies at all.

viscous · Boundary layer The chart, with one exact line on it. The friction factor of a pipe against Reynolds number, for five relative roughnesses. Every curve here except one is Colebrook's correlation, solved by iteration rather than read off a chart. The exception is the short straight line at the left: f = 64/Re is the laminar solution and it is exact. The curves flatten to the right because once the roughness pokes out of the viscous layer the Reynolds number has nothing left to change.

The roughness a wall cannot feel

A rough pipe and a polished one carry the same flow for the same pressure over three decades of Reynolds number, and then suddenly they do not. What changed is not the pipe. It is the thickness of the film of fluid at the wall, which is the only part of the flow that can see the roughness at all.

applied · Internal flow A wing's two drags, and the lift at which they are equal. Friction drag and induced drag plotted against lift coefficient, with their sum above them. Friction is flat, because a surface costs the same whatever the wing is doing; induced drag rises as the square of the lift. The total is least where the two are equal, and that is also the point of best glide.

The two drags a wing pays

A wing pays for having a surface, and it pays for making lift with a finite span. One of those bills falls as it flies faster and the other rises, so there is a speed at which the total is least — and the condition for it turns out to be that the two are equal.

viscous · Drag budget 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 The cliff a rough ball reaches sooner. The drag coefficient of a sphere against Reynolds number, on log axes. The smooth curve is Morrison's correlation, which is a fit to measurements and is drawn in the colour this site reserves for a borrowed claim. The other is the same curve shifted along the Reynolds axis by a factor of 6 — a stated model of what roughness does, which is to trip the boundary layer early, and not a measurement of any real ball.

The drag that falls as it speeds up

There is a band of speeds in which a smooth ball experiences less drag the faster it goes. Not a smaller coefficient — a smaller force. Dimples move that band down to where a golf ball actually flies, and they do it by making the friction worse.

applied · Sport ball 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 cost of going turbulent

A turbulent boundary layer costs several times the friction of a laminar one, and the multiple is not a constant — it rises with Reynolds number, because the two laws have different exponents. That is why laminar flow is worth more on a long fast surface than on a short slow one.

viscous · Drag budget A boundary layer with no x in it. The velocity profile over a porous wall with uniform suction: U(1 − e^{−Vy/ν}), exactly, at every station along the wall. The displacement thickness is ν/V, the momentum thickness is half of it, and the shape factor is two — all of them constants, none of them a function of distance. It is the cleanest demonstration there is that a boundary layer's thickness is a balance rather than an accumulation.

The layer that stops growing

Blasius' boundary layer thickens as the square root of distance and never stops. Suck fluid through the wall at a uniform rate and it stops immediately — the profile becomes a single exponential with no x anywhere in it, and the friction comes out exactly equal to the momentum of the fluid that was taken away.

viscous · Exact layer 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 One layer, and the four others inside and outside it. The velocity profile over a flat plate, and the temperature profile in the same layer at four Prandtl numbers: a liquid metal at 0.01, air at 0.71, water at 7 and a heavy oil at 100. The equations differ by one number and the profiles differ by a factor of twenty in thickness. At Pr = 1 the two are the same function — not similar, identical, to eight decimal places — because the equations and the conditions are then the same, which is what every statement called a Reynolds analogy rests on.

The other layer, and the one number that separates them

A wall in a stream carries two boundary conditions and grows two layers. Their thicknesses differ by a factor of twenty across ordinary fluids, and at exactly one Prandtl number the two profiles are not similar but identical.

regimes · Peclet The air a wing carries along, and how little of it there is. The Blasius profile, in the wing's frame, with the free stream at one. No slip says the air at the surface is at rest relative to the surface, which in the ground's frame means it is moving with the wing — but only exactly at the wall. The deficit, integrated across the layer, is the displacement thickness: at a Reynolds number of 1e+6 and a metre of chord it is 1.72 millimetres of air moving at flight speed, which is the whole of what is 'carried'. The step drawn on the axis is that same deficit as a solid slab. A wing does not drag a blanket of air with it; it leaves a boundary layer behind it, and the layer is made of air that keeps being replaced.

The air a wing does not carry

No slip says the air touching a surface moves with it, and the usual reading is that a wing drags a blanket of air along. The blanket is 1.7 millimetres thick per metre of chord, it is different air every instant, and the drag it costs falls as it gets thicker.

misconceptions · The no-slip condition 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 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 Four wall models, through the buffer layer. Van Driest's damped mixing length, Reichardt's fit, Spalding's implicit law, and a control with no buffer layer at all — the viscous sublayer joined straight to the logarithm where they cross, at y+ = 11.6. The three fitted models agree with each other to a per cent and a half; the control is nineteen per cent above them at y+ = 10.

Three buffer layers, one friction

Four wall models are put through a pipe. The one with no buffer layer at all is nineteen per cent wrong where the turbulence production peaks and two and a half per cent wrong in the friction; the three respectable ones are within two per cent of each other in the buffer layer and spread by seven in the friction. The answer is not where it was expected.

turbulence · Wall law 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 4 nodes and 2 saddles, and the difference is the shape. The skin-friction pattern on a body, drawn as the streaks a film of oil would leave, looking at the nose. Every critical point of the field is marked: 4 nodes, where the streaks converge or diverge, and 2 saddles, where two streaks cross. Their indices sum to 2, which is the Euler characteristic of a sphere, and nothing about the flow enters that number — the stream is at 0 and the crossflow at 0, and neither may change it. What the flow decides is where the points are and what they look like; what the body decides is how many there are of each.

The count computed on a body

The rule for a closed surface is usually quoted and seldom solved for. This solves for one, and what the computation adds is not confirmation — it is the discovery that the count survives two events in which the number of stagnation points falls, and that both of them have closed forms.

kinematics · Topology Four guesses at a boundary-layer profile. A straight line, a parabola, Pohlhausen's cubic and a quarter sine, each rising from zero at the wall to the free stream at the edge. Two of them also satisfy the conditions the true profile satisfies — no curvature at the wall, no slope at the edge — and two do not, which is what sorts them.

Four profiles, one drag

The momentum integral is exact and asks nothing about the shape of the velocity profile. Four guesses at that shape span twenty-three per cent in the drag they give — and the two that satisfy the conditions the true profile satisfies are within three, which says the freedom belongs to the family rather than to the constraint.

viscous · Boundary layer A keel pulls the best beat 12 degrees closer to the wind. Speed made good as a fraction of the wind against the course, in 3 m/s: for the boat with its keel solved as a wing, and for a boat whose λ is fixed at 25.69°, the value the keeled boat has on its own best course. The search puts the keeled boat's best beat at 45.83°; the fixed-angle rule puts it at 45° + λ/2 = 57.84°. Pointing higher loads the keel towards its best lift coefficient — 0.130 at 40° against 0.076 at 60° — and lowers its drag angle, so the curve peaks early and falls away faster on the far side. The fixed-angle curve promises 0.653 of the wind; the keeled boat can make 0.553, and only by sailing twelve degrees higher than the rule says.

A keel flies wherever the course puts it

A keel is a wing whose lift is whatever side force the rig happens to make, so its lift coefficient is chosen by the course and the wind rather than by its designer. Solved that way, the hull's drag angle stops being a property of the boat: it pulls the best beat twelve degrees closer to the wind, halves the reaching speed a fixed angle predicts, and finally charges for reefing.

applied · Sailing A cooled wall at β = 1: the linear relation is 158 K out inside the layer. The static temperature across a laminar layer at Mach 5 with an edge temperature of 220 K, over a wall whose total enthalpy is 0.5 of the edge's, at a pressure gradient β = 1: exact (thick) and from the Crocco–Busemann linear relation (thin), at a Prandtl number of one with constant properties. The wall is at 660 K in both. The exact profile peaks at 684 K and the linear one at 759 K; the largest difference, −158.3 K, is at η = 0.81.

The gradient the heat never hears

On a flat plate at a Prandtl number of one, a boundary layer's total enthalpy is a straight-line function of its velocity, whatever the wall's temperature. Put the same layer in a pressure gradient and the straight line fails everywhere except on an insulated wall, because the gradient enters the velocity's equation and not the enthalpy's — and a favourable gradient can leave a band of gas colder than the free stream above a wall three times hotter than it.

compressible · Recovery Three profiles that do not depend on the radius. The radial, azimuthal and axial velocities of the flow above a rotating disc, as functions of one similarity variable. The radial one is a jet: fluid thrown outward by the swirl it has picked up, peaking at 0.181 of the local disc speed a fifth of the way through the layer. The azimuthal one falls from the disc's own speed to nothing. And the axial one is the surprise — it does not vanish far from the disc but tends to a constant, so the disc draws fluid down onto itself at 0.8845 times the square root of the viscosity times the rotation rate, at every radius and for ever.

The solution that keeps its nonlinear term

Every exact solution before this one has been exact because the nonlinear term vanished. A rotating disc's does not vanish — at the wall it is the whole of the balance — and the reduction is exact anyway, because the radius divides out of all three momentum equations at once.

viscous · Exact layer One bowl, and a line across it for every way of flying. The least induced drag of a wing, a canard and a tail together, as a function of the canard's and the tail's shares of the lift, drawn as rings of equal drag about the bowl's foot, where each carries 1.9 per cent. A static margin of a tenth of a chord and a wing pitching moment draw a straight line of trimmed splits across the bowl. The tail aircraft trims where its line crosses the axis of zero canard load; the three-surface aircraft slides along the same line to the point nearest the foot.

A third surface is worth a square

A canard and a tail each leave an aircraft one free share of its lift, and the static margin spends it. Give an aircraft both and one share stays free after trim, so the split can slide along a line to the point of least induced drag. What that slide is worth turns out to be a square: half the bowl's curvature times the square of how far the tail aircraft already trims from twice the bowl's foot. For a cruising wing the saving is under one per cent and less than the canard's own skin friction; it pays only for a wing whose pitching moment is as large as a flapped section's.

circulation · Moment Holding the neutral point with a canard ahead takes a tail two and a half times the size. The tail area, as a fraction of the wing's, that puts the neutral point where the tail aircraft with a tail of 0.2 has it, against the tail's height above the wing's plane in spans: with no canard, with the canard, and with the canard but its wake at the tail removed. At the reference height the tail aircraft needs 0.2; with the canard it needs 0.496, of which 0.22 is for the canard's own lift ahead of the centre of gravity and 0.077 for its wake at the tail. A tail 0.25 span up needs 0.354 with the canard and 0.154 without.

The tail a canard needs costs more than the canard saves

A canard added to an aircraft as a second trimming surface saves induced drag only if the neutral point is held where it was, and holding it is not free. The canard's own lift ahead of the centre of gravity pulls the neutral point forward, its wake reaches the tail and weakens it, and the tail has to grow to put the neutral point back. For a canard of a tenth of the wing's area the tail grows two and a half times, and its skin friction is twenty to thirty counts against a saving of one to twenty. A T-tail escapes most of the wake and a third of the bill, and still only a flapped wing comes near to paying.

circulation · Moment A cold wall holds its layer on longer; a hot one lets it go sooner. The pressure-gradient parameter β at which the wall shear vanishes, against the wall's total enthalpy over the edge's. With the velocity deaf to the temperature it is −0.1988 for every wall. Coupled, it is -0.3264 for a wall at absolute zero, -0.2623 for a wall at half the edge's total enthalpy, -0.1988 at one, -0.1573 at one and a half and -0.1295 at two. Cooling delays separation by more than the same heating hastens it.

A cold wall holds its layer on longer

A compressible boundary layer's pressure gradient pushes on its density, and its density is its temperature. Solve the velocity and the enthalpy together and the wall's temperature reaches the velocity: a cold wall holds its layer on to an adverse gradient two-thirds steeper, a hot one lets it go a third sooner, the stagnation point's heat-to-friction ratio more than doubles between them, and the cold band an uncoupled layer put above a cooled nose disappears.

compressible · Recovery Above the logarithm the profile lifts away: the wake. The mean velocity in wall units against the distance from the wall, for a flat-plate layer at Reθ = 10⁴ (δ⁺ = 3484), from Spalding's inner law alone and with Coles's wake of strength Π = 0.3, 0.55 and 1. Below about a fifth of the layer the curves coincide on the logarithm; above it the wake lifts the profile by up to 2Π/κ, 2.68 wall units for a flat plate — a tenth of the edge velocity, and the part of the profile a log law cannot describe.

The wake is a tenth of the velocity and a third of the displacement

The logarithmic law describes a band in the middle of a turbulent boundary layer, and above it the profile lifts away by an amount Coles called the wake. It is a tenth of the edge velocity, so it looks like a correction. It is not: it carries a third of the layer's displacement, it is what turns the log law into a friction law for a boundary layer, and with it the friction comes out within a few per cent of a measured correlation that contains no logarithm at all.

turbulence · Wall law An adverse gradient hollows the layer from the outside in. The velocity profile of an equilibrium layer, as a fraction of the edge velocity against y/δ, at Clauser's β = 0, 2, 10 and 50, with the wake strength from Das's fit. The inner part keeps the wall law; the outer part's wake grows — Π = 0.476, 1.59, 4.67 and 15.2 — and the profile sags until most of the layer is moving slowly over a thin fast wall region.

An adverse gradient grows the wake, not the log law

Push a turbulent boundary layer up a rising pressure and its profile changes from the outside in. The wall region keeps its law; the wake grows, the layer hollows, the shape factor climbs and the friction falls. Clauser found the layers that stay similar while this happens, and their arithmetic has a surprise in it: there is a steepest pressure rise any such layer can climb, close to a free stream falling as the inverse fourth root of the distance, and the friction approaches zero only as the layer becomes all wake.

turbulence · Wall law

A layer released from a gradient overshoots the shape it left

A turbulent layer that meets an adverse gradient suddenly entrains only half what the gradient asks, because its large eddies were made upstream, and it separates sooner than a layer that answers at once. Release it — let the outer flow stop decelerating, as at a diffuser's exit or a wing's trailing edge — and the same lag runs the other way. The eddies made in the gradient keep entraining after it has gone. The layer arrives at the release with the worse shape and, twelve thicknesses later, has the better one, its shape factor falling past the flat plate's and its friction rising past the shape-only layer's.

turbulence · Wall law

A stagnation point can leave the flow through the wall, and the count is kept by halves

A no-slip wall is a sheet of zeros, so the index count that holds for a flow in open space has nothing to count there. Divide the velocity by the height and the wall's zeros shrink to the skin friction's nodes, foci and saddles; reflect the result through the wall and each of them is counted once while every point in the fluid is counted twice. That is why the points on a wall are halves, and an exact swirling flow onto a plate shows the halves doing their job: a stagnation point in the fluid can disappear without a partner by sinking into the wall, where the skin friction's focus changes from attachment to separation to pay for it. A section through the flow sees none of this, because in the section every one of the points is a saddle.

kinematics · Topology

A stagnation point off the axis can land only on a saddle or a node, and turns it inside out

In a swirling flow onto a plate every stagnation point lived on the axis, so one sinking to the wall always met the wall's single critical point. Without the symmetry it might seem free to land anywhere — on an ordinary stretch of skin-friction line, leaving a pair of new points behind. It cannot. Divided by the height, the velocity is continuous down to the wall, where it becomes the skin friction, so a point can only reach the wall where the skin friction already vanishes. It lands on a saddle or a node, and the count does the rest: the point it lands on turns from separation to attachment. An oil-flow picture of the wall, which shows the same pattern before and after, cannot tell which it is looking at.

kinematics · Topology

Named alongside it

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

Boundary layerModel limitSeparationBlasiusShape factorCorrelationInduced dragReynolds numberSimilarity solutionAdverse pressure gradientDisplacement thicknessHeat transfer

All concepts