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

Transport — where it appears

The carrying of a quantity by a flow, whether momentum, heat or a tracer. Its rate is usually written as a diffusivity times a gradient, which assumes the carrying parcels have travelled far enough to have forgotten where they came from.

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

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 The current at the surface is 45° from the wind, and nothing sets that angle. The Ekman spiral drawn as a hodograph: each point is the velocity at one depth, and depth runs along the curve. At the surface the flow is at exactly 45 degrees to the wind that drives it — not approximately, exactly, and independently of the wind, the viscosity and the latitude. By one Ekman depth the flow has turned another radian and lost 1/e of its speed; by three it is a hundredth of the surface value and pointing back the way it came. The angle is a property of the equation having two terms in it, and nothing else.

The layer that stops at a depth

Every other boundary layer grows. This one does not — rotation supplies a frequency, the balance against diffusion supplies a length, and the transport that comes out contains the stress on the surface and not the viscosity underneath it.

viscous · Rotating A material region, and the dye that stays inside it. The same fluid at four times, carried by an unsteady straining flow whose strain rate oscillates. The outline is a circle of the fluid at the first instant, tracked by integrating the velocity field; the shading is a blob of passive dye. The region is stretched to nearly seven to one and its area is unchanged to fifteen decimal places, because the flow is incompressible. The amount of dye inside it is unchanged to thirteen, because the dye is carried by the same fluid.

A rate of change that will not hold still

Three boxes drawn in one flow at one instant give three different answers to how fast the dye inside them is changing — one falling, one falling twice as fast, one rising. All three reconcile with a single material rate, and that rate is zero.

kinematics · Transport theorem 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 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 And the net drift, which is where they disagree. The Stokes drift plus each return flow. All three drift forward at the surface, because the Stokes drift there swamps any return current of the right size. Below that they part company completely: the uniform current has most of the column moving upstream, and the two that satisfy no slip have almost none of it. The reversal depths span 77 per cent of the water column.

The drift a closed box will not allow

A wave in a wave tank carries mass forward, and the tank has nowhere to put it. So a return current appears carrying exactly the opposite transport — exactly, from mass conservation and nothing else. Which fixes a total and leaves the answer anybody wants entirely open.

kinematics · Stokes drift 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 Six orbits that do not close. One parcel's path under a linear deep-water wave of steepness 0.05, at a fifth of a wavelength down, released at the phase that centres the orbit on its release depth. Each loop returns almost to where it began and not quite.

A drift made of two things that average to zero

Stokes drift is usually explained as a parcel spending longer in the forward half of its orbit. That is true and it is not a formula. The formula is a correlation between a displacement and a gradient, each of which averages to exactly nothing, and it splits into two halves that are equal to twelve figures.

kinematics · Stokes drift 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 The layer that grows from a change of surface. The internal boundary layer's height against distance downwind of a change in roughness, with the sublayer inside it that is genuinely in equilibrium with the new surface. The layer grows as the fetch to the four-fifths power and the equilibrium sublayer is a tenth of it.

How far downwind a surface is remembered

Walk from a field into a wood and the wind ten metres above your head is still the field's wind. It takes about a kilometre of trees before a ten-metre measurement is measuring the trees — a hundred times the height it is made at, and a great deal more than most masts are given.

turbulence · Roughness What the fluid at one height is listening to. The weight the fluid two millimetres above a moving wall gives to the wall's velocity a given delay earlier, in water. It peaks at two thirds of a second and has a tail that falls as the delay to the power minus three halves — so the fluid is responding to a broad stretch of the wall's past rather than to a moment of it.

The wall the fluid is listening to

Water two millimetres above a moving wall is responding to what the wall did two thirds of a second ago — most likely. Half of its response is older than four and a half seconds, a tenth is older than two minutes, and the average age of what it is responding to does not exist at all.

viscous · Exact layer The gas that does not stop at the wall. Channel flow profiles with and without slip, at a Knudsen number of a twentieth. The slipping profile does not reach zero at the wall: the gas there is moving, by an amount proportional to the mean free path times the velocity gradient.

Slip is a memory of one mean free path

A molecule arriving at a wall last collided about a mean free path away and carries the velocity from there. Averaged over arrivals and departures, that leaves the gas at the wall moving — by two per cent of the centreline speed at a Knudsen number of a hundredth, and sixty per cent more flow through a microchannel at a tenth.

regimes · Knudsen How long a fluid takes to forget it was not rotating. The fraction of solid-body rotation a container's interior has reached, against time, by the two available routes. The Ekman layers on the end walls pump fluid radially and carry angular momentum inwards in a hundred seconds; diffusion alone would need ten thousand.

How long a fluid takes to forget it was not rotating

Spin a container of water and the fluid inside reaches solid-body rotation in a hundred seconds rather than the three hours diffusion would need. The shortcut is the thin layers on the end walls, and the advantage they give is exactly the reciprocal of the square root of the Ekman number.

viscous · Rotating A float under swell on a rotating planet goes round instead of away. The track of a float at the surface over 1 inertial periods (16.9 hours) after a 8 s swell of amplitude 1 m arrives at latitude 45°, in kilometres, the waves travelling to the right. Without friction the float runs round a circle of radius Uₛ/f = 0.479 km and comes back to where it started every 16.92 hours. With a drag on the Eulerian current it spirals out into a steady drift veered to the right: 24.3 per cent of the drift at 76.0° for a drag of 0.25 f; 70.7 per cent of the drift at 45.0° for a drag of 1 f; 94.9 per cent of the drift at 18.4° for a drag of 3 f. In a non-rotating ocean the same float would have gone 3.0 km straight on.

The drift a rotating planet takes back

In a wave tank the Stokes drift is cancelled by a return current because the tank has walls. The open ocean has none, and the drift is cancelled anyway: the Coriolis force acts on the water's real motion, drives an Eulerian current that answers it, and leaves the depth-integrated transport exactly zero at every viscosity. A float under steady swell with nothing to stop it goes round a circle instead of away.

kinematics · Stokes drift Three profiles that do not depend on the radius. The radial, azimuthal and axial velocities of the flow above a rotating disc, as functions of one similarity variable. The radial one is a jet: fluid thrown outward by the swirl it has picked up, peaking at 0.181 of the local disc speed a fifth of the way through the layer. The azimuthal one falls from the disc's own speed to nothing. And the axial one is the surprise — it does not vanish far from the disc but tends to a constant, so the disc draws fluid down onto itself at 0.8845 times the square root of the viscosity times the rotation rate, at every radius and for ever.

The solution that keeps its nonlinear term

Every exact solution before this one has been exact because the nonlinear term vanished. A rotating disc's does not vanish — at the wall it is the whole of the balance — and the reduction is exact anyway, because the radius divides out of all three momentum equations at once.

viscous · Exact layer Six that are symmetries and five that look like them. Each transformation applied to an exact solution, with the Navier-Stokes residual recomputed from the transformed field by finite differences — nothing differentiated by hand. The six symmetries leave the residual at the differencing floor, a few parts in 10^8. The five near-misses leave between 0.048 and 4.3, which is six to nine orders of magnitude larger. The gap is what makes this a test rather than an illustration: a transformation that is nearly a symmetry does not exist here, and every one of the five is something a reader might reasonably believe.

Why the list is this long

Every textbook list of exact solutions of the Navier–Stokes equations is about a dozen long, and the usual explanation is that the equations are hard. It is not the reason. A similarity reduction is a solution invariant under a subgroup of the equations' own symmetries, so the catalogue of possible reductions is the catalogue of subgroups — and that is a finite, countable object.

viscous · Exact layer The kept transport spirals into nothing as the sea deepens. The net Lagrangian transport as a vector, scaled on the Stokes transport, traced as the water depth increases from a quarter of an Ekman depth to eight, for an 8-second swell with an eddy viscosity of 0.01 m²/s. Shallow water keeps the whole transport pointing with the waves, at the right-hand end. As the sea deepens the vector shortens and swings to the right, crosses the across-wave axis near two Ekman depths, and winds into the origin, which is the open ocean's exact cancellation.

The floor that gives the drift back

In the open ocean the Coriolis force drives a current that cancels a swell's Stokes transport exactly. Over a continental shelf the sea floor holds a stress, and whatever it holds is transport the rotation does not take back. How much survives depends almost only on the depth in Ekman depths; which way it points depends on the wave.

kinematics · Stokes drift The spreading per metre has a best speed; the spreading per second has none. The plate height — how much variance a slug picks up per unit length travelled — against the velocity, for four retention factors. Each curve falls as one over the velocity, because molecular diffusion along the tube matters when the flow is slow, and rises with it, because that is Taylor's dispersion. The minimum is what a chromatograph is run at, and Taylor's effective diffusivity has no such point anywhere.

The best a tube can do is its own radius

Taylor's dispersion coefficient rises without limit with velocity, so there is no best speed to run a tube at. Ask instead how much a slug spreads per metre travelled and an optimum appears at once — and the spreading there is the tube's radius over the square root of three, with no diffusivity and no velocity left in it.

regimes · Dispersion What an oscillation carries, against how fast it is asked. The extra axial diffusivity of a zero-mean oscillation, at fixed velocity amplitude, against the Womersley number. At low frequency it is exactly Taylor's dispersion evaluated at the mean square velocity — the profile has time to be Poiseuille's at every instant. At high frequency the shear is confined to a Stokes layer, the tracer in the core is never sheared, and the transport falls as the cube of the Womersley number.

Transport with nothing transported

Taylor's mechanism needs shear and diffusion and nothing else — and in particular it does not need a mean flow. Oscillate a tube about a fixed position and the tracer still spreads along it, by hundreds of times the molecular rate, which is how a patient is ventilated with a tidal volume smaller than the airway it goes down.

regimes · Dispersion At a coast the swell's water goes back, or goes along. The Lagrangian mean current — the Stokes drift plus the Eulerian current — under an 8-second swell arriving straight at a coast, in seas 8, 25 and 60 metres deep, each as a fraction of its own surface drift against depth as a fraction of the water's. Left, the part across the shore, positive onshore; right, the part along it. The cross-shore parts carry no net water. At 8 metres, inside the Ekman layer, the whole exchange is across the shore; at 25 the return is shared with a current along the coast; at 60 the open shelf's own rotation has already taken the transport back, and the coast has little left to do.

A coast sends the drift back, or sends it along

A shelf sea under a swell keeps part of the Stokes drift's transport, turned to the right by the rotating Earth, and near a coast some of it points at the land. The sea answers with a slope in its surface. In water shallower than an Ekman depth the slope drives the water straight back as an undertow, as it would without rotation; a few Ekman depths down the coast turns the transport into a current along the shore instead, and there the sea can slope down towards the land it faces.

kinematics · Stokes drift The current along the shore runs fastest three Ekman depths down. The depth-averaged Eulerian current along the shore, positive to the left looking shoreward, against the local depth in Ekman depths, for a swell arriving 20° off the shore-normal at the shelf edge and for one arriving head-on. Both make a jet to the right, fastest at 2.9 Ekman depths for the oblique swell (-1.42 mm/s at 37.6 m) and 2.6 for the head-on one (-0.948 mm/s), and a weak current the other way in water shallower than one Ekman depth. Rotation turns the transport the shelf returns, so even a swell with no alongshore component drives one.

A sloping shelf puts the swell's current three Ekman depths down

A swell's Stokes transport arrives at a coast and the sea must send it back or turn it aside. On a flat shelf one depth decided which. On a real shelf, deepening from the beach to its edge, every depth is present at once, and the sea does both: it sets up a few centimetres against the beach, dips a tenth of a millimetre where the depth passes two Ekman depths, and runs a current along the shore that is fastest about three Ekman depths down — wherever the friction puts it.

kinematics · Stokes drift

A thinning gas flows least at one free path across, and then more without end

Pump the gas out of a narrow channel and, for the same pressure gradient, less of it flows — until the channel is about one mean free path high. Thinner than that, more flows, and between infinite plates the flow rate climbs without limit, because the molecules carrying it are the ones skimming along the walls, whose flights only collisions with other molecules can end. Knudsen measured the minimum in 1909. The linearised kinetic equation puts it at a reduced flow rate of 1.538, at a Knudsen number of 0.83, and shows it is not a dip below the two closed forms but a place far above both.

regimes · Knudsen

Named alongside it

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

MeasurementModel limitModel validityBoundary conditionRegimeDiffusionEulerian and LagrangianMixingStokes driftDispersionEkman layerMass conservation

All concepts