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 16 essays across 4 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 two per cent of a window to one observer and twenty-five 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

Named alongside it

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

VorticityMeasurementModel limitVelocity gradientCirculationStagnation pointDissipationMixingObjectivityQ-criterionContour dynamicsCritical point

All concepts