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

Wake — where it appears

The region of disturbed fluid a body leaves behind it, carrying the vorticity the body created. It is the one disturbance that does not decay with distance, which makes it the only place a drag can be seen from far away.

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

Flow past a cylinder at Re 100. 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.

When the flow lets go

Every body asks the air behind it to slow down and climb back up to the pressure it started at. Sometimes the air cannot, and the moment it refuses is separation — the source of most drag, the cause of stall, and the reason a golf ball has dimples.

viscous · Separation 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 Stokes' theorem on a solved wake at Re 40. A rectangle drawn in a viscous flow that was solved on a grid. The circulation round its boundary is computed by walking the four sides and adding up the velocity along them; the vorticity inside it is computed by adding up the stored vorticity cell by cell. The two computations share no sample point and the theorem says they must agree.

Circulation is vorticity, added up

One of these two quantities is measured by walking round a loop and one by summing over the area inside it, and a theorem says they are the same number. That reconciles the site's most confusable pair — and explains how a flow with circulation can have no spin in it anywhere.

kinematics · Vorticity 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 sides of one ball, at different pressures. The surface pressure coefficient round a ball, measured from the front stagnation point, on the side the seam trips and on the side it does not. Up to separation both follow the exact potential-flow distribution 1 − (9/4)sin²θ. After it both take the same wake pressure, which is what a manometer measures rather than what the ideal theory predicts. The asymmetry is the shaded area between them, and integrating it gives a side force of 0.2949 towards the later-separating side.

A ball that swings without spinning

A cricket ball curves in flight with no spin about any useful axis. The mechanism is not the Magnus effect; it is a seam tripping the boundary layer on one side so that side lets go later. Which way the ball then goes depends on one borrowed number, and this site's own inviscid solver supplies the value that gets it wrong.

applied · Sport ball The street, as two rows of point vortices. The exact velocity field of a staggered double row of point vortices at the stable spacing ratio, with the row spacing and the circulation of one core printed from a line integral of the field rather than from the number that built it. This is a model of a wake. It contains no body, no viscosity and no mechanism that would shed anything, and the viscous stepper used here does not produce a street at any Reynolds number.

The street this site cannot draw

The alternating wake behind a cylinder is the most photographed structure in fluid mechanics, and this site's solver does not produce one. What can honestly be drawn instead is a model of it — and the model settles one thing exactly, which is the spacing.

turbulence · Wake Quasi-steady stops being true a long way before one. The magnitude of Theodorsen's function, which is the factor a quasi-steady lift calculation is wrong by, and the phase the lift lags the motion. Quasi-steady means C = 1, and the amplitude is already one per cent low at k = 0.0061 and fifteen per cent low at k = 0.1 — a reduced frequency at which nobody hesitates to call a flow quasi-steady. The lag is worse: it reaches a degree at k = 0.003, and a flutter calculation is decided by phase rather than by amplitude.

Slow enough to be steady

A wing moving slowly enough is assumed to carry the lift its instantaneous angle asks for. The reduced frequency has two thresholds — one where the apparent-mass and circulatory lifts are equal, and one where the quasi-steady answer stops being right — and they are a hundred and seventy-eight apart.

regimes · Reduced frequency 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 The Kirchhoff flow past a flat plate. A uniform stream meeting a flat plate held across it, with two streamlines leaving the edges and never returning. Between them is a wake of fluid at rest at a constant pressure. The equations solved are the same equations that give d'Alembert's paradox for a closed body, and this flow has a drag coefficient of 0.8798.

Drag in the theory that forbids it

d'Alembert's paradox is a theorem about flows that close behind the body. Stop requiring that, let two streamlines leave the edges and never come back, and the same equations — no viscosity, no vorticity — produce a drag coefficient of 0.8798.

inviscid · Free-streamline 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 Most of the fluctuation is the switching. What a hot wire at the edge of a shear layer records, taken apart. The conditional intensity — the fluctuation inside the turbulent fluid — is flat at 0.16 by construction. The measured intensity peaks 20 per cent above it, and the excess is the third curve: the variance of a signal that keeps switching between two mean velocities, which peaks where the switching is most even and is not turbulence at all. A model calibrated against the measured curve is being fitted to an artefact of averaging.

Turbulent some of the time

At the edge of a jet or a wake a probe is inside turbulent fluid for part of the time and in perfectly smooth flow for the rest, and an ordinary time average mixes the two. Most of the fluctuation it records there is not turbulence at all — it is the switching between two states, and it peaks where the switching is most even rather than where the turbulence is strongest.

turbulence · Intermittency Elliptic loading rolls up to πb/4, and it is π that puts it there. Three loadings on the same span carrying the same lift, with the rolled-up core positions marked below each. The core sits at the centroid of the shed vorticity, which is a quadrature over the loading and needs nothing about the roll-up itself. Elliptic loading gives πb/8 from the centreline, so the pair ends up 0.7854 of the span apart — the number every wake-separation rule is written against, and one of the few places in this subject where π turns up in an answer an engineer uses directly. The bell rolls up to 0.586 and a nearly rectangular loading to 0.978, because it sheds at the tips.

Where the wake ends up

The sheet a wing sheds rolls up into two cores within a few spans, and nobody can compute the roll-up cheaply. Nobody has to: what the cores conserve is fixed before they form, and for an elliptically loaded wing the answer contains π and comes out at 78.5 per cent of the span.

circulation · Tip vortex 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 Three bluff bodies, whose Strouhal numbers agree once the wake's width is used. The measured Strouhal number of each body, and Roshko's universal number formed with the wake's width and the speed on the free streamline that bounds it. The raw numbers span a factor of 1.462; the collapsed ones span 1.0011, with a mean of 0.16281. The shedding was never body-dependent — the length in the number was.

The frequency a wake chooses

Bluff bodies shed at Strouhal numbers from 0.145 to 0.212, and the spread is not a fact about shedding — it is a fact about which length went into the number. Change the length to the wake's own width and three bodies agree to a tenth of a per cent. Then let the body move, and the number stops deciding anything at all.

regimes · Strouhal 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 Morison's two terms over one cycle, at KC = 10. The drag term, in phase with the velocity and going as its square; the inertia term, in phase with the acceleration and ninety degrees ahead of it; and their sum, which is what a load cell records. The peak of the total is not the peak of either, and its position in the cycle is the only thing in the record that says how the two are divided.

Long enough to make a wake

A Reynolds number cannot ask whether an oscillating flow gets round a body before it turns and comes back, because it has no time in it. The Keulegan–Carpenter number can, and it decides which of Morison's two terms is the force. What it discards is the phase — and a peak force measurement cannot recover it.

regimes · Keulegan–Carpenter number 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 Every unstable mode of every profile, inside one circle. The complex phase speeds of the unstable modes, scaled so that each profile's own semicircle is the unit one. Howard's theorem says every one of them must lie inside — the centre and the radius are the mean and half-range of the velocity profile and nothing else — and every one of them does, with the closest approach at 0.915 of the radius.

Every unstable wave is inside one circle

Before solving anything, you know where the answer is. Howard's theorem says the complex phase speed of any growing disturbance in a shear flow lies inside a circle fixed by the fastest and slowest parts of the profile — and by nothing else about it at all.

turbulence · Instability 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 Captured, and not captured. The amplitude of the wake's oscillation over time, for a forcing inside the capture band and one outside it. Inside, the amplitude settles; outside, it beats at the difference between the two frequencies, because the wake is keeping its own time and the forcing is keeping its.

A wake told what to do

Force a shedding wake near its own frequency and it abandons its own and adopts the forcing's. The band over which it will do so is proportional to how hard it is pushed — and near the edge of the band it takes six times as long to make up its mind.

regimes · Strouhal The wake's memory, as a gain and a phase. Theodorsen's lift deficiency against reduced frequency. It is one at zero frequency — the quasi-steady limit, where the wake has had time to convect away — and falls to a half at high frequency, with a phase lag peaking near 15 degrees in between.

The lag that makes flutter possible

This site's own flutter model set the lift deficiency to one and recorded in its notes that doing so throws away the lag which stabilises the torsion mode. Putting the lag back moves the flutter speed from 80.8 metres a second to 131, and removes the need for the structural damping that was covering the artefact.

circulation · Flutter 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 What is left behind, and what it is made of. The two rolled-up vortices behind a large aircraft, drawn to scale against its span. They sit at 78.5 per cent of the span apart, each carries 508 square metres a second of circulation, and the pair descends at 1.72 metres a second under its own induction.

A wake that says what made it

Two vortices sitting behind an aeroplane carry its weight and its span in a form that can be read back out: circulation times separation times density times speed is the lift, exactly. The reading is exact for about a minute and worth nothing after four.

circulation · Tip vortex A blade passing the vortex it shed a passage ago. The lift the blade section feels as it passes a tip vortex at five per cent of the radius. The pulse is a doublet rather than a bump — upwash on one side and downwash on the other — so the blade is pushed one way and then the other in the width of a few chords.

A blade that flies through what it shed

A rotor blade meets the tip vortex the blade in front of it left, a fifth of a second earlier, at a few per cent of the radius. What it feels is a doublet — five and a half degrees of upwash and then five and a half of downwash within a few chords — and the peak goes as one over the miss distance.

circulation · Rotor A gust's lift keeps falling where a pitching wing's stops at a half. The magnitude of Sears' function, the lift a wing gets flying through a sinusoidal gust as a fraction of the quasi-steady value, against the reduced frequency on a logarithmic axis, beside Theodorsen's function for a wing that pitches or heaves. The two agree at low frequency. Above a reduced frequency of about a tenth they part: Theodorsen's levels off at one half, because a moving wing changes its whole boundary condition at once, while Sears' keeps falling as one over the square root of 2πk, because several wavelengths of gust lie along the chord and cancel.

The gusts that cancel along the chord

A wing that pitches keeps half its circulatory lift however fast it moves. A wing flying through a gust does not: once the gust is a few chords long, its ups and downs lie along the chord together and cancel, and the lift falls without limit. For an airliner that barely touches the root-mean-square gust load, and cuts the load spectrum at the wing's own torsion frequency to a quarter of the quasi-steady value.

regimes · Reduced frequency What a downstream blade sees going past. The axial velocity a blade in the second row meets, over one revolution, as it passes through the wakes of thirty upstream blades. Each dip is one wake, and the blade meets all thirty of them every time it goes round.

A row that meets the row before it

A compressor blade is loaded and unloaded thirty times a revolution by geometry it does not have. The forcing sits at the blade count of the row in front of it and at multiples of that, and how far up the harmonics it reaches is decided by how far apart the two rows are.

circulation · Cascade

Named alongside it

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

MeasurementCirculationModel validityDragModel limitRegimeMemory kernelSeparationBoundary layerUnsteady liftControl volumeInduced drag

All concepts