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

Mixing — where it appears

The process by which two parts of a fluid become indistinguishable, which requires molecular diffusion at the last step and stirring to make that step short. Stirring alone is reversible and mixes nothing, which is why a well-stirred flow can be unmixed exactly.

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

The same patch, 6 periods later, in two flows. A round patch of 848 marked particles, advanced 6 periods by the blinking flow and by a steady flow of the same strength. The steady flow has drawn the patch into a smooth ribbon along a streamline and every particle in it is still on the streamline it started on; the blinking flow has folded the patch through itself repeatedly and its particles are spread across the whole region. Neither flow has any diffusion in it and neither has lost a particle. The difference between them is that one depends on time.

No randomness, and it mixes anyway

A steady two-dimensional flow cannot mix, however fast it is stirred, because its trajectories are its streamlines. Switch two vortices on and off alternately and the same fluid, obeying an exact map with nothing random in it, folds a patch of dye through itself until neighbouring particles separate by a factor of a thousand in six periods.

kinematics · Advection A million times faster, and the constant is 48.0. The effective diffusivity along a pipe, divided by the molecular one, against the Péclet number — both logarithmic. Below Pe ≈ 7 the tracer simply diffuses and the curve is flat at one. Above it the dispersion is all Taylor's, rising as the square of the Péclet number, so a thousandfold Péclet number is a millionfold enhancement. The constant in D(1 + Pe²/48) is not quoted here: it is recovered from a numerical solution of the cell problem across the section, giving 48.0000 for a tube and 52.5 for a plane channel, which is Aris' 2/105.

Two slow things make a fast one

Shear stretches a slug of dye and mixes nothing, because it is reversible. Molecular diffusion is hopeless at any scale bigger than a hair. Put the two together in a pipe and the dye spreads along it with an effective diffusivity two million times the molecular one — which gets larger as the molecular one gets smaller.

regimes · Dispersion Four turns of a wave, and the parcel is not back. Two parcels traced through four periods of a linear deep-water wave of steepness 0.1, by integrating the exact velocity field. Each orbit is very nearly a closed circle and misses closing by a little, every time, in the same direction — that miss is the whole of the Stokes drift. The near-surface parcel advances 0.0357 of a wavelength over the four cycles and the one a tenth of a wavelength down advances 0.0122, a third as far — because the drift falls off twice as fast with depth as the orbit's own size does.

The drift in a wave that has none

The velocity at any fixed point under a passing wave averages to exactly zero, and every parcel of water in it moves steadily forward anyway. The orbits do not close, they miss by the same amount every time, and the missing amount is the square of the steepness times the wave speed.

kinematics · Stokes drift A cross at 30.0° for every wavelength it makes. A body oscillating at ω = 0.5N in a stratified fluid, and the four beams along which its energy leaves. The angle is arccos(ω/N) from the vertical — 30.0 degrees from the horizontal here — and it is the same for every wavelength the body excites, because the dispersion relation has no length in it. The short strokes are the crests, which lie along the beams rather than across them: the phase advances perpendicular to the energy, and the two are exactly at right angles. Raise the frequency and the cross closes towards the vertical; reach ω = N and it shuts entirely, because nothing above the buoyancy frequency propagates.

The number that stops the mixing

A fluid whose density falls with height resists being stirred, and the resistance has a threshold at exactly one quarter, from an energy balance with no fluid mechanics in it. The waves such a fluid carries are stranger still — their frequency decides the direction they travel in and says nothing about their wavelength.

turbulence · Stratification A millionth, doubling every three-quarters of a second. The separation of two trajectories started a millionth apart, on log axes against time. It grows as a straight line until it saturates at the size of the attractor, and the slope of that line is the largest Lyapunov exponent — measured here on the system by renormalising a nearby pair, not quoted. Sensitive dependence is what the straightness of the line means.

The randomness that is not in the equations

Turbulence is described in the language of statistics — means, variances, spectra, probability distributions — and none of that language appears in the equations it is a description of. The Navier–Stokes equations have no random term anywhere in them. What is random is the observer's ignorance of the initial data, and the flow's habit of amplifying it.

misconceptions · Randomness An hour for every tenfold, for ever. How long a forecast lasts, against how well the initial state is known. The relation is T = ln(tolerance/error)/lambda — exactly logarithmic — so improving the measurement by a factor of ten buys exactly the same extra time every time: ln(10)/lambda, which for this flow is 24.5 time units. It does not get harder and it does not get easier.

An hour for every tenfold

Turbulence is deterministic and unpredictable, and the exchange rate between those two is exact: measuring the initial state ten times better buys the same extra forecast time every time, for ever. A constant, and it belongs to the flow rather than to the instrument.

misconceptions · Randomness An oscillation with no mean, and the steady flow it drives. The steady second-order velocity through a Stokes layer, in units of U U′/ω. The first-order flow averages to zero at every height; the average of its own nonlinear term does not, and the pale curve is that forcing. Integrating it twice across the layer, with no slip at the wall and no stress at the top, gives a steady velocity that rises through the layer and settles at -0.749998 — Rayleigh's −3/4, which was not put in anywhere. Beyond about five layer thicknesses nothing more happens, which is why the number is a boundary condition for the flow outside.

An oscillation with somewhere to go

Shake a fluid back and forth over a body and it develops a steady circulation that never reverses. The driving flow has no mean at all; the mean of its own nonlinear term does, and integrating that twice across the oscillatory layer gives a slip velocity of exactly three-quarters of U dU/dx over the frequency.

viscous · Streaming 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 Ballistic first, diffusive later. The mean square displacement of a parcel carried by a fluctuating velocity, from Taylor's 1921 integral with an exponential correlation of time scale 1. Below the correlation time it follows the straight-line law u²t² exactly — the parcel has not yet changed its mind — and above it the curve joins the diffusion law 2u²T·t, offset by the head start the ballistic phase gave it. The eddy diffusivity is u²T = 1.000, and it is a property of the correlation rather than of any equation of motion.

How far a parcel gets

Turbulent transport is one integral. A parcel carried by a fluctuating velocity goes in a straight line while it still remembers its own motion and performs a random walk once it has forgotten, and the crossover is the correlation time — so an eddy diffusivity is not a property a flow has, it is the limit of a measurement that has run for long enough.

turbulence · Mixing 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 Five numbers, one flow, one name. A hyperbolic-tangent shear layer with a matching density profile, measured five ways, each of which appears in the literature as "the Richardson number". The minimum gradient value is what Miles' and Howard's theorem is about and is the only one the quarter belongs to. The bulk numbers depend on which thickness and which velocity difference were used; the depth-averaged one depends on how far from the layer the measurement extended, and grows without limit as it extends further.

Five numbers, one name

The Richardson number has a threshold at a quarter, and the quarter belongs to one of the five quantities that go by the name. On a single tanh shear layer they run from J to seventeen J, and the largest of them grows without limit as the measurement is extended further from the layer.

turbulence · Stratification One streamline, sectioned, in two steady flows. Every time a single streamline crosses the plane z ≡ 0 going upwards, a point is plotted. On the left the flow is integrable and the points lie on a curve, however long the trajectory is run. On the right one coefficient of the same exact solution has been changed and the same single streamline scatters over a sixth of the plane. Both flows are steady, both are incompressible to machine precision, and both are exact solutions of the Euler equations.

Steady, three-dimensional, and mixing anyway

A steady flow that solves the Euler equations exactly, with its vorticity equal to its velocity to six parts in ten thousand million — and one of its streamlines wanders through a sixth of the box while another, started nearby, lies on a curve for ever.

kinematics · Advection One rate per moment, and none of them the same. lambda_p = ln<l^p>/(p t) against p. As p goes to zero it is the Lyapunov exponent, the rate of the typical element; at p = 1 it is the rate of the average length, which is nearly twice as large. If ln l were exactly Gaussian this would be a straight line with the Lyapunov exponent as its intercept, and the departure from that line is the same multifractality the velocity increments have.

The stretching rate that is not one number

A material line in a flow gets longer, and there is a theorem saying its length grows at a definite exponential rate. There is also a rate at which the average length grows, and it is nearly twice as large — and a different rate for every moment of the distribution.

kinematics · Material lines 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 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 Five points, and where each one's fluid came from. Back-trajectories through six units of time in an unsteady double gyre. Each curve ends at the place the fluid now at the marked point started; nothing about that place can be read off the velocity at the marker.

A scalar is a record of where its fluid was

A conserved scalar has no value of its own. Its value at a point is whatever it was at the place that point's fluid started from, which makes a dye field a photograph of the past — and makes the map from now to then the only thing in the flow that carries the past at all.

kinematics · Advection Instability up to a quarter, and none past it. The fastest growth rate of a stratified shear layer against its Richardson number, on a profile whose gradient Richardson number is the same at every height. It falls smoothly towards zero and reaches it at a quarter: at Ri = 0.2499 the fastest mode still grows at 0.00106, and at 0.26 the solver finds no unstable mode at all.

Sufficient, and not necessary

A stratified shear layer whose Richardson number exceeds a quarter everywhere cannot go unstable. That is a theorem with an exact number in it. What it does not say — and what it is constantly read as saying — is that a layer below a quarter will.

turbulence · Stratification The stretching a window of history did, drawn as a field. The largest finite-time Lyapunov exponent over eight units of time in the double gyre, darkest where two neighbouring parcels were pulled furthest apart. The bright crest is a curve across the domain, and it is a property of the eight units rather than of any instant inside them.

A boundary that only exists over a window

The curve that separates fluid going one way from fluid going another is not in any snapshot of the flow. It is the crest of a field built from a stretch of history, it moves when the stretch is changed, and reversing the direction of time gives a different curve entirely — both of them real.

kinematics · Coherent structures 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 A duct that does not change, and a parcel that does. A contraction of area ratio four, with the parcel marked at five equal intervals of time. The walls do not move, the field at every point is the same at every instant, and the spacing of the markers grows because the parcel is carrying its own history through the duct.

How long the fluid has been in there

Age is the simplest thing a flow can remember. It obeys the shortest transport equation in the subject — its material derivative is one — and no instrument pointed at a steady flow can read it, because a steady flow's every field is constant and its fluid is getting older all the time.

kinematics · Acceleration There and back again. A blob of a hundred and twenty tracer particles at the start, after four time units of stirring, and after the same four run backwards. The third set is drawn over the first and the worst particle is 1.4·10⁻⁹ from where it began.

Reversible, and unusable

Ideal flow has no arrow of time in it. Run a stirring backwards and the dye comes back — here to 1.4 parts in a thousand million. Nudge the state by a hundred-millionth first and the same reversal returns a blob almost five hundred times further from home than the nudge was large.

inviscid · Vortex dynamics

Named alongside it

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

MeasurementModel validityTransportThe Lyapunov exponentModel limitFlow mapMemory kernelPredictabilityRegimeTurbulenceAveragingVorticity

All concepts