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

Navier–Stokes equations — where it appears

The momentum equations for a Newtonian fluid, complete, deterministic and known since 1845. Whether their three-dimensional solutions stay smooth for all time is unsettled, and the difficulty of solving them has nothing to do with anything missing from them.

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

Ten unknowns, four equations. What is left after the Navier–Stokes equations are averaged. The mean velocities and mean pressure were there before; the six Reynolds stresses are new, and they arrived from the one term that does not average away. Nothing in the count is an approximation — the averaged equations are exact — and that is what makes the gap uncomfortable.

What averaging costs

Split the velocity into a mean and a fluctuation, average the equations, and the result is exact. It is also short of six equations, because the one nonlinear term does not average away and leaves six new unknowns behind that nothing determines.

turbulence · Closure The one place the stretching argument closes. Burgers' vortex: an axisymmetric strain carrying vorticity inwards at exactly the rate viscosity spreads it outwards. The vorticity profile is a Gaussian of radius √(4ν/α), the swirl velocity peaks at 1.12 core radii rather than at the core radius itself, and the circulation reaches its full value by about two. The steady vorticity equation is evaluated on this profile by differencing it, not by re-deriving it.

The spin that feeds itself

Stretch a vortex tube and its spin rises in exact proportion, because the circulation round it cannot change and its area has fallen. Nothing in that argument sets a limit — and the one flow where the limit can be written down exactly puts it at a length of √(4ν/α).

kinematics · Vorticity 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 The pressure, relaxed rather than quoted. The pressure field round a cylinder, obtained by relaxing ∇²p = −ρ∇·(u·∇u) on a body-fitted polar lattice with the exact pressure on the surface and on a far circle. The interior was told nothing except the velocity gradients. It agrees with Bernoulli's closed form everywhere to under two parts in ten thousand of the dynamic pressure, which is the strongest statement this site can make that the elliptic equation is the pressure's own.

Pressure has no speed

Take the divergence of the momentum equation for an incompressible flow and the time derivative disappears, the viscosity disappears, and what is left is Poisson's equation. Pressure is not carried anywhere: it is whatever satisfies an elliptic equation everywhere at once, and that is a statement about a fluid nobody has.

inviscid · Pressure A cylinder in a uniform shear, K = 0.4. A stream whose velocity increases with height, meeting a circular cylinder. The oncoming profile is drawn at the left. The flow carries uniform vorticity −K, so it is a solution of Euler's equations and not of Laplace's, the pattern is no longer symmetric top to bottom, and the body feels a lift towards the fast side with no circulation anywhere.

Inviscid does not mean irrotational

Dropping viscosity gives Euler's equations. Assuming nothing is spinning gives Laplace's — one scalar, linear, unique. The second step is a separate hypothesis about the flow's history, and a flow that fails it is still an inviscid flow with exact solutions of its own.

inviscid · Euler rotational Seven shock structures, and the two states they all connect. The velocity through the shock for seven dissipation models — Prandtl numbers from a quarter to two, viscosities from constant to linear in temperature. Each curve is shifted so its midpoint sits at the origin. They start at the same speed, end at the same speed, and are nothing alike in between.

The jump does not ask what made it

Seven different dissipation mechanisms are made to smear the same shock. Their interiors are a factor of two and a third apart in thickness, their entropies overshoot the final value by between a quarter and a doubling, and the state they all arrive at agrees to seven parts in ten billion — because the end states are conservation and the interior is transport.

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

Named alongside it

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

Exact solutionBoundary conditionNonlinearityBoundary layerLaminar flowModel limitSimilarity solutionViscosityViscous diffusionBernoulli's equationKelvin's circulation theoremMeasurement

All concepts