Fluids at work

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.

Worth reading first: When the flow lets go.

Drag goes as the square of speed. Everybody knows this, it is right nearly everywhere, and it is built into every intuition about how things move through air.

There is one place where it is not merely inaccurate but backwards. Over a narrow band of speeds, a smooth sphere pushed harder experiences less force, not more. The drag falls as the speed rises, by nearly thirty per cent, and the reason is that the boundary layer on the ball changes character partway through.

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.
Fig. 1 The drag coefficient of a sphere against Reynolds number, on log axes. The smooth curve is Morrison’s correlation — a fit to measurements, 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 six, which is a stated model of what roughness does and not a measurement of any real ball.

What is borrowed here, stated before anything is built on it

This essay is the first in the field to rest on a correlation for its main quantity, and the terms should be clear before any argument is made on top of it.

This site computes no turbulent boundary layer at any Reynolds number. The reason is arithmetic rather than effort: resolving one needs a grid whose cost nobody can pay. So the separation point on a ball — which is the thing that decides everything below — cannot be predicted here, and any figure claiming otherwise would be inventing.

What is used instead is Morrison’s correlation for the drag of a smooth sphere: an algebraic fit, valid from Stokes’ law at the bottom to Re=106\mathrm{Re} = 10^6 at the top, which reproduces the plateau in the middle and the crisis near Re=2.6×105\mathrm{Re} = 2.6\times10^5. It is a fit. No part of it is derived and no assertion in this essay tests it.

What is computed, and what the assertions are about, is what follows given that curve: the band in which the drag force falls, located by bisection on the logarithmic slope; the trajectories, whose energy budgets close to a part in a million; and the comparison between the two balls.

The roughened curve is a second, weaker kind of statement. It is the smooth correlation shifted along the Reynolds axis by a stated factor, which is the simplest thing that captures what roughness does — trip the layer early, bring the crisis forward. It is a model of an effect and not a measurement of a golf ball, and it gets one thing wrong that the essay’s last section returns to.

The crisis, and what changes at it

Below the crisis, the boundary layer on the front of a sphere is laminar. It climbs the pressure rise on the back half, and a laminar layer climbing a pressure rise gives up early — at about eighty degrees from the front stagnation point, barely past the equator. The wake behind is as wide as the ball, the base pressure in it is low, and almost all the drag is pressure drag.

Above the crisis, the layer goes turbulent before it gets there. A turbulent layer carries much more momentum near the wall, because turbulence mixes fast fluid down into the slow region, so it can climb a much steeper pressure rise before separating — around a hundred and twenty degrees, well round the back. The wake is narrower, the base pressure higher, and the pressure drag much smaller.

The trade is the essay’s refutation in one line:

  • Friction goes up. A turbulent boundary layer has several times the skin friction of a laminar one at the same Reynolds number, which is computed elsewhere on this site for a flat plate and rises further with Reynolds number.
  • Pressure drag goes down, by more. On a bluff body, pressure drag is the dominant term, and narrowing the wake cuts it dramatically.

Total drag coefficient falls from about 0.5 to about 0.15. The smaller of the two contributions went up; the larger one collapsed.

It is worth pausing on why the pressure term is so large on a bluff body in the first place, because the answer is the site’s founding paradox. Ideal flow predicts no drag at all: the pressure recovers perfectly on the back of a sphere, the fore-and-aft distribution is symmetric, and the integral is exactly zero. Every bit of the pressure drag a real sphere has is the failure of that recovery, and the failure happens because the flow separates. So the pressure drag is not a separate mechanism from separation — it is separation, measured as a force, and anything that delays separation attacks it directly.

That also explains the size of the effect. Moving the separation point from eighty degrees to a hundred and twenty does not shave a percentage off the drag; it recovers a large part of the pressure on the back of the ball that ideal flow said should have been there all along.

The force, not the coefficient

A falling coefficient is unremarkable — almost every drag coefficient falls with Reynolds number somewhere. What is remarkable is that here the force falls, and that requires something much stronger.

Drag force is 12ρU2ACD\tfrac{1}{2}\rho U^2 A\,C_D, and at fixed size and fluid UReU \propto \mathrm{Re}, so the force goes as Re2CD\mathrm{Re}^2 C_D. For that to fall, CDC_D must fall faster than Re2\mathrm{Re}^{-2}: the logarithmic slope dlnCD/dlnRe\mathrm{d}\ln C_D/\mathrm{d}\ln\mathrm{Re} must be steeper than 2-2.

The one stretch where speeding up costs less. The drag force on a sphere, as Re²C_D on log axes, for a smooth ball and one whose boundary layer is tripped early. Drag normally rises with the square of the speed and this curve normally rises with slope 2. Through the crisis it falls: between Re = 2.40e+5 and 3.53e+5 the drag force on a smooth sphere is 29.2% lower at the top of the band than at the bottom, and the ball is going faster.
Fig. 2 The drag force itself, as Re²C_D on log axes. Everywhere else the curve rises with slope close to two, which is the ordinary quadratic law. Through the crisis it falls, and the ball is going faster while being pushed less hard.
Steeper than minus two, and only just. The logarithmic slope of the drag force, d ln(Re²C_D)/d ln Re, through the drag crisis. A drag force that grows with the square of speed sits at 2; a coefficient that merely falls pulls it down a little. Only where the slope crosses zero does the force itself fall, which needs the coefficient to fall faster than the square. It happens over about a sixth of a decade and nowhere else on the whole curve.
Fig. 3 The logarithmic slope of the force through the crisis, with the region below zero marked. The band is found by bisection on the slope rather than read off the curve, and it runs from Re = 2.40·10⁵ to 3.53·10⁵ — about a sixth of a decade, and nowhere else on the whole range.

The band comes out at Re=2.40×105\mathrm{Re} = 2.40\times10^5 to 3.53×1053.53\times10^5, with the force falling by 29.2 per cent across it. In air, for a 70 mm ball, that is roughly 51 to 76 metres per second: a sphere at 76 m/s is being pushed backwards less hard than the same sphere at 51.

Both halves of that are asserted. A curve reported as having a crisis must actually have a band where the force falls; and — the more useful half — a curve with no such band must not be reported as having one. That second refusal exists because the failure is symmetric: a correlation with the crisis term dropped still falls with Reynolds number everywhere, and would pass any test that merely asked whether CDC_D decreases.

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 12 — a stated model of what roughness does, which is to trip the boundary layer early, and not a measurement of any real ball.
Fig. 4 A more heavily roughened ball, drawn by shifting the correlation twelve-fold rather than six. The crisis has moved to Re ≈ 2·10⁴ and the post-crisis plateau is reached at a walking pace — which is why a badly scuffed ball behaves oddly at speeds where a new one is entirely conventional.

The flight, integrated

A ball in flight is a Reynolds number sweep. It launches fast, slows, and its drag coefficient walks along the curve while it does so — which is why a trajectory, rather than a coefficient, is the figure that settles the argument.

The same strike, from a smooth ball and a rough one. Two trajectories integrated at 70 m s⁻¹ and 12 degrees with 300 radians per second of backspin, differing only in whether the boundary layer is tripped. The rough ball spends its flight past the drag crisis and the smooth one does not, and the ranges come out at 192 and 411 metres. Both integrations close their energy budget to better than a part in a million, so the difference is the drag law and not the integrator.
Fig. 5 Two trajectories at 70 m/s and twelve degrees with 300 rad/s of backspin, differing only in whether the boundary layer is tripped. The tripped ball spends its flight past the crisis and the smooth one does not; the ranges come out at 411 and 192 metres. Both integrations close their energy budget to better than a part in a million.

Two assertions run on every trajectory, and both are about the integrator rather than the physics — which is the right place for them, because the comparison between two balls is only worth anything if the arithmetic is not leaking.

The Magnus force does no work. A force perpendicular to the velocity cannot, and the check evaluates the lift acceleration through the same code path that moves the ball and dots it with the velocity. It comes out at 101710^{-17} joules. The failure it catches is a sign or an index the wrong way round in the two components, which produces a perfectly plausible trajectory — and a sign error of exactly that kind has been caught by measurement rather than by looking four times in other fields on this site.

The energy budget closes. The kinetic and potential energy lost must equal the work done against drag, integrated along the path. It comes out at 3×1093\times10^{-9} of the launch energy. An integrator leaking a per cent would change the range by metres, which is the size of the effect being argued about.

The same strike, from a smooth ball and a rough one. Two trajectories integrated at 45 m s⁻¹ and 20 degrees with 150 radians per second of backspin, differing only in whether the boundary layer is tripped. The rough ball spends its flight past the drag crisis and the smooth one does not, and the ranges come out at 106 and 186 metres. Both integrations close their energy budget to better than a part in a million, so the difference is the drag law and not the integrator.
Fig. 6 The same comparison at a slower, steeper launch with less spin — a struck ball rather than a driven one. The advantage is smaller, because at 45 m/s a smooth ball is further below the crisis and a tripped one spends part of its flight below it too. Roughness is worth most where the flight sits across the transition, which is exactly where a golf ball’s does.

The second figure is the more instructive one. Roughness is not a free improvement; it is a transposition. It moves the crisis to a lower speed, which is worth a great deal if the ball flies near the shifted crisis and worth nothing if it flies far below it. A rough ball thrown slowly is simply a rough ball, with all the extra friction and none of the benefit.

Where the drag actually is

The split between the two contributions is worth seeing directly, because the folk explanation is entirely about the smaller one.

A wing's two drags, and the lift at which they are equal. Friction drag and induced drag plotted against lift coefficient, with their sum above them. Friction is flat, because a surface costs the same whatever the wing is doing; induced drag rises as the square of the lift. The total is least where the two are equal, and that is also the point of best glide.
Fig. 7 The two drags a body pays, in the setting where this site computes them exactly: friction, from Blasius’ solution, and the drag that comes from pressure. On a streamlined body at small incidence friction dominates and smoothness wins. On a bluff body pressure dominates, and that reverses every conclusion about surface finish.

That single distinction — streamlined or bluff — decides which way the trade goes, and it is why the dimple argument does not generalise. An aircraft wing is streamlined; its drag is mostly friction; a tripped boundary layer is a straightforward loss, and enormous effort goes into keeping laminar flow on one. A golf ball is bluff; its drag is mostly pressure; tripping the layer is worth a factor of three. The same modification, on two shapes, in opposite directions.

The one stretch where speeding up costs less. The drag force on a sphere, as Re²C_D on log axes, for a smooth ball and one whose boundary layer is tripped early. Drag normally rises with the square of the speed and this curve normally rises with slope 2. Through the crisis it falls: between Re = 2.40e+5 and 3.53e+5 the drag force on a smooth sphere is 29.2% lower at the top of the band than at the bottom, and the ball is going faster.
Fig. 8 The drag force for that same twelve-fold shift, against the smooth sphere’s. The falling band has moved with the crisis, and everywhere outside it both curves rise with the ordinary slope of two.

The same crisis, on everything bluff

A sphere is not special. Every bluff body has a drag crisis, at a Reynolds number set by its own geometry, and the mechanism is identical each time: the boundary layer goes turbulent before it has to climb the pressure rise, separation moves aft, the wake narrows.

A circular cylinder has one at about the same Reynolds number as a sphere, and it is the reason the site’s regime ladder puts a threshold at 2·10⁵ — that is where a cylinder’s drag coefficient falls from about 1.2 to about 0.3. The cylinder case is the one this site can draw the flow of at low Reynolds numbers, and it is a very different picture at each of them.

The consequences turn up wherever a bluff body moves fast enough:

  • Cycling. A rider’s limbs and body are bluff, and much of the speed range of a racing cyclist sits near the transition for a cylinder of a limb’s diameter. Textured skinsuit fabric on the arms is a trip wire, and it is legal in some codes and banned in others.
  • Ski jumping and speed skating. Same argument, same fabrics, same regulatory arguments.
  • Chimneys, cables and bridge decks. Here the crisis is a hazard rather than a benefit: a structure that sits near its own critical Reynolds number has a drag that changes sharply with wind speed, and one side of a circular member can be in a different regime from the other. Helical strakes on a chimney exist to force the issue.

None of these is a different physics. They are the same curve read at a different Reynolds number, which is the site’s oldest claim — the shape is not the question, the ratio is — applied somewhere with money on it.

The crisis is not a property of the ball

Everything above treats the critical Reynolds number as a number belonging to a sphere, the way its diameter does. It is not. It belongs to the sphere and the stream it is in, and the second half of that pairing moves it further than roughness does.

The reason follows from the mechanism. The crisis happens when the boundary layer goes turbulent before it reaches the separation point, and transition is triggered by disturbances — which may come from the surface, as a dimple or a trip wire, or may arrive already present in the oncoming air. A stream that is itself turbulent supplies them for free, so the layer transitions earlier, and the whole curve shifts left exactly as a roughened ball’s does. From the drag curve alone the two causes are indistinguishable.

The size of it is large. A sphere in genuinely quiet air — free flight, or a tunnel built for low turbulence — reaches the post-crisis state at about 3.9×1053.9\times10^5. The same sphere in a stream carrying a per cent or two of turbulence reaches it below 2×1052\times10^5, and in a deliberately stirred one lower still. That is a factor approaching two, which is a third of the six-fold shift this essay uses to represent a golf ball’s dimples, produced by nothing being done to the ball at all.

So it was turned into an instrument. Through the 1920s and 1930s a standard sphere was the accepted way of measuring how quiet a wind tunnel was: run it up, find the Reynolds number at which its drag coefficient passed a fixed value, and divide the free-air figure by the one the tunnel gave. The ratio was the tunnel’s turbulence factor, and every drag measurement made in that tunnel was corrected by it. Tunnels of the period ran from about 1.1 for the best to nearly 3 for the worst, and a model tested in one was effectively being tested at a higher Reynolds number than its speed implied.

That is a slightly vertiginous piece of history, and it has a direct bearing on this essay’s borrowing. The sphere’s own drag curve was for two decades the calibration standard for the instrument used to measure drag curves — so the correlation the whole argument rests on is a fit to measurements each of which carries its own stream’s disturbance level inside it. The band located above by bisection, 2.402.40 to 3.53×1053.53\times10^5, is therefore a property of a dataset rather than of spheres, and quoting it to three figures says more about the fit than about the physics.

None of that touches the argument. The falling force is real, the mechanism is real, and the trade between friction and pressure runs the same way whatever triggers the transition. What moves is the speed at which a particular ball, on a particular day, meets it.

Where the model stops

The shift is a model, and it is wrong in a knowable way. Shifting the smooth correlation predicts that a tripped ball reaches the smooth ball’s post-crisis coefficient of about 0.13. A real golf ball sits nearer 0.25. The difference is that dimples do not only trip the layer — they are themselves roughness elements with their own form drag, and the fully-rough asymptote of a sphere is higher than the smooth-sphere post-crisis value. The shift model captures the transition and not the penalty, and the trajectory figures therefore overstate the tripped ball’s range.

That is stated rather than corrected because correcting it would mean fitting a second parameter to golf-ball data, which would make the figure a curve fit rather than an argument.

No spin effect on drag. Spin changes the separation pattern and therefore the drag, not only the lift. The model here has spin entering the lift alone.

The lift model is linear and stated. CL=kωR/VC_L = k\,\omega R/V with kk chosen so that a driven ball comes out near the measured lift coefficient. Real lift coefficients saturate at high spin. This matters for the apex of the trajectory and not for the drag argument.

And nothing here is a solution. The separation angles, the base pressure, the correlation and the shift are all borrowed or stated. What this site contributes is the arithmetic around them, and the distinction is worth keeping sharp: the crisis is a measurement, and the falling force is a calculation.

Who found it, and when

Gustave Eiffel measured the drag crisis in 1912, dropping spheres down the tower that carries his name and later in the wind tunnel he built at Auteuil. Prandtl explained it almost immediately, and in 1914 demonstrated it in the neatest way anybody has: he fitted a trip wire — a thin ring of wire around the front of a sphere — and photographed the wake narrowing. A sphere with a wire on it had measurably less drag than the same sphere without.

That experiment is the whole essay in one object. Something was added to the surface, the friction went up, and the total drag went down, because the added thing changed where the flow let go.

The dimple came from the other direction entirely. Golfers in the 1840s noticed that old scuffed gutta-percha balls flew further than new smooth ones, which is an observation seventy years ahead of any explanation, and the deliberately patterned ball dates from the 1890s. It is one of the clearest cases in engineering of a practice that worked for half a century before anybody could say why, and of an explanation that, once available, turned out to be the opposite of the obvious one.

Where the ladder goes

This rung is about a symmetric ball, where moving the separation point changes the drag. The next is about an asymmetric one, where moving it on one side only changes the direction.

A cricket ball swings without spinning at all. The seam trips the boundary layer on one side, that side separates later, the pressure distribution is no longer symmetric, and the ball goes sideways. The calculation of the side force is the same pressure integral this essay’s mechanism implies — performed both by quadrature and in closed form — and it has a property that makes it unusually instructive: the sign of the answer depends on one borrowed number, and the site’s own inviscid solver supplies the value that gets it wrong.

What links here

Computed from the collection rather than written here: the essays that point at this one.

Reads more easily once this is understood

Essays that name this one as worth reading first.

Shares its objects with

Essays naming at least two of the same things, that neither author linked.

Named objects

A dashed tag is an object no other essay names yet.

Boundary layerCorrelationDrag crisisPressure dragReynolds numberRoughnessSeparationSkin frictionTransitionWake