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

Every essay — page 31

Page 31 of 61, continuing through the fields in the same order.

Flows and fields Ideal flow Circulation and lift Viscosity Regimes and numbers Compressible flow Transition and turbulence Fluids at work What is taught wrongly Series Concepts Regimes Refutations Search

Flows and fields

Streamlines, particle paths and the field that carries them. What is conserved, what a picture of a flow can show, and what it cannot.

4 nodes and 2 saddles, and the difference is the shape. The skin-friction pattern on a body, drawn as the streaks a film of oil would leave, looking at the nose. Every critical point of the field is marked: 4 nodes, where the streaks converge or diverge, and 2 saddles, where two streaks cross. Their indices sum to 2, which is the Euler characteristic of a sphere, and nothing about the flow enters that number — the stream is at 0 and the crossflow at 0, and neither may change it. What the flow decides is where the points are and what they look like; what the body decides is how many there are of each.

The count computed on a body

The rule for a closed surface is usually quoted and seldom solved for. This solves for one, and what the computation adds is not confirmation — it is the discovery that the count survives two events in which the number of stagnation points falls, and that both of them have closed forms.

7 figures Four kinds of critical point, and the curve that separates them. The invariants of a trace-free velocity gradient, with the discriminant curve 27R²/4 + Q³ = 0 drawn through them. Inside the two upper lobes the cubic has one real root and a complex pair, which is a spiral being stretched along its own axis on the left and squeezed on the right; below the curve all three roots are real and the point is a node with two saddle directions. Of 820 random incompressible gradients, 509 land in the spiral region and 311 in the real one. A plane flow is the vertical line R = 0 and nothing else, which is why a plane has two kinds and space has four. The marked points are the cases the calculation checks that fall inside this window; the two vortex cases it also checks sit at Q = 3.25 and |R| = 4.25, off the top corners, because a window wide enough to hold them would flatten the curve the figure is about.

Two kinds is a plane flow's privilege

A plane incompressible flow has a saddle or a centre and nothing else, and the proof is one line about a trace. The same line in three dimensions constrains three numbers instead of two, which is far less, and what it leaves is four kinds of point separated by a curve — with the one a plane cannot have being the structure the whole of turbulence is made from.

6 figures The line the dye actually draws. A streakline in an oscillating uniform stream: everything released from the origin over the last 3.4 units of time, drawn where it has got to at one instant. The parameter along it is release time, not distance and not time of flight, which is why it is a record rather than a curve of the field. The mean speed of 1 exceeds the amplitude of 0.6, so the fluid never reverses and the filament is single-valued in x.

The line the dye actually draws

The streakline is the curve most photographs really show, and of the three curves it is the hardest to compute, because it needs the whole history of the flow. This draws it, in a flow where all four curves have closed forms, and the third curve turns out to be a different kind of object from the other two rather than a third example of the same one.

6 figures 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.

6 figures 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.

7 figures Slip follows the stripes more closely than the shear does. Plan views of a striped surface, with the stripes running across each panel, for a shear at 0°, 30°, 54.7° and 90° to them. The faint arrow is the direction of the shear; the dark one is the slip velocity it produces, whose component along the stripes is the along-stripe slip length times the shear and whose component across them is half that. The slip is turned towards the stripes by 0.0°, 13.9°, 19.5° and 0.0°. It is largest, 19.47°, for a shear at 54.74°, where tan θ = √2. A surface with a tensor for a boundary condition can push a flow sideways, which a scalar slip length never can.

Twice as slippery along as across

A surface of alternating gas and solid stripes lets a liquid slip, and a flow far above it sees one number in place of the pattern — but the number depends on which way the flow goes. Along the stripes it is Philip's logarithm; across them it is exactly half, for a reason that takes one substitution to show. And the logarithm means that the slip is bought by the pattern's period rather than by how much of it is gas.

7 figures A spun cylinder carries its whole circulation at once, and hides it until the vorticity has left. The circulation round circles of radius r about a cylinder of radius a started spinning at once, as a share of the circulation of its own surface, 2πa²Ω, against r/a on a logarithmic axis, at νt/a² = 0.01, 0.1, 1, 10 and 100. At the surface it is the whole of it from the first instant, because no slip makes the fluid there turn with the cylinder. Just outside, the spin-up has laid down an equal and opposite ring of vorticity, so the circulation round a larger circle is only what has diffused past it: at two radii 0.000, 0.035, 0.611, 0.936 and 0.993 of the surface's at the same five times. The circulation a Magnus rotor needs is in the fluid the moment it spins; the far field learns of it only as fast as the counter-vorticity moves out.

A wall puts in exactly its own speed

A wall sliding in its own plane makes vorticity at a rate equal to its acceleration, with no viscosity in the rate. So however a wall is started, the vorticity it has put into the fluid is its speed, to the last digit; a wall that stops takes all of it back and leaves the fluid moving; and a spinning cylinder carries its whole circulation from the first instant, hidden behind an equal and opposite ring until viscosity carries the ring away.

8 figures A pair of points meets on R = 0, the one line an index can change on. The paths of the ABC flow's stagnation points across the (R, Q) diagram as C rises from 0.3 towards √2 with A = B = 1. The four points of index +1 share one path on the left and the four of index −1 its mirror image on the right. At C = 0.3 they sit at R = −0.088, Q = −1.045; at C = 1 they touch the discriminant curve at R = −0.707, Q = −1.500, where the strain has a repeated rate, and turn away from it without crossing; and as C approaches √2 they run in to R = -7.5e-3, Q = −2.000. Crossing into a lobe would have changed a node into a focus, which a Beltrami flow's stagnation point cannot be; reaching R = 0 is where each meets a partner of the other index.

The sign a stagnation point carries in space

In three dimensions a stagnation point's index is the sign of one determinant, and that determinant is minus the R of the invariant diagram — so the diagram's left and right halves are the two indices. An exact Euler flow in a periodic box has eight such points, four of each sign, never a spiral among them, and they can only disappear in pairs that meet on the one line where the sign is allowed to change.

7 figures A float under swell on a rotating planet goes round instead of away. The track of a float at the surface over 1 inertial periods (16.9 hours) after a 8 s swell of amplitude 1 m arrives at latitude 45°, in kilometres, the waves travelling to the right. Without friction the float runs round a circle of radius Uₛ/f = 0.479 km and comes back to where it started every 16.92 hours. With a drag on the Eulerian current it spirals out into a steady drift veered to the right: 24.3 per cent of the drift at 76.0° for a drag of 0.25 f; 70.7 per cent of the drift at 45.0° for a drag of 1 f; 94.9 per cent of the drift at 18.4° for a drag of 3 f. In a non-rotating ocean the same float would have gone 3.0 km straight on.

The drift a rotating planet takes back

In a wave tank the Stokes drift is cancelled by a return current because the tank has walls. The open ocean has none, and the drift is cancelled anyway: the Coriolis force acts on the water's real motion, drives an Eulerian current that answers it, and leaves the depth-integrated transport exactly zero at every viscosity. A float under steady swell with nothing to stop it goes round a circle instead of away.

7 figures The column in height and time: a falling interface, a rising shock, a fan. A batch settling test from a uniform φ₀ = 0.1, height above the bottom against time, both scaled on the column height and the single-particle settling time. The interface with clear water (thick) falls in a straight line at 0.4538; the sediment shock rises from the bottom at 0.1484 until the two meet at t = 1.661, height 0.2464; the thin lines are characteristics of the fan, each carrying one concentration between 0.317 and packing, and the interface bends as it crosses them. Dots are the finite-volume solve on 400 cells: the interface and the sediment front.

The column the chord rule cannot settle

A suspension settling in a closed column is a kinematic wave, and its flux curve bends both ways. At the top the chord rule works: clear water meets the suspension at a single falling front. At the bottom it does not, and the bed grows behind a shock that stops short of packing and a graded layer beneath it — so the interface, instead of arriving, slows for ever.

7 figures In the frame of the wave the walls stand still, and a bolus rides between them. Streamlines of a peristaltic channel of amplitude ratio 0.7 over two wavelengths, drawn in the frame moving with the wave, where the flow is steady and the walls are themselves streamlines. The time-mean flow is Θ = 0.5904 of the wave speed times the mean half-width, so the flow rate between centreline and wall in this frame is q = −0.4096 and the pressure rise per wavelength is 0.000 in units of μcλ/a². The centreline velocity changes sign at 0.106π and 0.894π, and the streamline through those points closes round a bolus holding 30.5 per cent of the fluid in each wavelength, which travels with the wave.

A wave on the wall is a pump

A channel whose wall only moves in and out, in a wave travelling along it, delivers a steady net flow with no part of the wall moving along the channel. In the frame of the wave the walls stand still and are streamlines, so continuity alone fixes how the laboratory flow rate follows the wall shape — and the momentum equation is needed only for one number, which also decides whether fluid rides along with the wave or leaks back against it.

7 figures Rolls turning side by side, with the fastest downwind water where they sink. The fastest-growing mode at a Langmuir number of 0.13, looking downwind, over two roll spacings of 2.89 decay depths and 4 decay depths down. The closed curves are streamlines of the overturning; the dashed curves are contours of the downwind velocity the rolls carry, positive under the lines where the water sinks. At the surface the cross-wind flow converges onto those lines, which is where floating foam and weed collect as windrows. The amplitude is arbitrary, as in any linear mode.

The drift that turns a current into rolls

A current carrying a Stokes drift feels a force the drift makes out of the current's own vorticity, and under a wind that force is unstable. It turns the surface layer into rolls lined up downwind, with windrows where they sink. The rolls need both the current's shear and the drift's; their growth rate sees only the product; and the split between the two decides which motion gets the energy.

8 figures