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

Wall shear — where it appears

The tangential stress a fluid exerts on a surface, equal to the viscosity times the velocity gradient at it. Its vanishing is the definition of separation, its integral is the friction drag, and its square root over the density is the velocity scale a turbulent wall layer is written in.

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

A wave that dies within one wavelength — 100 Hz in air. The velocity profile above an oscillating wall, at eight phases of one cycle, with depth in units of δ = √(2ν/ω). The motion is a wave travelling into the fluid, and its amplitude falls by 1/e in the same distance it turns by one radian — so it is dead within about one wavelength, and the fluid three δ up hardly knows the wall is moving at all. This is one of the very few exact solutions the Navier–Stokes equations have. The dashed line is one fixed distance above the wall in millimetres: in these units it climbs as the square root of the frequency, which is the whole of how far the motion reaches.

The wall that shakes

Slide a wall back and forth in its own plane and the fluid above it does not follow — a wave travels upwards into the fluid and dies within one wavelength. The depth it reaches is √(2ν/ω), it contains no length from the geometry at all, and the whole thing is one of the very few exact solutions the Navier–Stokes equations have.

viscous · Exact layer Two bills, and the radius that settles them. The cost of a vessel against its radius: the pumping power, which falls as the inverse fourth power, and the price of owning the fluid and the wall, which rises as the square. Their sum has a minimum, found here by golden-section search and agreeing with the closed form to eight figures. At that radius the pumping bill is exactly a third of the total — for every set of constants, because it follows from the two exponents alone.

The radius that costs least

A vessel that carries a flow costs two things to own — the power to push fluid along it and the price of the tissue itself. Minimising the sum gives a best radius, the best radius makes flow proportional to radius cubed, and the rule that follows is a statement about a photograph that came out of a cost function.

applied · Branching α = 13.0: a plug in the middle and everything happening at the wall. The velocity profile at eight phases of one cycle, at a Womersley number of 13.0 — the human aorta's at rest. The core moves almost as a solid plug, because viscosity cannot reach it within a cycle; all the shear is in a layer of thickness √(ν/ω) = 0.77 mm against a radius of 10.0 mm. At some phases the fluid near the wall is moving backwards while the core still moves forward, which is the reversal a Poiseuille profile can never show and which is routinely measured in arteries.

Too fast for a profile

A pipe carrying a steady flow has a parabolic profile. Make the pressure oscillate and one number decides whether it still does — and above about ten the core moves as a plug, a quarter of a cycle behind the pressure, with the fastest fluid in a ring near the wall rather than on the axis.

regimes · Womersley The wall makes vorticity at a rate with no viscosity in it. At a stationary wall the momentum equation collapses to ν ∂²u/∂y² = (1/ρ) ∂p/∂x, and the left-hand side is the diffusive flux of vorticity out of the surface. So the pressure gradient along the wall is the vorticity source, and the viscosity that made the no-slip condition necessary has cancelled out of what the condition produces. The curve is that flux across the Falkner–Skan family, computed from profiles solved by shooting and differenced at the wall; the straight line is the pressure gradient each of those flows has. They agree to 2.3e-14. At zero pressure gradient the flux is exactly zero: a flat plate creates no vorticity at all after its leading edge, and everything in its layer arrived from there.

Where vorticity comes from

Every scrap of vorticity in a flow past a body entered through its surface, and the rate at which it enters contains no viscosity at all — it is the pressure gradient along the wall. A flat plate makes none, and a closed body makes exactly as much of each sign.

kinematics · Vorticity The air a wing carries along, and how little of it there is. The Blasius profile, in the wing's frame, with the free stream at one. No slip says the air at the surface is at rest relative to the surface, which in the ground's frame means it is moving with the wing — but only exactly at the wall. The deficit, integrated across the layer, is the displacement thickness: at a Reynolds number of 1e+6 and a metre of chord it is 1.72 millimetres of air moving at flight speed, which is the whole of what is 'carried'. The step drawn on the axis is that same deficit as a solid slab. A wing does not drag a blanket of air with it; it leaves a boundary layer behind it, and the layer is made of air that keeps being replaced.

The air a wing does not carry

No slip says the air touching a surface moves with it, and the usual reading is that a wing drags a blanket of air along. The blanket is 1.7 millimetres thick per metre of chord, it is different air every instant, and the drag it costs falls as it gets thicker.

misconceptions · The no-slip condition 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 The plateau that is the log law. y⁺ du⁺/dy⁺ across a channel, at five Reynolds numbers. Millikan's argument says this quantity must be constant wherever neither the viscous length nor the channel width may appear, and its value there is 1/κ. At Re_τ = 180 there is no flat part at all; at Re_τ = 100,000 it is flat over 2.16 decades and gives κ = 0.4120. The log law is a statement about a limit, and this is the picture of the flow approaching it.

The layer with no length in it

The logarithm in a turbulent wall profile does not come from any model of turbulence. It comes from a region where neither of the flow's two lengths is allowed to appear, and where a velocity gradient therefore has nothing to depend on but the distance to the wall. The constant in it has never been derived from anything.

turbulence · Wall law The wall's condition, on its way to the middle. Five profiles across the half-channel, from just inside the entrance to fully developed, each drawn at the station where it occurs. The march starts from a slab of uniform flow and never assumes a shape: what arrives at the far end is a parabola, with a centre-line speed of 1.4979 times the mean against the exact 3/2 and a momentum flux of 1.1995 against 6/5. Notice what the middle does while the edges are being slowed: it speeds up, because the flow rate is held, and that acceleration is what the entrance's extra pressure drop pays for.

How far before a duct forgets what was fed into it

A pipe is always drawn with its answer already in place. Getting there takes a distance proportional to the Reynolds number, which means a more viscous fluid is done sooner — and the entrance costs a fixed number of dynamic pressures however long the pipe is.

viscous · Entrance Where the first grid point may go. The error in the friction a wall treatment infers, against the height of the first grid point, at Re_τ = 20,000. The velocity fed to each treatment is the closure's own, so what is plotted is the modelling of the boundary condition with nothing else in it. The log-law function is exact between y⁺ 30 and 10261 and is 64 per cent wrong at y⁺ = 1; the sublayer treatment is exact below y⁺ = 5 and hopeless above it; and the blend that most codes ship is within a few per cent everywhere and exact nowhere.

What a code says to a wall

A calculation that cannot afford to resolve the viscous sublayer has to tell the wall something else instead, and what it tells it is the law of the wall — an asymptotic result, applied at one grid point, on the assumption that the point lies in a region the calculation has not checked exists. Where it does, the answer is exact. Where it does not, the friction is out by tens of per cent, and refining the grid makes it worse.

turbulence · Wall law The wall shear, marched to the station where it stops. Howarth's linearly retarded outer flow, marched with an implicit finite-difference scheme from a Blasius profile. The wall shear falls, its slope steepens, and at x = 0.11983 it reaches zero — against Howarth's 0.1198, which is a quarter of a per cent. There is nothing downstream of it: the solution does not continue.

The singularity a layer makes for itself

March Prandtl's equations into an adverse pressure gradient and the wall shear reaches zero with an infinite slope at a finite station, and the solution cannot be continued past it. The singularity is real, it is not a numerical difficulty, and it belongs to the boundary condition rather than to the equations.

viscous · Boundary layer 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 A force that depends on the molecular scale only through its logarithm. The force per unit length needed to move a 30° contact line of water at 1 mm/s, out to 1 mm, against the slip length on a logarithmic axis, from a picometre to a tenth of a millimetre. Each tenfold change in the slip length moves the force by the same fixed amount, so eight decades of the most uncertain length in the problem change the answer by a factor of about twenty — and at zero slip length the line keeps rising without end. The dashed curve takes the exact wedge's angle factor with a sharp cutoff; the solid one the thin-film wedge with Navier slip.

The drop a no-slip wall would never let spread

Liquid touching a solid moves with it, and nearly everywhere that is as close to exact as anything in fluid mechanics. At the edge of a spreading drop it cannot be: the stress in the corner rises as one over the distance from the edge, and the force needed to move the edge is infinite. Something slips over a nanometre, and because the answer depends on that length only through its logarithm, a drop spreads at almost the same rate whatever the something is.

misconceptions · The no-slip condition

Named alongside it

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

Boundary layerBoundary conditionViscosityDimensionlessModel limitThe no-slip conditionSeparationAdverse pressure gradientClosureConservationDisplacement thicknessExact solution

All concepts