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

Discretisation — where it appears

The replacement of a continuous problem by a finite number of unknowns, on a grid or on a set of panels. Every numerical result carries the error of that replacement, and this collection states the resolution wherever the answer depends on it.

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

Every one of these is a solution, and they lift different amounts. Lift coefficient against the circulation the aerofoil was told to carry, with the Kutta condition's own answer marked. Each point is a complete solve: the sources were found for that circulation and the surface is a wall to within 1.3e-4 at every collocation point. Every member satisfies the equations of motion and the boundary condition, and the lift runs through them at exactly 2Γ/Uc — Kutta–Joukowski appearing as a property of the family rather than as a result about any member of it. Ideal flow round a closed body does not have a unique answer, and the Kutta condition is the extra sentence that picks one.

Nothing but the edge

Cut an aerofoil into panels, put a singularity on each, and require the surface to be a wall. The system that comes out has one more unknown than it has equations, and the row that is missing is not a bookkeeping slip — it is the fact that ideal flow round a closed body has no unique answer at all.

circulation · Panel The wall is what cancels the waves it made. The characteristic net of a minimum-length nozzle designed for Mach 2.4, with 18 waves. The pale lines run from the sharp throat down to the axis, reflect there by symmetry, and run back up to the wall; the wall turns through exactly the angle needed to cancel each one as it arrives, so nothing reflects back into the flow and the exit is uniform at Mach 2.400 and parallel to 0.0 degrees. The area ratio decides the Mach number and this net decides the shape, and the two agree on the exit height to 0.51 per cent at this resolution.

The wall that cancels its own waves

The area ratio of a supersonic nozzle fixes its exit Mach number and says nothing whatever about its shape. What fixes the shape is a wave-by-wave construction in which the wall turns through exactly the angle needed to absorb each expansion as it arrives — and getting it wrong leaves a stream full of oblique shocks at precisely the right Mach number.

compressible · Characteristics A patch of dye, folded. A circle of marker particles carried by the flow, at four times. It is stretched into a filament and folded through itself, and its area does not change at any point of that — which is not visible in the picture, and is the whole difficulty. A scheme that lost eight per cent an orbit would produce a picture indistinguishable from this one.

The area that must not move

A patch of dye in an incompressible two-dimensional flow keeps exactly the area it started with, for ever. Two respectable integrators are put on the same flow: one respects that identically at any step size, the other does not, and the pictures they draw are the same picture.

kinematics · Advection A factor of 1.8 between the extremes. The same numbers as a bar. The spread is exactly 96 over 160/3, which is 1.8 — and it is a spread in a quantity that a single correlation using the hydraulic diameter treats as a constant. A rectangle twenty times as wide as it is deep is 78 per cent more resistant than a circle of the same hydraulic diameter.

A number that is only the shape of the hole

The friction factor times the Reynolds number for laminar flow in a duct is a pure number that depends on the cross-section's shape and on nothing else — not the fluid, not the size, not the flow rate. It is exactly 64 for a circle and exactly 96 for a slot, and the correlation that treats them as the same is wrong by a third.

applied · Internal flow A uniform scalar in a fluid at rest, under two face rules. Nothing is flowing and the scalar starts at one everywhere. The swept-volume rule leaves it at one to the last bit, at every step of the two time units. The midpoint rule moves it by two parts in ten thousand, on a mesh motion that begins and ends in the same place, and the excursion looks exactly like a physical transient.

The mesh that makes its own mass

The transport theorem holds for a region moving at any velocity, which is what makes a moving-mesh calculation possible. Discretised carelessly it is not an identity but an approximation, and a fluid at rest with a uniform density then gains density from the motion of a grid — smoothly, plausibly, and looking exactly like a physical transient.

kinematics · Transport theorem The marginal Taylor number against the axial wavenumber. The smallest Taylor number at which a disturbance of a given axial wavenumber is neutral. Its minimum is 1707.757 at a wavenumber of 3.1158, which are the critical Rayleigh number and critical wavenumber of a layer of fluid heated between two rigid walls — the same numbers, because in the narrow-gap limit the two problems are the same sixth-order eigenvalue problem. This curve is computed by that essay's own solver.

A transition that needs a second number

Fluid between rotating cylinders goes unstable at a Taylor number of 1707.762 — which is the same number, to every digit, as a layer of fluid heated between two rigid walls. It is not an analogy. And two things the number cannot carry decide whether the transition happens at all and what it looks like when it does.

regimes · Taylor The growth rate a discretised sheet has, at every wavelength it can carry. Kelvin–Helmholtz gives a growth rate proportional to the wavenumber and without bound. A sheet represented by N point vortices has pi m (1 − m/N) instead — the same rate at long waves and half of it at the shortest wave the grid carries, with the fastest-growing mode at the grid scale itself. Smoothing the kernel over a length delta moves that mode back to a wavelength the physics chose.

A sheet that cannot stay a sheet

Let a shear layer's thickness go to zero and it becomes a surface across which the velocity jumps. The model is used everywhere in this subject, it is unstable at every wavelength, and the thing it does next is worse: it develops a singularity in its own shape, at a finite time, from a smooth start.

inviscid · Vortex sheet The wall shear, marched to the station where it stops. Howarth's linearly retarded outer flow, marched with an implicit finite-difference scheme from a Blasius profile. The wall shear falls, its slope steepens, and at x = 0.11983 it reaches zero — against Howarth's 0.1198, which is a quarter of a per cent. There is nothing downstream of it: the solution does not continue.

The singularity a layer makes for itself

March Prandtl's equations into an adverse pressure gradient and the wall shear reaches zero with an infinite slope at a finite station, and the solution cannot be continued past it. The singularity is real, it is not a numerical difficulty, and it belongs to the boundary condition rather than to the equations.

viscous · Boundary layer 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 patches' centroids, against the point-vortex circle. Two circular patches of uniform vorticity, advected by nothing but the velocity their own boundaries induce, over one full co-rotation. Their centroids stay within one per cent of a separation of the exact point-vortex orbit — which they must, because the exterior field of a circular patch is the point vortex's and a harmonic function's area average over a disc is its value at the centre.

What a point vortex is not

Two circular patches of vorticity move exactly as two point vortices do — the centroid velocity is the point model's with no correction of any order, and that is a theorem rather than an approximation. The trouble starts the moment they stop being circular, which is immediately.

inviscid · Vortex dynamics 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 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

Named alongside it

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

ConvergenceMeasurementModel limitToleranceCirculationDimensionlessOptimisationRegimeSimilarityConservationEigenvalueEulerian and Lagrangian

All concepts