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

Heat transfer — where it appears

The flux of thermal energy between a surface and the fluid past it, carried by the same layer that carries momentum. That shared mechanism is what ties the heat transfer to the friction, and the tie is exact when the Prandtl number is one.

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

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 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 number a duct settles at is an eigenvalue. The local Nusselt number against x⁺ = x/(D·Re·Pr), from a Crank–Nicolson march that contains no eigenvalue anywhere. It settles at 3.6568, which is λ₀²/2 for the Graetz eigenvalue problem — a completely separate calculation. The mark at x⁺ = 0.05 is the entry length every textbook quotes: it delivers a Nusselt number 1.45 per cent above the developed value, which is a perfectly reasonable tolerance and is never the one stated.

How far before the heat arrives

A duct's thermal entry length is quoted everywhere as x/(D·Re·Pr) = 0.05, with no tolerance attached. Working out what it delivers gives a Nusselt number 1.45 per cent above the developed value — and the developed value itself is not a term ratio at all but an eigenvalue, 3.6568, which is also the rate at which the duct forgets its inlet.

regimes · Peclet 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 temperature behind a shock, which is not the jump condition's. The static temperature along the flow behind a Mach 6 normal shock, with the frozen value the jump conditions give and the equilibrium value they give a long way behind. The gas arrives at 2382 K and settles at 2059, over about four tenths of a millimetre.

A gas that has not finished being shocked

The jump conditions give the state a long way behind a shock. Immediately behind it the molecules have not started vibrating yet, so the temperature is 2,382 K where the equilibrium answer is 2,059 — and the gas takes four tenths of a millimetre to get from one to the other.

compressible · Real gas A skin lags the air by a time its own thickness sets. A skin held at one flight condition, relaxing towards the adiabatic wall temperature. The approach is exponential with a time constant ρc τ / h — 8.37 s at 2 mm, 25.1 s at 6 mm, 50.2 s at 12 mm — so the temperature a steady recovery calculation gives is reached after several minutes rather than at once. The fluid supplies one number to this calculation, the adiabatic wall temperature, and the structure supplies everything else.

The skin that lags the flight

A wall can be told its temperature or told nothing, and both are solved problems. A real skin is told neither. It has heat capacity, so its temperature is a transient whose time constant is its own thickness divided by what the layer delivers — and the number the steady calculation returns is an upper bound a short exposure never collects.

compressible · Recovery Measuring the air temperature in flight is not measuring the air temperature. Three temperatures against Mach number in a stream at 220 K: the free stream itself, the stagnation temperature a perfect probe would read, and what a probe recovering 0.98 of the rise actually reads. At Mach 0.85 the rise is 31.79 K and the probe misses 0.64 K of it. The instrument is a small stagnation region with a thermocouple in it, so it obeys the same recovery arithmetic a wall does — a probe is a wall told nothing about its temperature, made small and put on a stalk, and its recovery factor is a calibration constant rather than a one.

The thermometer that heats itself

A total-temperature probe is a wall told nothing about its temperature, made small and put on a stalk, so it obeys the same recovery arithmetic an aircraft skin does. Its recovery factor is a calibration constant near 0.98 rather than a one — and the static temperature inferred from its reading amplifies that shortfall rather than inheriting it.

compressible · Recovery 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 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 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 A pipe at the gas's temperature does not keep the gas at it. Left, the static temperature of the gas along a pipe, as a fraction of the wall's, against the local Mach number, entered at Mach 0.1: with no heat transfer, with the heat transfer Reynolds' analogy gives, and with ten thousand times that. Right, the stagnation temperature. The isothermal model holds the static temperature at the wall's and needs the stagnation temperature to climb by 14 per cent; every real case keeps it within two per cent, cools, and runs on to Mach one.

A pipe cannot hold its gas at the wall's temperature

The textbook model of a long gas pipe in contact with the ground holds the gas at the ground's temperature and has it choke at 0.845 of the speed of sound. No pipe does either. A wall at the gas's temperature draws heat out of it rather than putting heat in, and with any strength of heat transfer at all the flow runs on to Mach one, within a tenth of a per cent of the length a perfectly insulated pipe would need.

compressible · Fanno rayleigh Where a carburettor's cold comes from. The temperature changes in a carburettor's throat, in kelvin: the cooling of the air as it expands to a hundred metres a second, the far smaller cooling of the metal it flows over, the cooling as petrol evaporates into it, and the warming as the water in humid air condenses back out. The fuel does five times what the venturi does to the air and fifty times what it does to the metal.

Carburettor ice is made by the fuel, not the venturi

Pilots are taught that a carburettor ices because air cools as it expands through the venturi. The expansion does cool the air, by five kelvin at a hundred metres a second, but the metal the ice grows on barely feels it. The cold comes from the petrol evaporating, and whether it makes ice is decided by the water in the air, which warms the mixture as it condenses. The worst day is not a cold one.

misconceptions · Venturi A fuel-wetted plate sits far below the stream that cools it. The temperature of a throttle plate wetted by petrol, against the day's temperature at 80 per cent humidity, beside the temperature of the stream flowing over it — the mixture after half the fuel has evaporated into it, at the plate's recovery temperature. A plate kept wet with fresh fuel sits nine to twenty-three kelvin below the stream, more on warmer days, and below freezing on any day up to about 34 °C. A plate whose film has lost its light ends sits one to six kelvin below it. The dashed line is freezing.

The throttle plate is a wet bulb for petrol

A carburettor's icing is usually worked out for the mixture — the fuel's latent heat spread through all the air — and that puts the warmest icing day near 15 °C, well short of the 30 °C the charts pilots use allow for. The ice does not grow in the mixture. It grows on the throttle plate, and a plate wet with evaporating petrol is a wet-bulb thermometer for fuel: nine to twenty-three kelvin colder than the stream flowing over it. Kept wet with fresh fuel it ices on a 30 °C day at full humidity and a 40 °C day at half — which is what the training handbook warns of. Once its light ends have gone it does not, and that difference is most of the story.

misconceptions · Venturi 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 A cold wall holds its layer on longer; a hot one lets it go sooner. The pressure-gradient parameter β at which the wall shear vanishes, against the wall's total enthalpy over the edge's. With the velocity deaf to the temperature it is −0.1988 for every wall. Coupled, it is -0.3264 for a wall at absolute zero, -0.2623 for a wall at half the edge's total enthalpy, -0.1988 at one, -0.1573 at one and a half and -0.1295 at two. Cooling delays separation by more than the same heating hastens it.

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 Each collapsed form of the plate's heat transfer reaches its limit decades away from air. The exact local Nusselt number of a vertical isothermal plate in its two collapsed forms against the Prandtl number: over (Gr·Pr²)^¼, which tends to Le Fevre's 0.600 as Pr → 0, and over (Gr·Pr)^¼, which tends to 0.503 as Pr → ∞. At Pr = 10⁻⁵ the first reads 0.5994; at 10⁵ the second reads 0.5023. The two limit laws cross at Pr = 0.494, where the exact value is 0.736 of either; air, at 0.72, reads 0.421 and 0.387.

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 A flat plate's layer is pushed along by the pressure its own warmth leaves, strongest at the wall. Across the layer of a horizontal plate warm side up, in air: the velocity scaled to its peak, the temperature excess, and the along-plate push −∂p′/∂x scaled to its wall value, against the similarity distance η = (y/x)·Gr^⅕. The push is the hydrostatic deficit under the warm fluid, (2/5)(ηθ − P), 0.694 at the wall in units of gβΔT·Gr^−⅕, falling to nothing at the layer's edge; the velocity peaks at η = 1.58 and the warm layer ends near 5.93.

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 Of the four common plumes, only the laminar one above a line speeds up as it rises. Centreline speed against height for four plumes, each scaled to its own speed at one reference height: laminar above a line source, rising as the fifth root of height; laminar above a point source and turbulent above a line, both constant; and turbulent above a point, falling as height to the minus one third. Each exponent follows from the heat flux being fixed and the width growing by diffusion, √(νx/u), when laminar, or in proportion to height when turbulent.

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 At a levitator's own thermal numbers every gas lifts within a few per cent of the isothermal film; only a thick, slow film lifts more. The mean pressure excess of a trapped film squeezed by a quarter of its gap, over the isothermal film's 0.0328 of ambient, against the thermal number ωh²/κ, for three ratios of specific heats; dots, the six gases of a 20 kHz, 20 µm levitator at one atmosphere: helium 1 at Ω = 0.29; air 1.002 at Ω = 2.4; argon 1.005 at Ω = 2.5; carbon dioxide 1.005 at Ω = 4.7; xenon 1.042 at Ω = 7.9; sulphur hexafluoride 1.003 at Ω = 16. At Ω = 300 a monatomic gas lifts 1.21 times the isothermal film, against the 1.69 a core on the walls' mean temperature would give.

A levitator lifts by its amplitude, not its gas

A squeeze-film levitator holds its load up by the mean of a pressure that swings symmetrically about nothing, and the gas in the gap was expected to decide how much: a gas that cannot shed its heat of compression in time is stiffer, so it should lift more. Solved to every order in the amplitude, with the heat crossing the gap, it barely does. At a levitator's own thermal numbers every common gas lifts within a few per cent of the isothermal film while its losses differ several-fold, and even a gas far too slow to shed its heat lifts only a fifth more, because the walls quietly pump heat out of it and cool its middle.

viscous · Squeeze film

Named alongside it

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

Model limitPrandtl numberBoundary layerSimilarity solutionRecovery factorAdiabatic wallMach numberMeasurementStagnation temperatureAerodynamic heatingBoundary conditionBuoyancy

All concepts