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

Similarity solution — where it appears

A solution that has the same shape at every station once the coordinates are stretched by the right power. The stretching turns a partial differential equation into an ordinary one, which is why so many exact results in viscous flow are of this kind.

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

The laminar boundary-layer profile. Speed against height through a laminar boundary layer on a flat plate, in the similarity variable that collapses every station along the plate onto one curve. The straight line is the slope at the wall, which is what the skin friction is proportional to.

How thick is thin

The boundary layer has no edge. It approaches the free stream and never arrives, so any thickness quoted for it is a convention — and the three conventions in use measure three different things, one of which is not a height at all.

viscous · Boundary layer The wall slope against pressure gradient, and where it runs out. How steeply the flow leaves the wall, plotted against the pressure gradient the layer is running into. A favourable gradient presses the profile against the surface and steepens it; an adverse one hollows it out. The curve reaches zero at a definite value, and beyond that there is no attached solution at all.

How much uphill a layer can take

A boundary layer running into rising pressure is climbing a hill on the last of its momentum. There is a definite steepness at which it can no longer do it, and the number is not a rule of thumb — it is where a family of solutions stops existing.

viscous · Separation 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 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 45° of corner, and an infinite number of eddies in it. The creeping flow in a corner of 45 degrees, drawn from Moffatt's similarity solution. Each eddy turns the opposite way to its neighbours and is 3.17 times smaller and 1.6e+3 times weaker than the one outside it. The contour levels are rescaled inside each eddy, because they differ in strength by three orders of magnitude per step and a single set of levels would show the first and nothing else — which is itself the reason nobody has seen the third. The dividing lines between eddies are drawn where the stream function changes sign, and the dots are the centres, both found from the solved field.

The eddies nobody stirs

A slow flow past the mouth of a sharp corner does not simply fail to get in. The corner fills with an infinite sequence of counter-rotating eddies, each about three times smaller and sixteen hundred times weaker than the last, and the whole structure is decided by one complex number.

viscous · Corner 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 core spreads and the outside never notices. The swirl velocity at four times a factor of four apart, with the free vortex Γ/2πr drawn behind them. Every curve leaves the free-vortex line at its own core radius and turns over into solid-body rotation inside it; outside the core all four are the same curve, to the precision of the plot. Viscosity has rounded off the singularity and changed nothing else. The peak swirl falls from 52.4 to 6.6 m/s across the four, and the circulation is identical for all of them.

What viscosity cannot take away

Leave a vortex alone in a viscous fluid and every local measure of it falls — the peak spin, the peak velocity, the enstrophy. The circulation round a large loop does not move at all, ever, and the far field is identical to the line vortex it started as.

viscous · Diffusion 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 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 The shock radius against time, from an equation that was given no exponent. The thin-shell energy balance integrated forward from a small initial radius, on logarithmic axes. The two-fifths power is not put in: the ordinary differential equation is Ṙ = √(E/AρR³), and the straight line is what it does. The fitted slope is 0.39983 and the fitted prefactor is 0.90721 against a closed form of 0.90702 — the small residuals being the integration's memory of where it started, which the similarity solution has no equivalent of.

A radius that gives the energy away

Four quantities, three dimensions, one group. The radius of a strong blast must be a constant times (Et²/ρ)¹⁄⁵, and nothing about the device, the chemistry or the initial size can appear. The exponent is free and the constant is not.

compressible · Blast wave 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 A vortex has no energy. The kinetic energy of a Lamb–Oseen vortex inside a circle, per metre of its length, against the logarithm of that circle's radius. It is a straight line and it does not stop: outside the core the swirl is Γ/2πr, the energy density falls as 1/r², and the area grows as r², so every decade of radius adds the same amount. There is no such thing as the energy of a line vortex without a stated cutoff, and no cutoff is physical.

The energy a vortex cannot have

A line vortex has infinite kinetic energy. Not a large amount — infinite, growing without limit as the logarithm of however far out the counting stops. And it is losing that energy at a rate that is finite, exactly known, and contains no cutoff at all.

viscous · Diffusion 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 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 10 per cent low 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 convergence exponent depends on the gas, which a dimensional exponent cannot. R ∝ (−t)^α for a converging shock, against the ratio of specific heats, for cylindrical and spherical symmetry. Guderley's exact values are marked and the agreement is to four figures. The Sedov blast's two-fifths is drawn beside them: it is the same for every gas, because it comes from dimensions and a conserved energy, and γ is dimensionless.

An exponent dimensions cannot give

A blast wave's radius goes as the two-fifths power of time, and the two-fifths is arithmetic: count the dimensions and it falls out. A shock converging on a point goes as the 0.717 power, and no amount of counting will produce that number — because it depends on the gas, and γ is dimensionless.

regimes · Similarity 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 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 The profile a diverging channel flattens into, and then cannot hold. Five purely outward profiles in a wedge of 0.2 radians, at rising flux, each normalised to its own centreline value. As the flux rises the profile flattens in the middle and steepens at the walls — and then it stops. The last one has zero slope at the wall, which is separation, and beyond it no purely outward profile of this form exists at all. Nothing was added to the equation to make that happen: the wall shear is the square root of a cubic and the cubic runs out.

One channel, one flux, two flows

Flow between two plane walls meeting at a line has an exact solution. Past a threshold that turns out to be a ratio of gamma functions, it has two — the same wedge carrying the same flux, once outward everywhere and once with the fluid running backwards along both walls, and nothing in the equations chooses.

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 Inside a strong blast: a shell of gas around a hot, nearly empty core. Left, the density, velocity and pressure inside a spherical strong blast wave at γ = 1.4, each as a fraction of its value just behind the shock, against the distance from the centre as a fraction of the shock radius. The density falls to nothing well inside the shock while the pressure levels off at 0.37 of its post-shock value. Right, the temperature, which is the pressure over the density and so rises without limit towards the centre.

The fireball is hollow

Inside a strong blast wave there is almost nothing. Half the air the shock has swept up lies in the outer four per cent of its radius; at half the radius the density is a hundredth of the air's; the centre is empty and, formally, infinitely hot. Integrating Sedov's equations to see this also shows that a published table of his constants was wrong at three of its four entries.

compressible · Blast wave A bow shock that grows as a line explosion does. The radius of the bow shock around a hemisphere-nosed cylinder, in body diameters, against distance behind the nose: the blast-wave analogy, R/d = 0.795 C^¼ (x/d)^½ with C, the drag coefficient, 0.919, which contains no Mach number, and Billig's correlation of measured bow shocks on spheres at Mach 5, 10 and 20, anchored at the nose and continued as a hyperbola to the Mach cone. Both grow as the square root of the distance; at Mach 20 the analogy's shock is a steady 0.71 of Billig's.

A hypersonic body leaves a line explosion behind it

A blunt body at hypersonic speed does work on the air at a rate equal to its drag, and each slice of air it passes through is struck once and left to expand. Seen from the ground, that is a line explosion, and Sedov's cylindrical blast wave gives the bow shock's width and the pressure on the afterbody without a Mach number in either. The analogy's classical constants come straight out of the blast solution. So do its limits: it holds only while its own shock stays strong, over a length that grows as the square of the flight Mach number, and it puts the body inside a core hotter than anything the flow can reach.

compressible · Blast wave Next to the body the gas remembers the nose. The gas temperature, over the free stream's, against distance from the axis of a hemisphere-nosed cylinder at Mach 15, at 2, 20 and 200 diameters behind the nose; the body's surface is at a half. The gas next to the body crossed the nearly normal part of the bow shock and carries its entropy: 16.7 times the free stream's temperature at 2 diameters, 9.29 at 200, falling outwards to the gas that crossed the weaker, oblique shock. The blast-wave analogy's core at 2 diameters, dashed, runs off the top of the frame on its way to infinity.

The gas beside a hypersonic body remembers its nose

The blast-wave analogy gets a blunt body's bow shock and afterbody pressure right and its temperature absurdly wrong: it puts the body in the empty core of an explosion, where the temperature has no bound. The real gas beside the body crossed the nearly normal shock at the nose and carries that crossing's entropy all the way down. Expanded to the afterbody's pressure, it is ten to twenty times the free stream's temperature at Mach 15, it never falls below a floor set by the nose alone, and the sheath that carries it is wider than the body.

compressible · Blast wave

The sheath outlasts the body

A blunt hypersonic body wraps itself in a sheath of hot, thin gas from its nose's shock, and its boundary layer grows by eating that gas from the inside. The usual picture has the boundary layer through the sheath within a few nose diameters. It is not: a laminar boundary layer takes about a third of a Reynolds number of diameters to swallow the sheath at Mach 15, which is thousands to millions of diameters, and a turbulent one hundreds. On any body of ordinary length the boundary layer never sees the cooler gas outside, and it is heated by the sheath's gas all the way down.

compressible · Blast wave

A cold wall holds its layer on longer

A compressible boundary layer's pressure gradient pushes on its density, and its density is its temperature. Solve the velocity and the enthalpy together and the wall's temperature reaches the velocity: a cold wall holds its layer on to an adverse gradient two-thirds steeper, a hot one lets it go a third sooner, the stagnation point's heat-to-friction ratio more than doubles between them, and the cold band an uncoupled layer put above a cooled nose disappears.

compressible · Recovery

A point force makes a jet only when it is strong

Push on a fluid at one point and there is an exact solution of the full Navier–Stokes equations for what follows, at any strength. Weak, it is Stokes's point force: fluid pushed forward over a whole hemisphere and drawn in behind. Strong, it is Schlichting's slender jet, with fluid drawn in from every direction outside a narrowing cone. One constant joins them, and it says how strong a push must be before the boundary-layer jet everyone uses is right — at the axis by a jet Reynolds number of a hundred, at the edges only by a thousand.

viscous · Free shear

A light hull slams softly, and a round one does not

A hull striking water slows as it sets the water moving, and the two share one momentum. For a V-shaped bottom that caps the slam: the force peaks when the water recruited weighs a fifth of the hull, and a hull as heavy as the water under its whole width feels a quarter of the force a constant-speed calculation predicts. For a round bottom it caps nothing, because the peak comes at the instant of contact, before the hull has slowed at all.

inviscid · Unsteady bernoulli

A jet along a wall keeps the flux of its momentum flux

A jet in open fluid keeps its momentum and gathers fluid as it goes. Blow the same jet along a wall and the wall takes its momentum away as drag, so it keeps neither. What it keeps is stranger: the flux, through each height, of the momentum flux carried above that height. That one invariant fixes everything else — the jet thickens as the three-quarter power of distance, slows as the inverse square root, and gathers fluid as the fourth root, more slowly than the free jet's cube root.

viscous · Free shear

Over a plate in mercury the warm layer flows as if inviscid

The exact solution for a warm vertical plate has been tabulated since 1953 and is usually shot from the wall, which fails in a liquid metal unless the shot starts within a thousandth of the answer. Posed instead as one algebraic system, it solves from a Prandtl number of a hundred-thousandth to a hundred thousand — and shows that in mercury the layer's speed is set without viscosity, that the integral method runs low there rather than high, and that air sits where neither limiting law is within a fifth.

regimes · Grashof

Laid flat, a warm plate is pushed by the pressure of its own layer

A warm vertical plate drives its layer with buoyancy acting along it. Lay the plate flat and buoyancy acts across the layer instead, so the only push along it is the pressure deficit its own warm fluid leaves under the cold — deeper where the layer is thicker. That push gives a heat transfer growing as the fifth root of the Grashof number, not the fourth, works only on the plate's upper face, and survives a tilt of a degree for a third of a metre.

regimes · Grashof

A jet striking a plate leaves a stronger wall jet the further it travelled

A plane jet that strikes a wall splits into two wall jets, and each keeps not the jet's momentum but Glauert's flux of momentum flux. Turned without loss at the stagnation point, each half arrives laid along the wall with its axis on the wall, and its invariant is the momentum it brought times the fluid it had gathered on the way: F = 3KQ/32. Hold the nozzle further away and the jet arrives slower but leaves a stronger wall jet, by the cube root of the distance — and the wall jet takes a distance of about the Reynolds number in jet widths to forget it was ever free.

viscous · Free shear

Above a warm wire the plume speeds up as it rises

A plume above a line of heat has no wall, no imposed speed and one conserved quantity, the heat it carries. That alone fixes a centreline speed growing as the fifth root of height — the only one of the four common plumes that accelerates — a temperature excess falling as the three-fifths power, and a warm core that is narrower than the moving sheet around it in any fluid more viscous than five ninths of its own diffusivity. At that Prandtl number the equations have a closed form, and the solver lands on it to three parts in a million.

regimes · Grashof

A scraper held flat is pushed off the wall it scrapes

A corner stirred from far away fills with an endless sequence of eddies. Drive the corner at the corner instead — let one of its walls slide under a still blade, as paint slides under a scraper — and the eddies vanish: every streamline turns once, and the stress rises into the corner without limit. G. I. Taylor solved it in 1962. The paint pushes the blade away from the wall with a force that grows as one over the square of the blade's angle, so a blade held flat rides up on what it was meant to remove, and the force is finite only because something too small to appear in the problem ends the corner.

viscous · Corner

A round jet hands its wall jet a strength proportional to its throw

A round jet striking a plate spreads into a radial wall jet, and like the plane jet before it, the wall jet keeps not the jet's momentum but Glauert's flux of momentum flux — the momentum the jet brought times the fluid it gathered on the way. A round jet gathers fluid at a rate set by the viscosity alone, eight π ν per metre of travel, whatever its strength. So the invariant it hands on grows in exact proportion to the distance it was thrown, where a plane jet's grew as the cube root: hold the nozzle four times further away and the wall jet is four times as strong and runs twice as fast.

viscous · Free shear

A bubble is throttled by its age, and a sheet by its length

A cavity cools the liquid it boils, and the cooling lowers the vapour pressure it needed. Whether that matters was a question about one number, the B-factor — vapour made over liquid cooled — and no real cavity's was known. Computed, the two kinds a pump has are nothing alike. A travelling bubble cools a thin shell that thins further as it stretches, so its B grows as the square root of its age and reaches the critical value at Brennen's thermal time: a millisecond in water at 60 °C, four microseconds at the boiling point. An attached sheet keeps renewing the liquid it cools, so its B is a few, fixed by its length and speed. In warm water the bubbles are suppressed forty degrees before the sheet.

applied · Cavitation

Named alongside it

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

Model limitBoundary layerEntrainmentHeat transferMeasurementMomentum fluxPrandtl numberSelf-similarViscosityExact solutionFalkner–SkanScaling

All concepts