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 19 essays across 7 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

Named alongside it

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

Boundary layerBlasiusCorrelationReynolds numberModel limitSeparationSimilarity solutionDimensionlessDisplacement thicknessDrag polarHeat transferInduced drag

All concepts