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

Boundary layer — where it appears

The thin region next to a surface in which viscosity is not negligible. Its thickness sets the friction, its survival decides whether the flow separates, and outside it the flow behaves as though the fluid had no viscosity at all.

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

Flow past a cylinder at Re 40. A real fluid past a circular cylinder. At low Reynolds number the flow closes up behind the body much as the ideal theory says; as it rises the flow separates and a region of reversed flow appears behind, which is where drag comes from.

Everything happens in a layer you cannot see

Air has so little viscosity that ignoring it works almost everywhere. Almost everywhere leaves out a film next to the surface, perhaps a millimetre thick, and that film decides drag, stall and whether an aircraft flies at all.

viscous · Boundary layer Two drag laws, one derived and one fitted. Flat-plate drag coefficient against Reynolds number, laminar and turbulent, on log axes. The laminar curve is Blasius' similarity solution, solved by shooting; the turbulent one is the 1/7-power correlation, which was fitted to experiment and is drawn dashed-in-kind to keep the difference visible. Their slopes differ — −1/2 against −1/5 — so the gap between them widens rather than staying put.

The number that is not a number

Transition Reynolds numbers are quoted to three figures and vary by two decades. That is not sloppiness in the measurement — it is the honest report of a quantity that depends on the laboratory as much as on the fluid, and knowing which part is which decides what may be designed on it.

turbulence · Transition A shear layer, and the point of inflection in it. The velocity profile U = tanh y across a layer of finite thickness, with the inflection point located by searching for a sign change in the second derivative rather than by reading it off the algebra. Rayleigh's theorem says an inviscid parallel flow can only be unstable if such a point exists — a necessary condition, not a sufficient one.

A layer with a kink in it

Rayleigh proved in 1880 that an inviscid shear flow cannot be unstable unless its velocity profile bends the other way somewhere. It is one line of algebra, it is necessary and not sufficient, and it ties instability to separation through the sign of a single derivative at the wall.

turbulence · Instability The chart, with one exact line on it. The friction factor of a pipe against Reynolds number, for five relative roughnesses. Every curve here except one is Colebrook's correlation, solved by iteration rather than read off a chart. The exception is the short straight line at the left: f = 64/Re is the laminar solution and it is exact. The curves flatten to the right because once the roughness pokes out of the viscous layer the Reynolds number has nothing left to change.

The roughness a wall cannot feel

A rough pipe and a polished one carry the same flow for the same pressure over three decades of Reynolds number, and then suddenly they do not. What changed is not the pipe. It is the thickness of the film of fluid at the wall, which is the only part of the flow that can see the roughness at all.

applied · Internal flow u⁺ = (1/κ) ln y⁺ + B, integrated rather than asserted. 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 Stokes' theorem on a solved wake at Re 40. A rectangle drawn in a viscous flow that was solved on a grid. The circulation round its boundary is computed by walking the four sides and adding up the velocity along them; the vorticity inside it is computed by adding up the stored vorticity cell by cell. The two computations share no sample point and the theorem says they must agree.

Circulation is vorticity, added up

One of these two quantities is measured by walking round a loop and one by summing over the area inside it, and a theorem says they are the same number. That reconciles the site's most confusable pair — and explains how a flow with circulation can have no spin in it anywhere.

kinematics · Vorticity Where the flow over the section first reaches Mach one. Two curves running towards each other. The falling one is the section's peak suction, solved exactly in incompressible flow and then deepened by the Prandtl–Glauert factor as the free-stream Mach number rises. The rising one is the pressure coefficient at which the local flow would be exactly sonic. Where they cross is the critical Mach number, and above it there is a pocket of supersonic flow on a subsonic aeroplane.

The pocket on top of the wing

An airliner cruising at Mach 0.85 has subsonic flow almost everywhere and a patch of supersonic flow over its wing. Closing that patch takes a shock, the shock separates the layer beneath it, and the correction that predicts all of this also predicts, in the plainest terms, where it stops being true.

regimes · Mach Pressure recovery along the upper surface at 6°. Surface speed and the local Falkner–Skan pressure-gradient parameter, plotted along the upper surface from the nose. The speed peaks near the leading edge and then falls, which is the layer climbing back up to the pressure it started at, and the parameter crosses the separation value where that climb becomes too steep.

Where the straight line stops

Ideal flow will report a lift coefficient at forty degrees of incidence without complaint. What ends the lift curve is the boundary layer refusing to follow the surface, and the estimate of when that happens joins two solvers that have nothing else in common.

viscous · Separation The cliff a rough ball reaches sooner. The drag coefficient of a sphere against Reynolds number, on log axes. The smooth curve is Morrison's correlation, which is a fit to measurements and is drawn in the colour this site reserves for a borrowed claim. The other is the same curve shifted along the Reynolds axis by a factor of 6 — a stated model of what roughness does, which is to trip the boundary layer early, and not a measurement of any real ball.

The drag that falls as it speeds up

There is a band of speeds in which a smooth ball experiences less drag the faster it goes. Not a smaller coefficient — a smaller force. Dimples move that band down to where a golf ball actually flies, and they do it by making the friction worse.

applied · Sport ball Two sides of one ball, at different pressures. The surface pressure coefficient round a ball, measured from the front stagnation point, on the side the seam trips and on the side it does not. Up to separation both follow the exact potential-flow distribution 1 − (9/4)sin²θ. After it both take the same wake pressure, which is what a manometer measures rather than what the ideal theory predicts. The asymmetry is the shaded area between them, and integrating it gives a side force of 0.2949 towards the later-separating side.

A ball that swings without spinning

A cricket ball curves in flight with no spin about any useful axis. The mechanism is not the Magnus effect; it is a seam tripping the boundary layer on one side so that side lets go later. Which way the ball then goes depends on one borrowed number, and this site's own inviscid solver supplies the value that gets it wrong.

applied · Sport ball Every wavelength grows, and the shortest grows fastest. Growth rate against wavenumber for a vortex sheet of zero thickness and for a layer of thickness δ. The sheet's curve rises without limit, which is the model announcing that it has no shortest scale in it; the layer's turns over and dies at kδ = 1, the neutral wavenumber computed from Rayleigh's equation rather than fitted.

The gradient that does both

One line of the boundary-layer equations at the wall says the profile's curvature there equals the pressure gradient. That single sign causes separation and causes instability, and it causes the instability a long way before it causes the separation.

viscous · Separation Two drag laws, one derived and one fitted. Flat-plate drag coefficient against Reynolds number, laminar and turbulent, on log axes. The laminar curve is Blasius' similarity solution, solved by shooting; the turbulent one is the 1/7-power correlation, which was fitted to experiment and is drawn dashed-in-kind to keep the difference visible. Their slopes differ — −1/2 against −1/5 — so the gap between them widens rather than staying put.

The cost of going turbulent

A turbulent boundary layer costs several times the friction of a laminar one, and the multiple is not a constant — it rises with Reynolds number, because the two laws have different exponents. That is why laminar flow is worth more on a long fast surface than on a short slow one.

viscous · Drag budget 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 One more condition, and the price of it. The model problem ε u″ + u′ = 0 with a condition at each end, at three values of ε. The outer solution is the flat line at one — that is the whole of the answer when ε is zero, and it is a first-order equation that can meet one condition, so it meets the one at the far end and misses the one at the wall by the whole range. Restoring ε restores the second condition and pays for it with a layer of thickness ε, inside which the gradient is of order 1/ε. The product of those two — which is what a stress is — does not depend on ε at all: it is 1.0000 at every value tried, to nine decimal places. Drag does not vanish as viscosity does. It converges.

How many things a flow must be told

The equations of motion do not have one answer. They have as many as the conditions on the edge allow, and the number of those is decided by the order of the equation — which is why viscosity does not make the same problem harder, it makes a different problem.

kinematics · Boundary conditions A boundary layer with no x in it. The velocity profile over a porous wall with uniform suction: U(1 − e^{−Vy/ν}), exactly, at every station along the wall. The displacement thickness is ν/V, the momentum thickness is half of it, and the shape factor is two — all of them constants, none of them a function of distance. It is the cleanest demonstration there is that a boundary layer's thickness is a balance rather than an accumulation.

The layer that stops growing

Blasius' boundary layer thickens as the square root of distance and never stops. Suck fluid through the wall at a uniform rate and it stops immediately — the profile becomes a single exponential with no x anywhere in it, and the friction comes out exactly equal to the momentum of the fluid that was taken away.

viscous · Exact layer α = 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 One layer, and the four others inside and outside it. The velocity profile over a flat plate, and the temperature profile in the same layer at four Prandtl numbers: a liquid metal at 0.01, air at 0.71, water at 7 and a heavy oil at 100. The equations differ by one number and the profiles differ by a factor of twenty in thickness. At Pr = 1 the two are the same function — not similar, identical, to eight decimal places — because the equations and the conditions are then the same, which is what every statement called a Reynolds analogy rests on.

The other layer, and the one number that separates them

A wall in a stream carries two boundary conditions and grows two layers. Their thicknesses differ by a factor of twenty across ordinary fluids, and at exactly one Prandtl number the two profiles are not similar but identical.

regimes · Peclet 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 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 A slotted flap at 30°, in a flow with no viscosity anywhere. Streamlines through a main element and a flap, computed by a two-body panel solve. Each element carries its own circulation and its own Kutta condition, and the two interfere through their velocity fields and through nothing else — there is no boundary layer here, no wake, no mixing region and no high-energy air. The system's lift coefficient is 2.757 against 0.698 for the main element alone at the same incidence, and the main element itself is carrying 3.98 times the circulation it carries by itself.

A slot is not a nozzle

The gap between a wing and its flap is supposed to blow fast air into a tired boundary layer. A solver with no boundary layer in it at all — no viscosity, no wake, no mixing — produces most of the lift increment anyway, and produces it on the element nobody moved.

circulation · Slot 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 What the skin settles at, before anything is done to it. The adiabatic wall temperature against Mach number, in air at 216.7 K, with the stagnation temperature above it. The gap between the two is the recovery factor, which is 0.8417 here and stays there at every Mach number — it is a property of the Prandtl number and not of the speed. At Mach 2 the skin sits at 363 K, at Mach 3 at 545 K, and at Mach 5 at 1128 K, which is past what aluminium will do. Nothing has been burnt and nothing has been rubbed: the air was brought to rest, and this is where its kinetic energy went.

The wall that heats itself

A surface told nothing about its temperature does not settle at the air's. It settles most of the way to the stagnation temperature, and the heat flux is driven from that invented temperature rather than from the free stream's — so a wall hotter than the air can be being heated by it.

compressible · Recovery 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 Two profiles, and they are the same profile. The chordwise and spanwise velocity profiles in the boundary layer of a yawed flat plate, each as a fraction of its own edge velocity. They are computed by different code — the chordwise one by shooting a third-order nonlinear equation, the spanwise one by a single pass through a linear second-order one — and they agree to 2.1e-8 over the whole layer. They are the same function of η, because the two equations reduce to the same equation. A swept flat plate has no crossflow at any sweep angle, and that is the independence principle in the only form that has no wriggle room in it.

The wind a swept wing feels

Sweeping a wing back is usually justified by saying it meets a slower wind. It does not meet a slower wind. The equations split exactly in two, and the flow along the span is a passenger that exerts no force and changes nothing — until a pressure gradient breaks the split, and then it becomes the reason a swept wing is a different problem rather than a harder one.

circulation · Sweep 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 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 things moving by a million, and their product standing still. Viscosity, the squared velocity gradient at the dissipation scale, and their product, across six decades of Reynolds number at a fixed large-scale flow. The viscosity falls by a factor of 1e+6; the squared gradient rises by exactly the same factor, because η falls as Re^(−3/4) and u_η as Re^(−1/4); and the dissipation ν(u_η/η)² does not move at all. That is the dissipation anomaly stated as arithmetic: the limit of the dissipation as viscosity vanishes is not the value it takes when viscosity is zero.

The limit that is not the value

Dissipation is viscosity times the square of a velocity gradient, so it ought to vanish as the viscosity does. It does not. The gradient rises by exactly the factor the viscosity falls by, the product stands still, and a fluid with no viscosity at all would dissipate nothing — which is why the limit and the value are different numbers.

turbulence · Dissipation A curved shock, and the entropy each streamline picks up crossing it. A parabolic bow shock ahead of a blunt nose at Mach six, with the streamlines drawn arriving horizontally and a marker at each crossing whose size is the total pressure lost there. The streamline through the nose crosses a normal shock and keeps three per cent of its total pressure; one four nose radii out crosses at fourteen degrees and keeps ninety-four per cent. Every streamline gets a different entropy, and the stagnation enthalpy is the same on all of them.

The spin a shock leaves behind

A curved shock gives every streamline a different entropy rise and the same stagnation enthalpy. Crocco's theorem then forces vorticity into a flow with no viscosity anywhere — and it scales as the inverse of the shock's radius of curvature, exactly, so a straight shock makes none.

compressible · Crocco 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 The wall the outer flow is really solving for. A flat plate, the edge of its boundary layer, and the line the outer flow behaves as though the plate were on. The displacement thickness is the mass deficit divided by ρU — checked here against the profile's own integral rather than quoted — and moving the wall out by that much reproduces exactly the flow rate the viscous layer lets past. It is a third of the visible thickness of the layer and it is the only part of the layer the outer problem knows about.

The body the outer flow actually sees

A boundary layer lets less fluid past than an inviscid one would. The outer flow can be given exactly the same reduced flow rate by leaving the fluid inviscid and moving the wall out — so the potential flow that matters is not the flow past the body, but the flow past the body plus a thickness the boundary layer computes.

inviscid · Interaction A fifth of the drag power is not heat yet. The energy account of towing a flat plate a metre long through air at thirty metres a second. The power it takes is the drag times the speed. The heat made inside the boundary layer is half the free-stream energy times the energy thickness, and it is less — 78.6 per cent of what was paid. The rest has not been destroyed: it is kinetic energy still in the wake, which will become heat somewhere downstream, in fluid that is no longer touching the plate.

The third thickness

A boundary layer has no edge, so every thickness quoted for it is an integral of the profile against some weight. Two of them are famous. The third answers a question the other two cannot — how much of the power spent towing a plate has actually become heat by the time the fluid leaves it — and the answer is 78.6 per cent.

viscous · Energy thickness Three answers for the lift-curve slope, with three different signs. Thin-aerofoil theory has no thickness term at all and returns 2π for every section. The exact potential solution rises: 2π(1 + 0.766 t/c), measured off the Joukowski map. Real sections do the opposite, because the boundary layer thickens towards the trailing edge and decambers the section. Two of these curves are computed here; the third is what measurement says.

The half that carries nothing

Thin-aerofoil theory splits a section into a camber line that carries all the lift and a thickness distribution that carries none — at any incidence, exactly none. That is very nearly true, and what it discards decides the peak suction, the critical Mach number and where the boundary layer gives up.

circulation · Thickness Nu/√Re against the Prandtl number, over eight decades. The whole of the flat plate's heat transfer, as one curve. It is not a power law: it goes as Pr^½ at the bottom, where the thermal layer is far thicker than the viscous one, and as Pr^⅓ at the top, where it is buried inside it. The Pr^⅓ everybody quotes is the upper half. The two asymptotes are drawn beside it, and the low one is √(Pr/π) in closed form.

The number that is an answer

Almost every dimensionless group is a hypothesis: somebody sets the speed, the size and the fluid, and the number licenses a model. The Nusselt number is not. It is what the experiment produces, it sits on the left of the equals sign, and a regime diagram drawn on it is a category error.

regimes · Nusselt Stokes' drag stops being the answer almost at once. The damping on a millimetre sphere oscillating in air, divided by Stokes' steady drag on the same sphere, against frequency. The extra term is the sphere's radius over the layer thickness, so it takes over as soon as the layer is thinner than the body — which for a millimetre sphere in air is below one hertz. At a kilohertz the damping is fifteen times the steady value, and the exponent is a half rather than zero.

What a fluid takes out of a swing

The damping a body feels from the air around it is not Stokes' drag, and stops being it far earlier than anybody expects — a millimetre sphere in air is already forty-six per cent above the steady answer at one hertz. Past that the damping rises as the square root of the frequency, and the fluid it is fighting is a shell a fraction of its own size.

viscous · Damping Nu/Gr^¼ against the Prandtl number, with the exact solution's points on it. The closed form 0.508 Pr^½(20/21 + Pr)^−¼, over eight decades, with Ostrach's exact similarity values marked. The integral method is two to eight per cent high from Pr = 0.7 upwards and 27 per cent high at Pr = 0.01 — which is where the thermal layer is ten times the momentum layer and giving them one thickness stops being an approximation to anything.

A speed nobody imposed

Every regime number in this collection contains a velocity somebody chose. Natural convection has none: a warm plate makes its own flow, and the Grashof number is what is left when the speed is taken out. The Reynolds number of the result — six thousand, on an ordinary radiator — is an output of the solution rather than a setting on an apparatus.

regimes · Grashof 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 exact solution and its three approximations, at ε = 0.02. The outer solution is excellent everywhere except in a layer of width ε at the left, where it is wrong by a whole unit. The inner solution is excellent inside that layer and wrong everywhere else. The composite is their sum less the part they agree about, and it is within order ε of the exact solution across the whole interval — which is the entire content of matched asymptotics, drawn.

One formula for both ends

Two limits, each with its own description, neither valid everywhere. The composite is the sum less the part they agree about, and it is uniformly good — but the overlap region that justifies the construction does not exist at ε = 0.01, and the composite is still accurate to two per cent there.

regimes · Crossover 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 Where a rotor's pressure rise comes from, as the radius moves. The static pressure rise across a rotor, split into the two terms rothalpy gives it. The diffusion term is held at the de Haller limit throughout — the blade is being asked to slow the relative flow as hard as a boundary layer will allow — so it is a flat 11558.4 Pa at every radius ratio. Everything above that line is the centrifugal term, which costs no diffusion and has no limit of its own. At a radius ratio of 2 it supplies 49.92 per cent of the rise and at 3, 72.66 per cent. An axial machine, at a ratio of exactly one, gets none of it.

What a turning frame keeps

Euler's equation prices the work and says nothing about where the pressure comes from. In the frame turning with the blades — which is accelerating, and carries two fictitious forces — a Bernoulli-like quantity survives both of them, and it splits the pressure rise into a term a boundary layer limits and a term that is free if the radius moves.

applied · Turbomachine Four guesses at a boundary-layer profile. A straight line, a parabola, Pohlhausen's cubic and a quarter sine, each rising from zero at the wall to the free stream at the edge. Two of them also satisfy the conditions the true profile satisfies — no curvature at the wall, no slope at the edge — and two do not, which is what sorts them.

Four profiles, one drag

The momentum integral is exact and asks nothing about the shape of the velocity profile. Four guesses at that shape span twenty-three per cent in the drag they give — and the two that satisfy the conditions the true profile satisfies are within three, which says the freedom belongs to the family rather than to the constraint.

viscous · Boundary layer The eighths nobody chose. Four physical statements — the inner layer sits in the classical one's shear, its inertia balances its own viscous stress, the pressure is of the order of that inertia, and the displacement it makes produces that pressure — are a linear system in four exponents. Solving it gives three eighths, five eighths, one eighth and a quarter, exactly.

The length the limit invents

Prandtl's equations are parabolic, so nothing at one station can depend on anything downstream of it. Every experiment shows the pressure rising ahead of a shock or a step. The resolution is a region three eighths of a power of the Reynolds number long, which the limit that produced the equations was supposed to have removed.

inviscid · Interaction 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 In clean water a bubble rises nearly three times as fast as the same bubble in tap water. The terminal rise speed of an air bubble in water at 20 °C against its radius, from buoyancy balanced against drag: with a clean, shear-free surface using Moore's law, and with a surface immobilised by contamination using the rigid-sphere correlation. At 0.3 mm the clean bubble rises at 13.0 cm/s against 6.7; at 0.5 mm at 31.0 against 11.2, a factor of 2.76. Beyond a radius of 0.47 mm the clean bubble's Weber number passes one, its shape flattens, and a spherical calculation stops describing it; that region is shaded. Nothing about the bubble's size, gas or liquid changes between the two curves — only whether its surface can move.

The vorticity a clean surface cannot refuse

A clean bubble's surface cannot hold a shear stress, and it is easy to conclude that it makes no vorticity. On a curved surface it must carry exactly 2κu — three times the speed over the radius at a sphere's equator, whatever the Reynolds number. That is so much weaker than a rigid wall's that the flow stays irrotational to leading order, and the bubble's drag is the dissipation of that irrotational flow: 48/Re, three to ten times below a rigid sphere's.

kinematics · Boundary conditions Three modes carry heat up to a ceiling of three; the rolls they were cut from keep going. The Nusselt number — heat carried across the layer over what conduction alone would carry — against r, the Rayleigh number over its critical value, for Lorenz's three-mode truncation, 1 + 2(r − 1)/r, and for steady rolls at the same wavenumber computed with 6 and with 42 temperature modes. At r = 2, r = 5 and r = 30 Lorenz gives 2.000, 2.600 and 2.933; the 42-mode rolls give 2.143, 3.323 and 5.970. Lorenz's value can never exceed three whatever the Rayleigh number; the rolls' keeps rising. The ceiling is not in the convection. It is in the three modes, which have only one way to thin the thermal layers, and half of it is used by the time the layer is twice past onset.

The truncation that cannot carry three times the heat

Lorenz's three modes give a convecting layer's heat flux in one line — one plus twice (r − 1) over r — and it can never reach three times what conduction carries. The same rolls computed with forty-two modes agree with that line exactly at onset, carry 7.1 per cent more heat at twice the critical Rayleigh number, and twice as much at thirty times it. The ceiling is not in the convection; it is in having one sine to draw the temperature with.

turbulence · Convection Two external flows that agree where it matters and nowhere else. The velocity just outside the boundary layer, for two pressure distributions, against distance along the surface. They cross at the half-way station with the same value and the same gradient, and they have nothing else in common: one accelerates steadily and the other does most of its accelerating at once.

A layer that is an integral of everything upstream

Two surfaces are given external velocity distributions that agree exactly at one station — the same speed and the same gradient. The boundary layers there differ by 38 per cent in momentum thickness, and the two surfaces separate five per cent of their length apart.

viscous · Boundary layer 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 cooled wall at β = 1: the linear relation is 158 K out inside the layer. The static temperature across a laminar layer at Mach 5 with an edge temperature of 220 K, over a wall whose total enthalpy is 0.5 of the edge's, at a pressure gradient β = 1: exact (thick) and from the Crocco–Busemann linear relation (thin), at a Prandtl number of one with constant properties. The wall is at 660 K in both. The exact profile peaks at 684 K and the linear one at 759 K; the largest difference, −158.3 K, is at η = 0.81.

The gradient the heat never hears

On a flat plate at a Prandtl number of one, a boundary layer's total enthalpy is a straight-line function of its velocity, whatever the wall's temperature. Put the same layer in a pressure gradient and the straight line fails everywhere except on an insulated wall, because the gradient enters the velocity's equation and not the enthalpy's — and a favourable gradient can leave a band of gas colder than the free stream above a wall three times hotter than it.

compressible · Recovery 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

Named alongside it

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

SeparationModel limitSimilarity solutionSkin frictionReynolds numberViscosityWall shearBoundary conditionTransitionCorrelationFalkner–SkanAdverse pressure gradient

All concepts