Concept

Eddy viscosity — where it appears

A coefficient introduced to make a turbulent momentum flux look like a viscous stress, being a velocity times a length taken from the flow rather than the fluid. It is a modelling assumption, it is orders of magnitude larger than the molecular value, and it is not a property of anything.

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

u⁺ = (1/κ) ln y⁺ + B, integrated rather than assumed. The velocity profile in wall units, produced by integrating the mixing-length closure outward from the wall. The straight portion is the log law and the constants beside it were least-squares fitted to the integrated curve over 50 < y⁺ < 500 — so the 1/κ is a measurement on the drawing rather than the number that was fed in. The viscous sublayer u⁺ = y⁺ comes out rather than being pasted on.

A guess with a constant in it

Prandtl's mixing length is one line — an eddy near a wall can only be as big as its distance from the wall. Integrate it and the whole structure of a turbulent wall profile falls out, sublayer and log region and all. That is a fact about the assumption, and the essay is careful about which.

turbulence · Closure
The jet has to be fed from the sides. The velocity field of the plane jet with streamlines integrated through it. The seven central streamlines run down the jet and spread; the ten started at the top and bottom edges bend inwards and join it, which is entrainment and is a consequence of the solution rather than an addition to it. The dashed lines are the half-speed edges, widening as x^{2/3}. The transverse velocity far from the axis is 5.70e-3 m/s at this station, inward on both sides — a jet is a sink as seen from a distance, which is why two parallel jets pull together.

What a jet keeps, and what it collects

A jet leaving a nozzle into still fluid has no boundary anywhere and one conserved quantity. Its momentum flux is exactly the same at every station downstream; its mass flux is not conserved at all and grows without limit, because a jet is a machine for acquiring fluid it did not start with.

viscous · Free shear
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 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
Four wall models, through the buffer layer. Van Driest's damped mixing length, Reichardt's fit, Spalding's implicit law, and a control with no buffer layer at all — the viscous sublayer joined straight to the logarithm where they cross, at y+ = 11.6. The three fitted models agree with each other to a per cent and a half; the control is nineteen per cent above them at y+ = 10.

Three buffer layers, one friction

Four wall models are put through a pipe. The one with no buffer layer at all is nineteen per cent wrong where the turbulence production peaks and two and a half per cent wrong in the friction; the three respectable ones are within two per cent of each other in the buffer layer and spread by seven in the friction. The answer is not where it was expected.

turbulence · Wall law
One shape, at every station and for every jet. The axial velocity of Schlichting's round jet at four distances, each scaled on its own centreline speed and its own half-width. The four curves are one curve: the shape function (1 + eta²/4)⁻² has no parameter in it at all, and every station of every jet — laminar or modelled, weak or strong — collapses onto it exactly.

Exactly similar, and one number short

The round jet has an exact solution of the Navier–Stokes equations, and a turbulent jet is modelled by the same formula with the viscosity replaced by a fitted constant. The shape is identical. Nothing a measurement of the profile can do will tell the two apart.

viscous · Free shear
The production jumps and the dissipation does not. Production and dissipation against time, through a step change in the strain rate. The production follows the strain immediately — it is the strain squared times an eddy viscosity — and the dissipation moves by less than one per cent at the instant of the step, because it is set by a cascade that has not been told yet.

A dissipation that lags its production

Change the strain rate on a patch of turbulence and the production of energy follows instantly — it is the strain squared. The dissipation moves by less than one per cent, because it is the far end of a cascade that has not been told yet, and the two are out of balance by a factor of six for the next turnover.

turbulence · Dissipation
The stress a closure predicts, and the stress there is. The Reynolds-stress anisotropy through a step change in the strain rate, against what an eddy viscosity gives — which is the equilibrium value at every instant. The real stress takes about a turnover to get there, and during that turnover the closure is wrong by up to sixty per cent.

A closure with no memory at all

An eddy viscosity says the Reynolds stress is the mean strain rate times a number, now. The stress it is standing in for takes a turnover to arrive, so the closure is the zero-frequency limit of a response that has a lag in it — and the curve it is the limit of is the same shape as an aerofoil's lift deficiency.

turbulence · Closure
Two wakes with one drag. The two initial velocity deficits: a slab, which is roughly what a bluff body leaves, and a pair of separated lobes, which is roughly what a body with a jet through the middle of it leaves. Their integrals are identical, so the two bodies have exactly the same drag.

A wake that keeps the drag and forgets the body

Two very different wakes with the same momentum deficit converge to the same profile, because the deficit is conserved and everything else diffuses away. The convergence is a power law rather than an exponential, so it takes two hundred widths for twenty per cent agreement and nine hundred for five.

turbulence · Wake
Rolls turning side by side, with the fastest downwind water where they sink. The fastest-growing mode at a Langmuir number of 0.13, looking downwind, over two roll spacings of 2.89 decay depths and 4 decay depths down. The closed curves are streamlines of the overturning; the dashed curves are contours of the downwind velocity the rolls carry, positive under the lines where the water sinks. At the surface the cross-wind flow converges onto those lines, which is where floating foam and weed collect as windrows. The amplitude is arbitrary, as in any linear mode.

The drift that turns a current into rolls

A current carrying a Stokes drift feels a force the drift makes out of the current's own vorticity, and under a wind that force is unstable. It turns the surface layer into rolls lined up downwind, with windrows where they sink. The rolls need both the current's shear and the drift's; their growth rate sees only the product; and the split between the two decides which motion gets the energy.

kinematics · Stokes drift
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
Where the eddies outconduct the molecules. The ratio of turbulent to molecular heat diffusivity across the pipe at Reτ = 2000, for five fluids, on logarithmic axes. Wherever it is above one the eddies carry more heat than conduction does. For air it passes one inside the buffer layer and reaches 135; for water, earlier and higher. For liquid sodium it never reaches one anywhere: at its peak, halfway to the axis, the eddies carry just over half what conduction carries, and the temperature profile is set by conduction across the whole pipe.

The heat the eddies do not carry

In a laminar layer the temperature and the velocity have different thicknesses in every fluid but one. In a turbulent pipe the eddies carry both, and the difference nearly vanishes — for air, water and oil alike. It does not vanish for a liquid metal, whose molecules conduct heat faster than the eddies can, and the boundary between the two behaviours is a Péclet number of about four hundred at every Reynolds number.

regimes · Peclet
Three pairs released at different distances end on one line. The mean square separation of two parcels against time, both in Kolmogorov units, for pairs released one, thirty and a thousand Kolmogorov lengths apart in turbulence whose integral scale is a million of them. Each starts flat — the pair has not yet moved — and each joins the same line, g ε t³ with g = 0.5, after which nothing about where it started survives. Beyond the integral scale all three become ordinary diffusion, growing as t.

The further apart, the faster they part

Release two specks of smoke close together in turbulence and the distance between them does not diffuse the way a single speck wanders. It grows faster the larger it already is, because larger separations are pulled apart by larger eddies — as the cube of the time, forgetting where it started. The law is ninety years old, its constant is still argued about, and it needs a very large flow to be seen at all.

turbulence · Mixing
The Earth turns a straight river into a very gentle bend. The transverse circulation through the depth of a straight river running at 1.2 m/s at latitude 50°N, in thousandths of the main flow, against height above the bed — and the circulation of the same river without rotation but bent to the left with a radius of 2V/f, 21 kilometres. The two nearly coincide: the Earth's rotation pushes the fast surface water right and the slow bed water left, exactly as a bend does, and for a river it is a bend of about twenty kilometres' radius everywhere at once.

The Earth bends every river a little

In 1860 Karl Ernst von Baer claimed that the Earth's rotation makes rivers in the northern hemisphere cut into their right banks. The Coriolis force on a river does drive a helix that throws surface water towards the right bank, the same helix a bend drives. For a mountain stream it is a ten-thousandth of what the stream's own bends do. For the great slow rivers of the far north it is a fifth to two-fifths, and there Baer's claim is not absurd.

misconceptions · Coriolis scale
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 entrance grows with the Péclet number, peaks, and then shrinks. The distance from a wall-temperature step at which the local Nusselt number has come within 5 per cent of its developed value, in diameters, against the Péclet number, for turbulent pipes at Reτ = 500, 1000, 2000 and 5000, with the slug's laminar line L/D = 0.029·Pe. The liquid metals sit on the line; the entrance peaks at 11.6 diameters at Pe ≈ 4700 for Reτ = 1000 and falls to about three and a half diameters for water.

The longest thermal entrance belongs to neither metal nor gas

A liquid metal carries its heat across a turbulent pipe by conduction, as a laminar flow would, and so its thermal entrance should be long. It is long only when its Péclet number is large. Below a few hundred the entrance grows in proportion to the Péclet number with the constant of a fluid moving as a solid block, and it leaves that scaling close to the threshold at which the eddies first match conduction. The longest entrance in a turbulent pipe belongs to the fluids between the metals and the gases, and how long it is depends on a number the measurements have never pinned down.

regimes · Peclet
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
Let its shear stress lag, and a layer reaches its new equilibrium shape in about seventeen thicknesses. The shape factor downstream of a step from zero gradient to an outer speed falling as the power −0.2 of distance, in layer thicknesses at the step: with the shear stress lagging, by Green's lag-entrainment method; with the entrainment fixed by the layer's shape, as in Head's method; and Clauser's equilibrium layer at the lagged layer's own local β, through Coles's profile with Das's wake strength. Twenty thicknesses on, the lagged layer's shape factor is 1.479 and the shape-only layer's 1.445; eighty on, 1.439 and 1.432.

A layer whose stress lags the gradient separates sooner

A turbulent boundary layer that meets an adverse pressure gradient suddenly cannot answer it at once: the large eddies that carry its shear stress were made upstream, on a flat plate, and they draw in too little fluid for the new gradient for the first twenty or so layer thicknesses. The deficit does not protect the layer. It lets the gradient pile momentum deficit into it, so the layer reaches its new equilibrium shape quickly, its friction falls through the equilibrium value on the way, and in a strong gradient it separates in little more than half the distance a calculation without the lag predicts.

turbulence · Wall law

Named alongside it

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

Model limitClosureCoriolisBoundary layerEkman layerMeasurementMixing lengthModel validityPressure gradientRegimeRotating frameStokes drift

All concepts