Viscosity

The surface that moves with the flow

A clean gas bubble feels two-thirds of the drag a rigid sphere of the same size would, and the formula has no density in it anywhere. What buys the third is that the bubble's surface is free to move — and real bubbles in ordinary water do not get it, for a reason that is a millionth of a per cent of the water by mass.

Worth reading first: The world with no inertia · How many things a flow must be told.

Every problem in this collection until now has put a rigid wall in the fluid. No slip, no flow through, and the surface goes where it is told. That is a boundary condition and it is a modelling choice, and the interesting question is what changes when it is dropped.

A liquid drop or a gas bubble has a surface the fluid outside can drag along — which is the same freedom a swimming sheet uses on purpose, arriving here as something the drop has no say in. The shear the outside applies is transmitted to the inside, which circulates; the interface moves; and the fluid outside is sheared less than it would be against a rigid wall of the same shape.

The consequence is a drag that is smaller, and the amount is exact.

Inside a drop, at a viscosity ratio of 0.5. Creeping flow past a spherical drop whose surface is free to move. The outer streamlines are fore-and-aft symmetric, as every creeping flow is; the inner ones are a closed circulation, driven by the shear the outside applies to the interface. The interface itself is moving at 33 per cent of the free stream, which is what a rigid sphere forbids and is the whole of why the drag is lower.
Fig. 1 Creeping flow past a drop whose surface can move. The outer streamlines are fore-and-aft symmetric, as every creeping flow is. The inner ones are a closed circulation driven by the shear the outside applies. The interface is moving at a third of the free-stream speed, which is what a rigid sphere forbids.

Why it is a boundary condition and not a force

Before the arithmetic, it is worth being precise about what has changed, because the usual summary — “a bubble has less drag because it is not solid” — hides the mechanism.

Nothing about the equations has changed. The fluid outside a bubble obeys the same Stokes equations it obeys outside a marble. What has changed is what is imposed where the fluid meets something else, and there are two conditions at a clean interface where a rigid wall has one.

At a rigid wall the velocity is given: it is the wall’s. At a clean interface the velocity is not given — it is whatever the two fluids agree on — and what is given instead is that the tangential stress must be continuous across it, since an interface with no thickness has no mass and cannot carry a net force. Two conditions replace one, and the extra unknown is the interface’s own motion.

That is why the answer depends on the viscosity ratio and on nothing else: the ratio is what decides how much motion a given stress buys on each side. And it is why the whole effect can be destroyed by a monolayer, since a monolayer does have a tangential stress of its own and restores the velocity-is-given condition without changing anything about either bulk fluid.

The drag

Hadamard and Rybczynski solved it independently in 1911, and the answer is

F=6πμaU23+κ1+κ,κ=μinμout.F = 6\pi\mu a U \cdot \frac{\tfrac23 + \kappa}{1 + \kappa}, \qquad \kappa = \frac{\mu_{\text{in}}}{\mu_{\text{out}}}.

At κ\kappa \to \infty the interior is infinitely viscous, the surface cannot move, and the factor is one — Stokes’ law, recovered exactly. At κ=0\kappa = 0 the interior offers no resistance at all and the factor is two-thirds.

The range is 2/3 to 1 and nothing lies outside it. There is no drop, of any material, in any liquid, that has less than two-thirds of a rigid sphere’s drag by this mechanism.

Two-thirds to one, and nothing outside it. The drag on a spherical drop as a fraction of Stokes' drag on a rigid sphere of the same size, against the ratio of the viscosities. It runs from exactly two-thirds for a clean gas bubble to exactly one for a rigid particle and takes no value outside that range, however light or heavy the drop is — the formula has no density in it anywhere. A drop of water in air and a drop of mercury in water sit at opposite ends of it.
Fig. 2 The drag factor across the whole range of viscosity ratios. A drop of water in air sits at the rigid end, because water is fifty times as viscous as air and its interior barely moves. An air bubble in water sits at the mobile end. The formula contains no density anywhere on this axis.

Where the heat goes, and the ten per cent

This is a phase about dissipation, and the drop offers a split that a rigid sphere cannot.

The total dissipation is F ⁣ ⁣UF\!\cdot\!U, exactly — the balance established for a translating sphere, now over two regions instead of one, and closing to a part in a hundred thousand on a shell integration out to two thousand radii. What is new is that it can be apportioned.

The share destroyed inside the drop works out at

DinDtotal=κ(1+κ)(2+3κ),\frac{D_{\text{in}}}{D_{\text{total}}} = \frac{\kappa}{(1+\kappa)(2+3\kappa)},

which is zero at both ends and has a maximum of 10.10 per cent at κ=2/3=0.8165\kappa = \sqrt{2/3} = 0.8165. At κ=1\kappa = 1 it is exactly a tenth.

How much of a drop's heat is made inside it. The share of the total dissipation destroyed within the drop, against the ratio of the two viscosities. It is zero at both ends and has a maximum of 10.1 per cent in the middle, at a ratio of √(2/3). A bubble makes none of its heat inside because there is nothing viscous in there to make it; a very viscous drop makes none because its interior barely moves. Nine tenths of the heat a falling drop makes is made in the fluid it is falling through, at every viscosity ratio there is.
Fig. 3 The share of the heat made inside the drop, against the viscosity ratio. A bubble makes none of it inside, because there is nothing viscous in there to make it. A very viscous drop makes none, because its interior barely moves. And in between the most it ever reaches is a tenth — so nine-tenths of the heat a falling drop makes is made in the fluid it is falling through, at every viscosity ratio there is.

That ceiling is not obvious and it is worth saying why it exists. The interior circulation is driven by the exterior; its velocity scale is set by what the interface is doing, which is at most half the free stream; and its volume is small compared with the region of the exterior that is being disturbed, because the exterior disturbance extends for many radii. Both factors work against the interior, and the product of them never exceeds a tenth.

The rise speed that clean water does not deliver

Put the drag law into a buoyancy balance and a clean air bubble in water should rise 49 per cent faster than Stokes’ law says — the factor is 1/(2/3+κ)/(1+κ)1/(2/3 + \kappa)/(1+\kappa) evaluated at κ=0.018\kappa = 0.018, and it is a constant, so the whole line shifts rather than tilting.

A bubble that should rise half as fast again. The terminal rise speed of a clean air bubble in water, computed both ways: with Stokes' rigid-sphere drag and with Hadamard and Rybczynski's mobile-interface drag. The mobile answer is 49 per cent faster at every size, because the factor is a constant. Small bubbles in ordinary water are measured rising at the rigid speed, not this one — the reason is surfactant, and it is what the picture cannot show.
Fig. 4 The two predictions for a rising bubble in water, over a decade and a half of size. The mobile-surface answer is a constant factor above the rigid one at every size. Measurements of small bubbles in ordinary tap water lie on the lower line.

They lie on the lower line. A bubble smaller than about a millimetre in ordinary water rises at Stokes’ speed, as though its surface were rigid — and the discrepancy was a genuine puzzle for decades before it was resolved.

What immobilises a surface

The resolution is surfactant, and the argument is worth following because it is a beautiful piece of reasoning about a boundary condition.

Any real water contains trace surface-active material — parts per billion is enough. It adsorbs onto the bubble’s surface. As the bubble rises, the surface flows from the front to the back, and it carries the adsorbed material with it, which accumulates at the rear.

A gradient of surface concentration is a gradient of surface tension, and a gradient of surface tension is a stress in the interface: it pulls from where the tension is high towards where it is low, which is from the back towards the front, which is against the motion of the surface.

So the surface stops. Not because anything is holding it, but because the contamination it swept backwards is pulling it forwards, and the balance is reached when the surface velocity is essentially zero. The bubble then presents a rigid boundary to the fluid outside, and its drag is Stokes’.

The bubble’s boundary condition is set by a contaminant present at parts per billion, and the factor it costs is a third of the drag. That is as sharp an example as this subject offers of a boundary condition mattering more than an equation.

The order of the equation is the number of conditions. What each model of a fluid allows to be said at a wall. Euler's equations are first order in the wall-normal direction and take one condition — the flow may be told not to go through the wall and may not be told anything about going along it, which is why an inviscid body has no friction and no drag. Navier–Stokes is second order there and takes two, and the second one is no slip. Viscosity does not make the same problem harder, it makes it a different problem, with one more thing that has to be true at every wall. The boundary-layer equations are the parabolic middle case: two conditions at the wall and a matching rather than a value at the outer edge.
Fig. 5 Where such a failure sits. The interior equations and the conditions on their boundaries are separate claims, and a model can be exactly right about the first while being wrong about the second — which is harder to notice, because it shows up as a wrong answer rather than as a wrong equation.

The energy version of the same argument

The surfactant story can be told in this essay’s own currency, and doing so makes it sharper.

A clean bubble and a contaminated one of the same size, rising at their own terminal speeds, are both in a steady state, so in both cases the buoyancy power equals the total dissipation. The clean bubble rises faster, so it is releasing potential energy faster — 49 per cent faster — and therefore dissipating 49 per cent more per second.

That is the opposite of what “less drag” suggests. A clean bubble is not a lower-loss object; it is a faster one, and it makes more heat per unit time precisely because it goes faster. What it is is more efficient per unit height climbed: the energy released per metre of rise is the buoyancy times a metre either way, so both bubbles dissipate exactly the same amount of energy getting to the surface.

The mobility does not save energy. It saves time. Which is the right way round for the applications: a degassing furnace and a flotation cell care about how quickly bubbles arrive, not about how much heat they make on the way.

How much of a drop's heat is made inside it. The share of the total dissipation destroyed within the drop, against the ratio of the two viscosities. It is zero at both ends and has a maximum of 10.1 per cent in the middle, at a ratio of √(2/3). A bubble makes none of its heat inside because there is nothing viscous in there to make it; a very viscous drop makes none because its interior barely moves. Nine tenths of the heat a falling drop makes is made in the fluid it is falling through, at every viscosity ratio there is.
Fig. 6 And where that heat is made in the two cases. A contaminated bubble makes all of its heat outside itself, because its interior is not moving. A clean one makes 0.86 per cent inside — which is small, and is the entire difference between a boundary condition that is a wall and one that is not.

The surface that is rigid at one end

The account above has two states — clean and contaminated — and a real bubble is usually in neither. The reason is in the mechanism itself: the surfactant is swept backwards, so it does not coat the bubble evenly. It piles up at the rear and leaves the front bare.

What that produces is a stagnant cap: a region around the rear stagnation point where the surfactant is dense enough for the Marangoni stress to hold the interface still, and a leading region where there is not enough to matter and the surface moves freely. The boundary between them sits at some angle from the rear, and that angle — not the presence or absence of surfactant — is what the drag depends on.

The problem is harder than either limit, because the boundary condition changes type partway along the surface: velocity given over the cap, stress given over the rest, with the dividing angle an unknown fixed by the surfactant’s own balance of convection against adsorption. It has an exact solution anyway, obtained by Sadhal and Johnson in 1983, and the drag it gives rises monotonically from the Hadamard–Rybczynski value at zero cap angle to Stokes’ law at a cap covering the whole sphere. Both limits of this essay are the two ends of one curve, and the intermediate cases are not a smudge between them but states a bubble genuinely occupies.

Three observations become explicable at once.

Rise speed varies continuously with water quality rather than switching between two values. Adding surfactant to a clean system slows the bubbles gradually, which is what the cap angle doing so predicts and what two discrete states cannot.

A bubble slows down as it rises. Released clean into water that is not, it starts near the mobile speed and decelerates over the first part of its travel while the cap grows to its equilibrium extent — so a measurement in a short column and one in a tall column disagree, and the disagreement is not experimental error.

And size decides which limit wins. A bubble’s surface area grows as the square of its radius while the surfactant available near it grows as the cube of the distance swept, so a large bubble presents more surface than the trace impurity can cover and keeps a bare front. That is the quantitative version of the paragraph below about large bubbles staying mobile, and it says the threshold is a comparison of two rates rather than a size anybody could look up.

So the honest statement of the boundary condition is not that it is one thing or the other. It is that the interface carries a field of its own, and where that field is strong enough the fluid outside sees a wall.

When the surface is mobile after all

Three cases, and they are the cases where the effect is used.

Large bubbles. A bubble big enough has a surface area growing faster than the surfactant can adsorb onto it, and the leading portion stays clean. Bubbles above a few millimetres in water do show mobile-interface behaviour — along with a great deal else, since they are no longer spherical and no longer creeping.

Very clean systems. Distilled and degassed water, carefully handled, gives bubbles that rise at the Hadamard–Rybczynski speed. It is difficult to arrange and it is the experiment that settled the question.

Liquid metals. Molten metals have no surfactant chemistry of the kind water has, and drops in them behave as the theory says. Bubbles in liquid steel — which matters, since that is how steel is degassed — are genuinely mobile.

The same story in a suspension

The rigidity assumption dropped here is the one Einstein’s coefficient is built on, and dropping it changes that number by exactly the same mechanism.

Taylor gave the emulsion version in 1932: the coefficient is not 5/2 but

52κ+25κ+1,\frac{5}{2}\cdot\frac{\kappa + \tfrac25}{\kappa + 1},

running from 1 for clean bubbles to 5/2 for rigid particles. A dilute foam is thickened by only ϕ\phi rather than by 2.5ϕ2.5\phi — and the same surfactant argument applies, so a real foam behaves like a suspension of rigid particles and is thickened by the full 5/2.

That is the same result twice, in two different measurements, with the same trace impurity deciding both.

Where Einstein's line stops being the measurement. The viscosity of a suspension of rigid spheres relative to the liquid's, against the volume fraction. Einstein's 1 + 5φ/2 is exact for one sphere and holds while the spheres cannot feel one another, which is up to about five per cent by volume. Batchelor and Green's two-sphere term takes it a little further; beyond about a fifth nothing derived works and the curve drawn is a fit, which diverges at a maximum packing that is itself a measurement.
Fig. 7 The rigid-particle version, for comparison. Every point on that curve is for particles that cannot deform and whose surfaces cannot move. An emulsion of clean drops sits below it by a factor that depends only on the viscosity ratio, and an emulsion of dirty ones sits on it.

What “mobile” means quantitatively

The interface speed is worth a number, because the phrase “the surface moves” invites a picture of the surface travelling with the flow, and it does not.

At the drop’s equator, the surface moves at

U2(1+κ)\frac{U}{2(1+\kappa)}

as a fraction of the free stream. For a clean bubble that is very nearly half; for equal viscosities a quarter; for a rigid particle zero. The most a free surface ever moves is half the free-stream speed, and it does so only in the limit of an inviscid interior.

That ceiling explains the drag ceiling. The exterior fluid is sheared between the free stream and the surface, and halving the velocity difference at the wall does not halve the shear — the disturbance adjusts over the same length scale — so the drag falls by a third rather than by a half.

It also explains why the effect is so much smaller than people expect it to be. A bubble is, from the point of view of the fluid outside it, a slightly slippery sphere rather than a hole. The whole mechanism is worth 33 per cent, and the difference between a bubble and a lead shot of the same size falling in the same oil is entirely in the buoyancy, not in the drag.

Inside a drop, at a viscosity ratio of 0.02. Creeping flow past a spherical drop whose surface is free to move. The outer streamlines are fore-and-aft symmetric, as every creeping flow is; the inner ones are a closed circulation, driven by the shear the outside applies to the interface. The interface itself is moving at 49 per cent of the free stream, which is what a rigid sphere forbids and is the whole of why the drag is lower.
Fig. 8 The clean-bubble limit drawn out. The interior circulation is at its most vigorous relative to the exterior, the interface is moving at close to half the free stream, and the total dissipation is two-thirds of what a rigid sphere of the same size would make at the same speed.

What the picture cannot show

The drop is spherical and stays spherical. That requires surface tension to dominate both viscosity and inertia — a small Weber number and a small capillary number. A drop moving fast enough deforms, and what it deforms into is not a tear.

The interface is clean, and almost none are. Everything above is what happens when the surface carries nothing. The surfactant argument is described here and not computed; doing it properly needs an adsorption isotherm and a surface transport equation, which is a different subject.

There is no Marangoni flow of the other kind. A temperature gradient along an interface also produces a surface-tension gradient and drives a flow, and a bubble in a temperature gradient migrates towards the warm side with no gravity involved at all. None of that is here.

The dissipation is the incompressible one. Both fluids here are treated as having no divergence anywhere, so the second viscosity contributes nothing — exactly rather than approximately, which is worth saying for a bubble, whose interior is a gas.

And the flow is creeping. The formula holds at Reynolds numbers well below one. The bubbles it is usually applied to are not, which is why the terminal-velocity figure marks the points that have left the regime.

Two places the third of the drag is worth money

Flotation. Mineral processing separates valuable particles from waste by attaching them to rising bubbles, and the rate of the whole plant is the rate at which bubbles rise through the cell. A thirty-three per cent difference in rise speed is a thirty-three per cent difference in throughput, and it is decided by the water chemistry rather than by anything mechanical. Flotation cells are operated with deliberate surfactant, because a bubble that is too mobile does not hold its particle.

Steel degassing. Argon is bubbled through molten steel to strip dissolved gases, and the contact time is set by the rise speed. Liquid metals have no adsorbed organic layer, so the bubbles are genuinely mobile and the design uses the two-thirds factor rather than Stokes’ law — one of the few industrial processes where the clean-interface answer is the right one.

Between them those make the point that a third of a drag is not a rounding error in any process whose economics is a residence time.

Who found it, and when

Hadamard and Rybczynski published within months of one another in 1911, unaware of each other’s work, and the result carries both names. The discrepancy with measured bubble speeds was noticed almost immediately and was not settled until Frumkin and Levich’s surfactant argument in the 1940s, with the experimental confirmation — clean bubbles rising at the predicted speed — coming later still.

The surprising connection is with a completely different kind of contamination problem. The mechanism here — a swept-along impurity accumulating at a stagnation region and generating a stress that opposes the sweeping — is the same mechanism by which a surfactant stabilises a foam, by which a tear film keeps an eye wet, and by which the Marangoni effect drives a droplet of alcohol across a water surface. In every case a vanishing amount of material at an interface governs the behaviour of a great deal of fluid, because an interface has no thickness and therefore no dilution.

Where the ladder goes next

Beside this rung is Einstein’s coefficient, which is the same boundary condition in a bulk property, and a swimmer that cannot go backwards, which is the other essay in this collection about a surface that moves on purpose.

Below it is how many things a flow must be told, which is where the arithmetic of boundary conditions is set out, and the world with no inertia, which is the regime this all holds in.

What links here

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

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 conditionBubbleCreeping flowDissipationDragDropMobile interfaceSurfactantTerminal velocityViscosity ratio