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

Strain rate — where it appears

The symmetric part of the velocity gradient, describing how a fluid element is being stretched and squashed. It is what viscosity acts on to produce stress, and it is the half of the velocity gradient that has nothing to do with rotation.

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

A circle of fluid, 0.90 of a gradient time later. A material circle in a uniform velocity gradient, carried by the exact matrix exponential of that gradient. It becomes an ellipse — always an ellipse, for every gradient — and the axes it stretches along are the eigenvectors of the symmetric half. In an incompressible flow the area is unchanged however extreme the distortion, which is the statement that stretching in one direction is squashing in the other.

What a parcel does in the first instant

Drop a circle of dye into a flow and it becomes an ellipse. The velocity gradient that did it splits into a stretch and a spin in exactly one way, the split is not a convention, and one half of it is the reason a fluid has any stress in it at all.

kinematics · Deformation The one place the stretching argument closes. Burgers' vortex: an axisymmetric strain carrying vorticity inwards at exactly the rate viscosity spreads it outwards. The vorticity profile is a Gaussian of radius √(4ν/α), the swirl velocity peaks at 1.12 core radii rather than at the core radius itself, and the circulation reaches its full value by about two. The steady vorticity equation is evaluated on this profile by differencing it, not by re-deriving it.

The spin that feeds itself

Stretch a vortex tube and its spin rises in exact proportion, because the circulation round it cannot change and its area has fallen. Nothing in that argument sets a limit — and the one flow where the limit can be written down exactly puts it at a length of √(4ν/α).

kinematics · Vorticity One flow, two observers, two pictures. The same ideal flow past a circular cylinder, drawn in the frame of the tunnel and in the frame of the undisturbed air. The two are related by subtracting one constant velocity. On the left the flow arrives from infinity, divides at a stagnation point on the nose and closes at another on the tail. On the right the air is at rest far away, the body pushes through it, the streamlines are closed loops, and there is no stagnation point anywhere in the field. Every force, every pressure and every measurement either observer can make is identical.

The picture belongs to whoever is watching

Photograph the flow past a cylinder from the tunnel and it has two stagnation points. Photograph the same flow from a frame moving with the air and it has none at all, and its surface speed is exactly the free stream at every angle. Both pictures are correct and no measurement distinguishes them.

kinematics · Frames Two flows with the same rate of strain, doing different things to a blob. A circle of fluid carried by two flows chosen to have exactly the same rate-of-strain magnitude, drawn at four times. Pure strain pulls it into an ellipse whose axes stay put; simple shear pulls it into an ellipse whose axes rotate as fast as they stretch. Both have zero divergence, so both preserve the area. The difference between them is not the strength of the straining but what the rotation does to the direction being stretched.

Longer, with nothing pulling it

Two flows with exactly the same rate of strain. In one a line of fluid grows by a factor of 148 in five time units; in the other it grows by 10. Turn the straining axes faster than the strain rate and no line grows at all, however hard the fluid is being strained.

kinematics · Material lines Q and the vorticity along a radius of one vortex. Two candidate measures of where the vortex is, along a radius of a Lamb–Oseen vortex. The vorticity is a Gaussian: positive at every radius, so a threshold on it puts the edge wherever the threshold is put. Q — the excess of rotation over strain — changes sign exactly once, at 1.121 core radii, and that radius is a property of the flow rather than of the person drawing it. Inside it, 71.5 per cent of the circulation.

Where a vortex stops

Four criteria decide where a vortex ends, and in two dimensions three of them are the same criterion. The fourth is a knob. And the one that is not a knob is not objective: a co-rotating pair of vortices occupies four per cent of a window to one observer and more than half of it to another.

kinematics · Coherent structures A normal stress the closure makes negative. The Boussinesq closure's first normal stress in a plane strain, against the strain measured in units of the turbulence's own time scale. It crosses zero at S k/eps = 1/(3 C_mu) = 3.704 — eleven per cent above the value the constant was calibrated at — and goes on falling. A variance below zero is not a small error; it is a statement that cannot be true.

The constant that makes a variance negative

Every engineering turbulence calculation in the world rests on one number, C-mu equals 0.09. It is not a property of turbulence. It is the assertion that a particular ratio is ten thirds, which is true in one flow — and eleven per cent above that flow the same closure reports a mean square below zero.

turbulence · Closure One pair of strainings, two orders, two lengths. The stretch of the most-stretched material direction against time, for a simple shear followed by a pure strain and for the same two in the other order. The two curves are identical until the swap and separate afterwards, ending a factor of 2.16 apart.

Two strainings, and the order they came in

A material line is stretched by a shear and then by a pure strain, and then by the same two in the other order. Every instantaneous measure of how hard the fluid was being worked is identical in the two cases. The lengths at the end differ by a factor of 2.16.

kinematics · Material lines Kirchhoff's rotation rate, which a point vortex does not have. A patch of uniform vorticity bounded by an ellipse turns rigidly at omega a b/(a+b)², a rate that depends on the shape alone. It is largest for a circle, where it is unobservable, and falls away as the patch is drawn out. A point vortex has no shape and therefore no entry on this axis at all.

The shape a vortex keeps

Outside a circular patch of uniform vorticity the flow is exactly the point vortex's — not nearly, exactly — so replacing one by the other looks free. It is not. The patch has a shape, the shape has a rotation rate of its own, and there is a strain above which no shape exists at all.

inviscid · Vortex patch The two dissipations, side by side. The strain form on the left and the enstrophy form on the right, for the same field, on the same scale. They have their maxima in different places — 0.46 apart on a box of side 2 pi — and neither is a smoothed version of the other. One says the dissipation is in the strained regions and the other says it is in the rotating ones, which is nearly a complete disagreement about what a turbulent flow is doing.

Equal on average, and nothing else

The rate at which a fluid turns motion into heat can be written two ways, and every textbook says the two are equivalent. Their averages are equal to fourteen decimal places. Point by point they are uncorrelated, and their maxima are in different places.

turbulence · Dissipation The scalar spectrum, with its two ranges. A model scalar spectrum at a Schmidt number of two thousand — dye in water. Below the Kolmogorov wavenumber it is Obukhov and Corrsin's five-thirds, inherited from the velocity; above it there is no turbulence left and the spectrum is Batchelor's minus one, which contains no velocity spectrum at all.

The scalar has its own cascade

Below the Kolmogorov scale there is no turbulence left, and a dye stirred into the flow goes on cascading anyway — on a spectrum whose exponent is minus one and whose amplitude contains no velocity spectrum at all. Resolving it costs the three-halves power of the Schmidt number, which for dye in water is a factor of ninety thousand.

turbulence · Mixing 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 Four kinds of critical point, and the curve that separates them. The invariants of a trace-free velocity gradient, with the discriminant curve 27R²/4 + Q³ = 0 drawn through them. Inside the two upper lobes the cubic has one real root and a complex pair, which is a spiral being stretched along its own axis on the left and squeezed on the right; below the curve all three roots are real and the point is a node with two saddle directions. Of 820 random incompressible gradients, 509 land in the spiral region and 311 in the real one. A plane flow is the vertical line R = 0 and nothing else, which is why a plane has two kinds and space has four. The marked points are the cases the calculation checks that fall inside this window; the two vortex cases it also checks sit at Q = 3.25 and |R| = 4.25, off the top corners, because a window wide enough to hold them would flatten the curve the figure is about.

Two kinds is a plane flow's privilege

A plane incompressible flow has a saddle or a centre and nothing else, and the proof is one line about a trace. The same line in three dimensions constrains three numbers instead of two, which is far less, and what it leaves is four kinds of point separated by a curve — with the one a plane cannot have being the structure the whole of turbulence is made from.

kinematics · Topology Four fluids under the same stress, and the four profiles that result. Each fluid carries the same linear stress and answers it differently: a Newtonian fluid with a parabola, a shear-thinning one with a blunter profile, a shear-thickening one with a sharper one, and a Bingham plastic with a plug in the middle where the stress is below its yield point.

The stress a pipe knows

A capillary viscometer measures a pressure drop and a flow rate and reports a viscosity. The first half of that inference is a force balance and is exact for any fluid there is; the second needs the slope of a whole flow curve, which is the experiment the instrument was bought to avoid.

viscous · Non-newtonian The one curve of the four that measures a gradient. A row of particles released along a line at one instant, seen at 4 later times, in a layer profile. Each line has been displaced by the local velocity times the elapsed time and by nothing else, so dividing the displacement back out returns the profile exactly — checked here to machine precision at every one of two hundred heights. A streamline, a pathline and a streakline each report where fluid went; this one reports how fast neighbouring fluid was going relative to its neighbour, which is the quantity a boundary layer is made of and the one the other three never show.

The curve that measures a gradient

Three of the four curves drawn through a flow answer the same question — where did the fluid go. The fourth answers a different one. A line of particles released together is displaced by the local velocity and by nothing else, so its shape is the velocity profile, and dividing the elapsed time back out returns that profile exactly rather than approximately.

kinematics · Flow curves Four solved flows, and the minimum is on the surface in every one. Sampling the whole exterior of each body on a grid and comparing the lowest pressure found there with the lowest found on the surface. The surface wins by a margin that is not marginal — between 0.18 and 0.56 in pressure coefficient — and it wins for a reason rather than by luck: the pressure of an irrotational flow is superharmonic, and a superharmonic function has its minimum on a boundary.

The lowest pressure is on the body

In an ideal flow the minimum pressure is always on a surface — not usually, not for the shapes people draw, always. The proof is an identity about the velocity gradient, and the identity says exactly which flows are exempt.

inviscid · Ideal flow A pair of points meets on R = 0, the one line an index can change on. The paths of the ABC flow's stagnation points across the (R, Q) diagram as C rises from 0.3 towards √2 with A = B = 1. The four points of index +1 share one path on the left and the four of index −1 its mirror image on the right. At C = 0.3 they sit at R = −0.088, Q = −1.045; at C = 1 they touch the discriminant curve at R = −0.707, Q = −1.500, where the strain has a repeated rate, and turn away from it without crossing; and as C approaches √2 they run in to R = -7.5e-3, Q = −2.000. Crossing into a lobe would have changed a node into a focus, which a Beltrami flow's stagnation point cannot be; reaching R = 0 is where each meets a partner of the other index.

The sign a stagnation point carries in space

In three dimensions a stagnation point's index is the sign of one determinant, and that determinant is minus the R of the invariant diagram — so the diagram's left and right halves are the two indices. An exact Euler flow in a periodic box has eight such points, four of each sign, never a spiral among them, and they can only disappear in pairs that meet on the one line where the sign is allowed to change.

kinematics · Topology Every gradient in one plane, and where each one goes. The two invariants of a trace-free velocity gradient, R across and Q up. Under the restricted Euler equation the combination 27R²/4 + Q³ never changes, so every trajectory is one of these curves, and R never decreases, so every curve is travelled from left to right. Vieillefosse's line, where the combination is zero, has two branches: the left one runs into the origin and is the only way to reach it, and the right one is where every other trajectory ends, at infinity, in finite time.

A gradient left to itself

Follow a parcel's velocity gradient with nothing acting on it but its own square and the part of the pressure a single parcel can know about. Every starting gradient but a set of measure zero reaches infinity in a finite time, and it gets there as a sheet with its spin lying in the sheet.

kinematics · Deformation Which bends of a vortex pair grow, and how fast. The square of the growth rate of a sinusoidal bend of a trailing vortex pair, in units of one e-fold per descent time, against the bend's wavenumber times the spacing, for a core of elliptic loading's size. Where the square is positive the bend grows. The symmetric mode, in which the two vortices bend as mirror images, grows over a long-wave band and fastest at a wavelength of 8.54 spacings, by 0.827 e-folds in the time the pair takes to sink one spacing. The antisymmetric mode, bending in step, is stable across the same band.

A wake ends by bending, not by fading

The two vortices behind an aeroplane do not simply weaken until they are gone. In still air each one bends in the other's strain, the bend grows by a factor of e every thirty-three seconds behind an airliner, and after a couple of minutes the pair has pinched itself into a chain of rings. The wavelength it chooses is eight and a half times the spacing, and the vortex cores hardly enter.

circulation · Tip vortex Every particle leaves a flow that is a vortex at every instant. Eight fluid particles starting on a small circle in Haller's rotating-saddle flow, followed for 2.4 time units. At every instant the velocity gradient has complex eigenvalues and a positive Q — the flow is a vortex by every criterion that reads a snapshot — and every particle spirals outward, its distance growing as e to the time. The fluid is not held; it is flung.

Every snapshot says vortex, and every particle leaves

There is a flow in which the velocity gradient has complex eigenvalues at every point and every instant — a vortex by every criterion that reads a snapshot, in the frame the flow is measured in — and in which every fluid particle is flung away exponentially. The snapshot is not wrong about the gradient. It is wrong about the fluid, and it is wrong exactly when the strain's axes turn.

kinematics · Frames Two shapes for each strain, and one of them holds. The strain rate, over the vorticity, at which an elliptical patch of aspect ratio λ stands still — Moore and Saffman's relation — rising to its maximum of 0.1501 at λ = 2.89 and falling again. Below the maximum there are two steady shapes for each strain: a rounder one, on which every disturbance computed here stays bounded, and an elongated one, which comes apart. Above it there is no steady shape at all.

A strained vortex holds until it has no shape to hold

A patch of vorticity in a strain has two steady shapes for every strain below 0.150 of its vorticity, a rounder one and a longer one, and none above. The longer one comes apart at the slightest nudge. The rounder one, computed with disturbances of two, three, four and five lobes, never does: it nods and holds right up to the strain at which it ceases to exist. So the existence limit is the real limit, and past it a vortex is not shattered but stretched — lingering first near the shape it has lost, for a time that grows as the fourth root of how close the strain is to the limit.

inviscid · Vortex patch Remembered pressure turns the runaway into a teardrop. The velocity gradient's two invariants, Q against R, both scaled by the mean enstrophy, for an ensemble of parcels whose pressure remembers one Kolmogorov time of deformation, at a memory of a tenth of a large-eddy time. The restricted Euler equation sends every parcel off to infinity along the right-hand branch of Vieillefosse's curve, dashed. With the remembered pressure the parcels stay: they crowd along that branch in the strain-dominated quadrant and above the axis on the left where vortices are being stretched — the teardrop measured in turbulence.

The pressure a parcel remembers keeps it finite

The restricted Euler equation follows a parcel's velocity gradient with the pressure's shape thrown away, and every gradient it follows blows up. Give the parcel back a pressure that remembers how its neighbourhood was deformed over the last Kolmogorov time — nothing more — and no gradient blows up at all. The ensemble settles into the teardrop measured in turbulence, its vorticity lines up with the middle strain axis, and its intermittency grows as the memory shortens. What the memory cannot do is keep the one identity homogeneity demands, and that miss says where the rest of the pressure lives.

kinematics · Deformation Every parcel goes straight while the streamlines wave. A uniform stream across the page carrying a pattern of transverse velocity with it, v = 0.6 sin(x − t): the streamlines at one instant, which wave, and the paths of five parcels over the next six time units, which are straight lines. Each parcel keeps the transverse velocity it started with, because the pattern moves with it; the streamlines are a snapshot of a pattern that is sliding past, and no parcel ever follows one.

A parcel goes straight while the streamlines curve

In a steady flow every parcel accelerates while the picture never changes. The opposite also happens: a flow whose picture changes all the time while no parcel accelerates at all, because the local and convective halves of the acceleration cancel exactly. Such a flow has no pressure gradient anywhere, and that is a severe demand. An incompressible flow can meet it only as a pure shear at every point; a cloud that can compress meets it freely and pays later, when its straight paths cross in a caustic.

kinematics · Acceleration Moving the saddles never lowers the quarter; fast large wobbles raise it. The smallest Stokes number at which a cloud of particles in the oscillating lattice folds within sixty time units, against the oscillation's amplitude ε, at frequencies of 1, 3 and 10 times the saddles' strain rate. With ε = 0 it is the steady lattice's 0.251. Moving the lattice never takes it below a quarter; at ε = 1 it is 0.283, 0.311 and 0.324 at the three frequencies.

Wobbling saddles keep the quarter, and raise it

Heavy particles in a steady lattice of vortices fold onto the cell walls only above a Stokes number of a quarter, the converging saddles' own threshold. Random flows fold at every Stokes number, so it was natural to suspect the quarter of belonging to steady saddles. Set the lattice wobbling and the particles cross from cell to cell, seventy per cent of them in forty time units, and still none folds below a quarter. The threshold belongs to the strongest strain the flow has; moving the saddles only moves particles out of reach of it, and raises the threshold.

kinematics · Acceleration A stronger strain folds particles below a quarter only if it lasts. The Stokes number above which particles in a lattice with pulsing strain 1 + a sin ωt fold onto its walls, against the pulse's frequency, for amplitudes of 0.15, 0.3 and 0.5, with the steady lattice's quarter dashed. A slow pulse folds them at the peak strain's quarter, 1/(4(1 + a)): 0.194 at a = 0.3 and ω = 0.05, against 0.192. A fast one folds them only at the mean's: 0.249 at ω = 5.

A stronger strain folds a cloud only if it lasts

Heavy particles in a lattice of vortices fold onto the walls above a Stokes number of a quarter of the saddles' strain. Let the strain pulse and the threshold moves — down to the quarter of the peak strain when the pulse is slow, back up to the quarter of the mean when it is fast — because a particle responds to the strain it has been in for about four strain times, not to the strain at the instant. In a random strain the same rule makes the folds arrive in bursts, at the excursions that are strong and long at once.

kinematics · Acceleration

A cloud folds where a saddle is strong, not when

When every saddle of a lattice strengthens and weakens together, a cloud of heavy particles folds in one step, at the strain's first strong excursion. Give each saddle its own random strength and the folding spreads out into a steady rate — and falls, because a particle now has to be at a saddle while that saddle is strong. The rate follows the Gaussian tail of the strain, at about a tenth of the time the strain spends past its fold condition, and a longer-lived strain helps only up to a few strain times.

kinematics · Acceleration

Named alongside it

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

Model limitVorticityVelocity gradientMeasurementQ-criterionSaddleCirculationEigenvalueObjectivityParticleStagnation pointThreshold

All concepts