Viscosity

Everything happens in a layer you cannot see

Air has so little viscosity that ignoring it works almost everywhere. Almost everywhere leaves out a film next to the surface, perhaps a millimetre thick, and that film decides drag, stall and whether an aircraft flies at all.

Air is about fifty times less viscous than water and water is not very viscous either. Over almost any flow anybody cares about, viscosity contributes so little that leaving it out entirely is accurate to a few percent.

Leaving it out entirely also predicts that nothing has any drag, which is wrong by all of it.

Both statements are true, and the reconciliation is the single most important idea in the subject: viscosity is negligible everywhere except in a very thin layer next to a solid surface, and that layer decides the quantities the rest of the flow cannot.

Flow past a cylinder at Re 40A real fluid past a circular cylinder. At low Reynolds number the flow closes up behind the body much as the ideal theory says; as it rises the flow separates and a region of reversed flow appears behind, which is where drag comes from.separatedrecirculation 0.56 Dviscous flow, solved on a coarse grid — the bubble is under-resolvedRe = 40
Fig. 1 A real flow past a cylinder. Outside a thin region next to the surface it looks much like the inviscid solution; inside that region it is doing something the inviscid solution cannot represent.

Air sticks

Start with the fact that makes the layer exist. At a solid surface, the fluid does not slide. It sticks: the velocity of the air in contact with a wing is exactly the velocity of the wing.

That is the no-slip condition, and it is not obvious. It was argued about through most of the nineteenth century, and it is not a consequence of anything more basic — it is an experimental fact about the interaction of a fluid with a solid, and it holds for every ordinary fluid and surface.

The consequence is immediate. A metre away from the wing the air is moving at flight speed. At the surface it is moving at zero. Between those two there must be a region where the velocity climbs from one to the other, and in that region there are large velocity gradients — which is exactly what viscosity responds to.

How thin

Thin enough to be startling.

The layer’s thickness grows along a surface roughly as the square root of distance, and scales inversely with the square root of the Reynolds number. For an airliner wing in cruise the boundary layer at the trailing edge is a few centimetres; over most of the wing it is millimetres. For a model aircraft it is thinner still in absolute terms.

So the region that decides the drag of a three-hundred-tonne aircraft is, over most of the wing, about the thickness of a coin.

That is why the ideal theory works as well as it does. Outside the layer, the flow genuinely is very nearly inviscid, and the exact solution is a good description of it. The theory is not wrong about the bulk of the flow at all. It is wrong about a film.

Why a film matters so much

Two reasons, and the second is the important one.

Friction. The velocity gradient at the wall is steep, and a steep gradient in a viscous fluid means a shear stress. Integrated over the wetted area, that is skin friction drag — for a streamlined body, most of the total.

Separation. This is the one that changes everything. Fluid inside the boundary layer has already lost energy to friction, and when it is asked to climb a rising pressure — as it always is over the rear half of any body — some of it cannot. It stops, reverses, and the whole outer flow lifts away from the surface.

Once the flow separates, the pressure never recovers over the rear of the body, the fore-and-aft cancellation that produced d’Alembert’s zero drag is destroyed, and what is left is pressure drag — which for a bluff body dwarfs the friction.

What the ideal theory predicts, and what happensThe same cylinder, the same free stream. On the left the exact inviscid solution, closing up behind the body and exerting no drag at all. On the right the real flow at the same conditions, separated, with a wake and therefore with drag.ideal flow — closes up, no dragreal flow at Re 100 — separatedleft: exact closed form · right: solved on a gridRe = 100
Fig. 2 The consequence, side by side. The ideal solution closes up behind the body; the real one does not, and the difference is a wake and a drag.

Two regions, two problems

Prandtl’s move was structural rather than computational, and it is worth stating in the form he did.

Outside the layer, solve the inviscid equations. They are exact, they are tractable, and they describe the bulk of the flow to a few percent. This is where the pressure distribution comes from.

Inside the layer, solve a simplified viscous equation. Because the layer is thin, pressure does not vary across it — the outer solution’s pressure is simply imposed — and the equations reduce to something one-dimensional in the cross-stream direction.

Match them at the layer’s edge, where the inner solution’s velocity has to approach the outer one’s.

That decomposition turns one intractable problem into two tractable ones. It is the original and still the clearest example of matched asymptotic expansions, a technique that has since spread across applied mathematics, and it was invented to explain why a wing has drag.

Ideal flow past a cylinderA uniform stream past a circular cylinder in a fluid with no viscosity. The solution is exact and closed-form: streamlines part at a stagnation point, run round the surface and close up perfectly behind, and the pressure recovers to exactly what it was in front.ideal flow — inviscid, irrotational, steadyno circulation
Fig. 3 The outer problem: the pressure distribution from the exact inviscid solution. The boundary layer takes this as given, and its job is to survive the climb up the rear half of it.

The most surprising part of Prandtl’s idea is that the thin layer is not merely a correction confined to itself.

By separating, the boundary layer changes the shape the outer flow sees. A cylinder with a separated wake presents the outer flow not with a circle but with something closer to a semi-infinite body, and the outer solution adjusts accordingly. The pressure distribution over the entire front of the body is different because of what happened in a film at the back.

So the coupling runs both ways: the outer flow sets the pressure gradient the layer must climb, and the layer, by giving up, resets the outer flow. That two-way coupling is why boundary layer theory is subtle rather than a simple additive correction, and why it took until 1904.

The pressure gradient is the whole story

Whether a boundary layer survives depends on one thing: which way the pressure is going.

Over the front of a body the flow is accelerating and the pressure is falling. A favourable gradient like that pushes the boundary layer along, thins it, and makes it stable. Nothing goes wrong here.

Over the rear the flow is decelerating and the pressure is rising. That is an adverse gradient, and it acts against the fluid nearest the wall — which is precisely the fluid with the least energy left. The velocity profile near the surface flattens, then develops an inflection, then reverses.

The moment the velocity gradient at the wall reaches zero, the flow has separated. Everything downstream of that point is a different flow.

So the design question for any body is not “how much friction does it have” but “how gently can the pressure be brought back up”. That is what a streamlined tail is for, and it is why the tail matters more than the nose.

What it costs to be blunt

The numbers make the point better than the argument does.

A circular cylinder in a typical flow has a drag coefficient near 1.2. A streamlined strut of the same frontal area — same thickness, but with a long tapering tail — has a drag coefficient near 0.06. Twenty times less drag for the same width.

The strut has considerably more surface area and therefore more skin friction. It wins because it has almost no pressure drag: the flow stays attached nearly to the tail, the pressure recovers, and the fore-and-aft cancellation that ideal flow predicts is very nearly restored.

Which is a useful way to see what streamlining is doing. It is not making the flow slippery. It is arranging for the inviscid theory’s beautiful, wrong answer to become nearly right.

The real flows here come from the vorticity–streamfunction solver in lib/flow.js, stepping the incompressible equations on a grid with no-slip imposed at the body.

What it resolves, and what the essays may therefore claim: the transition from attached flow to separated flow. At Reynolds number 1 and 10 there is no reverse flow anywhere behind the cylinder. At 40 and 100 there is, and the recirculating region grows with Reynolds number. Those are measured from the solved field by looking for negative streamwise velocity along the centreline, not asserted.

What it does not resolve, and what nothing here claims: the detailed structure inside the layer, the transition from laminar to turbulent, or turbulence itself. The grid is coarse enough that the computed recirculation length is shorter than the measured value — about 0.85 diameters at Reynolds number 40 against a laboratory figure near 2.2. The qualitative behaviour is right and the number is under-resolved, and a figure quoting a bubble length says so.

Getting even this far required fixing a solver that detonated: the wall vorticity was being computed from a streamfunction that had not converged, and the field grew by an order of magnitude every few steps. It looked like a timestep problem and was not.

How thick, in numbers

The scaling is worth having explicitly, because it is what makes “thin” quantitative.

For a laminar layer on a flat plate, the thickness grows as

δxRex\delta \sim \frac{x}{\sqrt{Re_x}}

where xx is distance from the leading edge and RexRe_x is the Reynolds number built on that distance. So the layer grows like the square root of distance, and shrinks like the square root of speed.

Put numbers in it. An airliner wing at cruise has a chord Reynolds number around 2×1072 \times 10^7. The laminar estimate gives a layer a fraction of a millimetre thick over the front of the wing; in practice the layer goes turbulent early and reaches a few centimetres by the trailing edge.

For a model aircraft with a Reynolds number of 10510^5 the layer is proportionally much thicker — which is one reason model aerodynamics is not simply scaled-down full-size aerodynamics, and why a trailing edge sharp to the eye may not be sharp to the flow, and why sections that work at model scale often behave differently at full size.

The general lesson is that “thin” is a statement about the Reynolds number and not about millimetres. At Reynolds number 10 the layer is comparable to the body, and the whole two-region picture stops being useful.

Skin friction, and where it goes

The other product of the layer is friction, and it is worth separating from separation.

The shear stress at the wall is the viscosity times the velocity gradient there:

τw=μuyy=0\tau_w = \mu \left. \frac{\partial u}{\partial y} \right|_{y=0}

Integrate that over the wetted surface and the result is skin friction drag. For a streamlined body with attached flow it is nearly the whole drag; for a bluff body with a big wake it is a small fraction.

A turbulent layer has a much steeper wall gradient than a laminar one — the mixing brings fast fluid close to the surface — so it has considerably more skin friction. That is the cost of the separation-resistance that makes turbulent layers useful, and it is why laminar-flow aerofoils, which try to keep the layer laminar over as much of the chord as possible, are worth the considerable trouble they cause.

Flow past a cylinder at Re 10A real fluid past a circular cylinder. At low Reynolds number the flow closes up behind the body much as the ideal theory says; as it rises the flow separates and a region of reversed flow appears behind, which is where drag comes from.attachedviscous flow, solved on a coarse grid — the bubble is under-resolvedRe = 10
Fig. 4 Reynolds number 10, where the layer is not thin at all. The two-region picture that makes boundary layer theory work has not yet come into existence at this scale.

Why it was not obvious

The century and a half between d’Alembert and Prandtl is worth explaining, because with hindsight the idea seems small.

The difficulty is that the layer is a singular perturbation. Setting viscosity to zero does not change the answer slightly; it changes the order of the governing equation, and therefore how many boundary conditions it can satisfy. The inviscid equation can enforce no flow through the wall but not zero flow along it. So the no-slip condition is precisely the one that gets dropped, and the whole of the physics that depends on it disappears with it.

That is why a small parameter produced a qualitative failure rather than a small error, and why nobody found the resolution by refining the existing theory. It needed somebody to say: solve two different problems in two different regions and match them at the boundary.

Prandtl’s eight pages

Ludwig Prandtl presented the idea at a mathematics congress in Heidelberg in 1904, in a paper of about eight pages. It was not immediately recognised — the audience was small, and the paper did not appear in a fluid dynamics journal because there was not really one to appear in.

Its consequences are hard to overstate. It made the drag of a body calculable, made stall comprehensible, and reunited hydrodynamics with the practical engineering that had given up on it. Prandtl’s students — Blasius, von Kármán, Munk, Betz, Schlichting — went on to build most of twentieth-century aerodynamics on top of it.

The paper is also a model of how to rescue a theory: not by adding terms until it fits, but by saying precisely where it applies and solving something else where it does not.

The layer is where the vorticity is

There is a second way to describe the boundary layer that is worth having, because it connects to what the ideal theory threw away.

Ideal flow is irrotational: no fluid element is spinning. That is what makes it solvable in closed form, and it is a genuine property of flow that started from rest and never touched a wall.

The no-slip condition manufactures rotation. Fluid at the wall is stationary and fluid just above it is moving, and a velocity that varies across a small distance is exactly what vorticity measures. So the boundary layer is the region where vorticity lives, and the wall is the source that generates it.

That reframes separation rather usefully. Attached flow keeps its vorticity in a thin sheet wrapped tightly round the body, where it disturbs the outer flow very little. Separation is the moment that sheet is released into the interior of the flow, where it forms the wake — and the outer flow, no longer irrotational, stops being described by the exact theory.

Circulation round a wing has the same origin. The bound vortex that lifts an aircraft is vorticity generated at its surface and organised by the sharp trailing edge.

What Prandtl’s decomposition costs

The two-region picture is so successful that its limits are easy to forget, and they are worth stating.

It assumes the layer is thin, which fails at low Reynolds number and near a leading edge, where the layer has not had distance to grow and the local Reynolds number is small.

It assumes the outer flow imposes the pressure, which fails once the layer separates and the outer flow starts responding to the wake instead of to the body. Beyond separation the decomposition has to be iterated, and it may not converge.

It assumes the two regions can be matched, which requires an overlap region where both descriptions hold. That is usually available and is not guaranteed.

So the model that rescued the subject is itself an approximation with a domain, which is the pattern this whole site keeps finding: every model here is excellent somewhere and silent elsewhere, and knowing which is which is most of the skill.

Where the model stops

The layer is thin only at high Reynolds number. At Reynolds number of order one there is no distinction between a layer and the rest of the flow, and the whole picture collapses. That regime is an entirely different world.

Laminar and turbulent layers behave differently. A turbulent boundary layer is thicker, has far more friction, and is much more resistant to separation. The last of those is why a golf ball is dimpled.

The solver here is coarse. It resolves separation and does not resolve the layer’s internal structure.

Nothing here treats transition. Where a layer becomes turbulent is a hard problem, sensitive to roughness, noise and pressure gradient, and it is not attempted.

What it means for a designer

The thin layer is not only an explanatory device. It sets the shape of the design problem.

Wetted area is expensive. Skin friction scales with the surface the flow touches, so a design that minimises frontal area at the cost of a long body may lose. This is why very long streamlined shapes have an optimum rather than getting indefinitely better.

The rear matters more than the front. Where the flow separates is decided by the pressure recovery over the after-body, so a blunt nose with a fine tail beats a fine nose with a blunt tail — which is why raindrops, airships and submarines look the way they do and not like darts.

Surface finish matters, but only in places. Roughness in a favourable gradient does very little. The same roughness near the point where the layer is struggling can trigger early separation, or prevent one.

What the ideal theory predicts, and what happensThe same cylinder, the same free stream. On the left the exact inviscid solution, closing up behind the body and exerting no drag at all. On the right the real flow at the same conditions, separated, with a wake and therefore with drag.ideal flow — closes up, no dragreal flow at Re 100 — separatedleft: exact closed form · right: solved on a gridRe = 100
Fig. 5 The whole design problem in one comparison: the left-hand flow is what a shape would achieve if its boundary layer never gave up, and the right is what happens when it does.
Flow past a cylinder at Re 40A real fluid past a circular cylinder. At low Reynolds number the flow closes up behind the body much as the ideal theory says; as it rises the flow separates and a region of reversed flow appears behind, which is where drag comes from.separatedrecirculation 0.56 Dviscous flow, solved on a coarse grid — the bubble is under-resolvedRe = 40
Fig. 6 And the same story at a lower Reynolds number, where the layer is relatively thicker and separation is milder. The regime decides how severe the problem is.

The ladder from here

Next rungs: what separation is and when it happens; the layer’s thickness derived rather than quoted; skin friction against pressure drag, and which bodies are dominated by which; laminar and turbulent layers compared; and the dimpled golf ball, which is a deliberate manipulation of all of it.

Then across to the Reynolds number, which decides how thin the layer is and therefore whether any of this applies.