Viscosity

What a plate takes with it

Pull a plate out of a bath and it comes out wet. How wet is not set by the plate, the bath or how much liquid there is, but by a competition in a region a fraction of a millimetre long that nobody looking at the plate can see — and the film goes as the two-thirds power of the speed, never as the first.

Worth reading first: The price of a gradient · The size a drop is allowed.

Take a clean glass plate out of a bowl of water and it comes out wet. The film it carries is about eight microns thick if the plate was drawn out at the speed of an unhurried hand, and the interesting question is what decided that number.

Not the plate: it has no length that enters. Not the bath: its depth does not appear. Not how long the plate was in there. The film’s thickness is set by an argument that takes place entirely in the few tenths of a millimetre where the flat film joins the curved surface of the bath, and it is decided there before the plate has cleared the water.

The region that decides how thick the coat is. The dynamic meniscus, in units of the film's own thickness horizontally and vertically. Far up the plate the film is flat and of thickness one; going down towards the bath it curves away and becomes a parabola, which is what a static meniscus looks like from close up. The whole of the physics is in the join, and its curvature far down the plate is the number the coating thickness is made of.
Fig. 1 The region in question, in units of the film’s own thickness both ways. Far up the plate the film is flat and of thickness one. Going down towards the bath it curves away and becomes a parabola, which is what a static meniscus looks like from close up. The whole of the physics is in the join.

Two things pulling in opposite directions

The plate drags liquid upwards by viscosity: no slip means the fluid touching it moves at the plate’s speed, and it takes its neighbours with it. Surface tension pulls the liquid back down, because a thick film on a plate has more area than a thin one and a curved surface has a pressure jump across it.

The ratio of the two is the capillary number, Ca=μU/γ\mathrm{Ca} = \mu U/\gamma — viscosity times speed over surface tension, dimensionless, and small for almost everything. Water at a centimetre a second gives Ca=1.4×104\mathrm{Ca} = 1.4\times10^{-4}.

The only length in the problem is the capillary length, c=γ/ρg\ell_c = \sqrt{\gamma/\rho g}, which is 2.7 mm for water and is the size at which surface tension and gravity balance. It is the size of a raindrop on a window, the height water climbs in a wettable tube, and — since it is the only length available — it must be what the film’s thickness is measured against.

So the answer has to have the form h=c×f(Ca)h = \ell_c \times f(\mathrm{Ca}), and the whole calculation is finding ff.

The exponent is not fitted

Here is the argument in the form Landau and Levich gave it in 1942.

Far up the plate the film is flat, moves with the plate, and is not being sheared at all. Down at the bath the surface is a static meniscus: gravity against surface tension, no motion, a shape with a definite curvature 2/c\sqrt2/\ell_c where it meets the flat wall.

Between them is a dynamic meniscus in which both matter, and lubrication theory with a capillary pressure applies. Scaling the film by its own thickness h0h_0 and the distance along the plate by h0Ca1/3h_0\mathrm{Ca}^{-1/3} — which is the only combination that makes the viscous and capillary terms the same size — reduces the whole thing to one equation with no parameters at all:

η=3(η1)η3.\eta''' = \frac{3(\eta-1)}{\eta^3}.

Integrate it downstream from the flat film and its curvature approaches a constant. Set that constant equal to the static meniscus’s, undo the scalings, and the thickness falls out. The two-thirds comes from the scaling and nothing else; it was fixed before the equation was solved.

A constant arriving. The curvature of the transition profile against distance down the plate, from three start amplitudes four orders of magnitude apart. All three converge on the same value, which is the check that the number belongs to the equation rather than to where the integration was started — the standard way a shooting calculation produces a confident wrong answer. The three curves are the same curve shifted, because the start amplitude only moves the origin.
Fig. 2 The constant arriving. Three integrations started at amplitudes four orders of magnitude apart, all converging on the same curvature — which is the check that the number belongs to the equation rather than to where the integration was started. The start amplitude only shifts the origin, which is why the three curves are the same curve translated.

The number

The curvature comes out at 1.337577, and matching to a static meniscus of curvature 2/c\sqrt2/\ell_c divides it by 2\sqrt2:

h0=0.94581cCa2/3.h_0 = 0.94581\,\ell_c\,\mathrm{Ca}^{2/3}.

Nothing in that was quoted. The 1.3376 is a Runge–Kutta march, checked for independence of its own starting condition, and the eigenvalue it starts on — the growing solution near the flat film, with k3=3k^3 = 3 — is checked too. The other two roots of k3=3k^3 = 3 are complex with negative real part, which means the approach to the flat film is oscillatory going upstream: there is a ripple ahead of the meniscus, and it is real, and it is why a freshly dip-coated plate has a faint corrugation near the waterline.

Twice as fast is not twice as thick. The film a plate carries out of water, against the speed it is withdrawn at, both logarithmic. The slope is exactly two-thirds, so doubling the speed thickens the film by 58.7 per cent and never by more. The open marks are outside the range the derivation holds in — above a capillary number of about a hundredth the film is no longer thin against the capillary length, and the measured thickness leaves this line.
Fig. 3 The two-thirds law across four decades of withdrawal speed, for water. Doubling the speed thickens the film by 58.7 per cent and never by more. The open marks are outside the range the derivation holds in — above a capillary number of about a hundredth the film is no longer thin against the capillary length, and the measurement leaves this line.

Where the work went, which is not where the film is

This essay sits in a phase about dissipation, and the reason is a fact about the film that its appearance conceals completely.

The flat film costs nothing. It is moving with the plate, uniformly, with no velocity gradient across it at all — a rigid translation, which by the null space of the dissipation function is exactly free. A metre of coated plate carrying eight microns of water is not being sheared anywhere.

Every joule of the extra force needed to withdraw the plate is spent in the dynamic meniscus, in the few tenths of a millimetre where the film is being formed. Integrating the wall stress through that region gives a force per unit width of γCa2/3\gamma\,\mathrm{Ca}^{2/3} times a pure number of about 3.8 — so the withdrawal force has the same two-thirds scaling as the thickness, and it is all incurred in a region a reader would never think to look at.

The film is free; the join is not. The extra force per unit width needed to withdraw a plate, over and above holding the static meniscus, against the capillary number. It goes as Ca to the two-thirds, like the thickness, and every bit of it is spent in the transition region: the flat film above moves with the plate, is not sheared at all, and dissipates nothing. A picture of a coated plate shows the film and shows nothing of where the work went.
Fig. 4 The extra force per unit width against the capillary number. It is not quite a constant times Ca2/3\mathrm{Ca}^{2/3}, and the drift is honest rather than numerical: the integrand falls only as one over the square of distance down the plate, so the integral has a tail that depends on where the dynamic meniscus is taken to end — which is where the static one begins, and that is itself a function of Ca.

The same constant in a completely different apparatus

The strongest evidence that 1.3376 belongs to the join rather than to the plate is that it turns up somewhere with no plate in it.

Drive a long bubble along a liquid-filled tube and it leaves a film on the wall behind it. There is no gravity in the problem, no bath, and no capillary length — the geometry supplies its own length, the tube’s radius. Bretherton solved it in 1961 and the answer is

hR=1.3376Ca2/3.\frac{h}{R} = 1.3376\,\mathrm{Ca}^{2/3}.

The same constant. What the two problems share is a transition region between a flat film and a curved static surface, governed by the same parameterless equation; what they do not share is anything else at all.

The same constant, in a completely different apparatus. A long bubble driven along a tube leaves a film on the wall, and its thickness relative to the tube's radius is the same 1.3376 times the capillary number to the two-thirds. Nothing about the two problems is shared except the transition region between a flat film and a curved static surface — one matches a plate to a bath, the other a tube wall to a bubble cap — which is the clearest evidence that the constant belongs to the join and not to either apparatus.
Fig. 5 Bretherton’s film. The matching is to a hemispherical cap of curvature 2/R rather than to a meniscus of curvature √2/ℓ_c, so the two coefficients differ by a factor that is pure geometry — and the number in the middle is identical.

Why there is no length in the answer

The absence of the plate’s own dimensions from the result is worth dwelling on, because it is the feature that makes the law useful and it is not obvious that it should hold.

A plate a metre wide and a plate a centimetre wide carry the same film. A plate withdrawn from a bucket and a plate withdrawn from a lake carry the same film. A plate that has been in the water for an hour and one dipped for a second carry the same film. None of that is true of, say, the amount of liquid that drips off afterwards, which depends on the plate’s height and on how long it is left.

The reason is that the film is decided in a region whose size is h0Ca1/3h_0\mathrm{Ca}^{-1/3} — a fraction of a millimetre for water at ordinary speeds — and everything outside that region enters only through one number, the curvature the static surface presents. The bath is summarised by a curvature and nothing else. That is the whole content of a matched asymptotic expansion, and this is one of the cleanest examples of one in the subject: an outer problem that is entirely static, an inner problem with no parameters at all, and a single number passed between them.

It is also the reason the answer is so robust. Deepen the bath, widen the plate, warm the room: none of those changes the curvature at the wall, and none of them changes the coat.

Six things coming out of six liquids

The law is easy to apply and the range it covers is wider than the arithmetic suggests, because the capillary number is a product of three quantities that each span decades.

Six things coming out of six liquids. The same law applied to ordinary objects, on a logarithmic axis of film thickness. The range is four decades and the objects span it not because they are pulled at very different speeds but because the capillary number does: honey is five thousand times as viscous as water and molten zinc has ten times the surface tension. A finger out of water carries about eight microns; a spoon out of honey carries most of a millimetre.
Fig. 6 Six ordinary objects, spanning four decades of film thickness. A finger out of water carries about eight microns; a spoon out of honey carries most of a millimetre, because honey’s viscosity is five thousand times water’s. A galvanising line pulls steel out of molten zinc at two metres a second and carries about twenty microns, which is the coating weight the whole process is specified by.

Two of those are worth a sentence each. The galvanising line is the case where the law is used in anger: coating weight is a purchased specification, it is controlled by line speed, and the two-thirds exponent means a ten per cent speed change moves the coating by six and a half per cent rather than by ten. And the finger is the case where the law is nearly always misapplied — a finger is not flat, and the curvature of a finger is comparable with the capillary length, so the number above is indicative rather than right.

The energy, done properly

The claim that the flat film is free deserves the same treatment as everything else here, which is to be computed rather than asserted, and doing so turns up a comparison worth having.

Three energies are in play per unit area of plate coated.

The work done against the extra viscous force is γCa2/3\gamma\,\mathrm{Ca}^{2/3} times about 3.8, which for water at a centimetre a second is 0.7 millijoules per square metre.

The surface energy created is γ\gamma times the new area — one new interface, so about 72 millijoules per square metre. A hundred times larger.

And the potential energy of the raised film is hogh2/2 ho g h^2/2, which is a nanojoule per square metre and is negligible against both.

So coating is overwhelmingly a surface-energy transaction with a small viscous overhead, and the overhead is what sets the thickness. That is the resolution of a paradox that would otherwise be awkward: the film’s thickness is decided by a competition between viscosity and surface tension, and the energy is almost entirely surface tension, because the viscous term appears at the two-thirds power of a small number.

Doubling the withdrawal speed raises the viscous overhead by 59 per cent and the surface term not at all. There is no speed at which coating becomes viscously expensive in the sense that a bearing is.

The motions that cost nothing. Four velocity fields and what each of them costs. A uniform translation and a solid-body rotation deform nothing, so their rate of strain is exactly zero and so is their dissipation, at any speed and any spin rate. The last two cost exactly the same as one another — a simple shear of rate γ and a pure strain of rate γ/2 have the same rate of strain — even though the shear's velocity gradient is √2 larger. All of that difference is rotation, and rotation is free.
Fig. 7 And the reason the film itself contributes nothing to the first of those three. A film moving with the plate is a rigid translation, its rate of strain is identically zero, and Φ = 2μe:e is zero with it — at any thickness, at any speed, for any fluid.

What the picture cannot show

The plate is perfectly wetted. The derivation assumes the liquid meets the solid at zero contact angle, so that the static meniscus runs smoothly onto the film. A partially wetting liquid has a contact line, the whole matching argument is different, and below a critical speed the plate comes out dry — a transition with no counterpart in the equations above.

The film is thin against the capillary length and the derivation says nothing beyond that. Above Ca102\mathrm{Ca}\approx10^{-2} the measured thickness leaves the two-thirds line, gravity begins to drain the film while it is still being formed, and the correct treatment carries a gravity term the scaling above dropped.

The flow is creeping. The film Reynolds number is returned by the solver and checked, and every number quoted here is from a case where it is below one — the same restriction this collection’s other slow-flow essays run under.

The film is uniform and real ones are not. A withdrawn film is unstable to long-wavelength disturbances and to drainage, and what is measured a second later is not what was deposited. The number here is the thickness at the moment of formation.

And there is no evaporation, no solvent and no drying. Every industrial coating is a solution that loses most of its volume after deposition, and the relation between the wet film computed here and the dry one that is sold is a chemistry problem.

The depth a puddle cannot exceed. Far from its edge a wide sheet of liquid is flat, so one of its two curvatures vanishes and the equation reduces to sinφ dφ = z̄ dz̄ — which integrates with no size in it at all. The depth is 2ℓsin(θ/2), and for water on a perfectly non-wetting surface that is 5.4 mm however much is poured on. The marks are the integrator run on wide ridges at six contact angles; they agree with the closed form to a part in four thousand.
Fig. 8 The static half of the matching, drawn on its own. A pendant shape held by surface tension against gravity has a curvature at the wall that is set by the capillary length, and that curvature is the only thing the dynamic region needs to know about the bath. Everything else about the bath is irrelevant to how wet the plate comes out.

What a coater controls, and what it cannot

The two-thirds exponent has one consequence for process control that is worth stating plainly, because it runs the other way from intuition.

A coating line’s thickness is set by speed, and speed is the variable a line operator has. Because the exponent is less than one, the film is less sensitive to speed than it would be if the relation were proportional — a one per cent speed error gives a two-thirds of one per cent thickness error. That is a helpful direction to be wrong in and it is the reason dip coating survives as a process at all.

The unhelpful direction is viscosity, which enters at the same two-thirds power and is far harder to hold. A bath that warms by ten kelvin over a shift loses a quarter of its viscosity, and the coating thins by seventeen per cent — an enormous drift, invisible to anyone watching the line, and the reason coating baths are temperature-controlled to a fraction of a degree.

And surface tension enters at minus two-thirds through the capillary number and at plus one-half through the capillary length, so it appears overall as γ1/6\gamma^{-1/6} — a sixth power, which is almost no dependence at all. Contamination that halves a bath’s surface tension changes the film by twelve per cent. Of the three properties, the one a chemist would worry about first matters least.

The ceiling on the one variable a line has

The two-thirds law says thickness rises with speed, and it says nothing about how fast a line may be run. That limit exists, it is sharp, and it comes from the one part of the problem this derivation assumed away: the contact line.

The derivation above takes the liquid to wet the solid perfectly, so there is no contact line — the static meniscus runs continuously onto the film. A partially wetting liquid has one, where solid, liquid and air meet, and it has to travel along the plate at the withdrawal speed. It manages that by deforming: the apparent contact angle grows with the capillary number, because the liquid wedge near the line has to accommodate a viscous stress that diverges as the wedge thins.

An angle cannot grow past 180°. When it reaches that, the line can no longer keep up, and what happens is that the air wins the race: a film of it is dragged in under the advancing liquid. That is air entrainment, and it is the failure mode every coating process is really limited by. The entrained film breaks into bubbles, and a coating with bubbles in it is scrap.

So the practical maximum speed of a coating line is a wetting instability rather than anything in the film equation, and it is why the two-thirds law is used to choose a speed inside a window whose upper edge the law says nothing about. Higher-viscosity liquids entrain at lower speeds, which sets the awkward trade of the whole business: viscosity is the cheapest way to thicken a coating and the surest way to lower the ceiling on how fast it can be laid.

Two responses to that are worth knowing, because both are counter-intuitive and both are used.

Change the gas. The entrained thing is air, so its properties matter. Running the coating station at reduced pressure, or in a helium atmosphere, delays entrainment substantially and buys real line speed — an intervention with no liquid, no plate and no film in it.

Or stop dipping. Slot-die and curtain coating do not rely on a free contact line advancing against a static bath; a curtain arrives with momentum of its own and helps the line along. Both reach speeds dip coating cannot, and both give up the property that makes dip coating attractive in the first place — that its thickness is set by a law with nothing adjustable in it.

Who found it, and when

Landau and Levich published in 1942 — five years before the boundary-layer energy equation this collection uses in the third thickness, and by an entirely separate route in the Soviet physics literature, and Derjaguin gave an equivalent treatment the following year, which is why the result carries three names in different places. Bretherton’s tube version came in 1961 and the recognition that the two share a constant came later still.

The surprising connection is with a problem in a different subject that shares the exponent as well as the mathematics. The advance of a wetting front, the spreading of a droplet, the shape a drop takes when it is falling, the rise of a film in a corner and the coating of a plate are all governed by a balance between viscous stress and capillary pressure in a thin wedge, and all of them produce exponents of the form Ca1/3\mathrm{Ca}^{1/3} or Ca2/3\mathrm{Ca}^{2/3} for the same reason: the only way to make a viscous term and a third derivative of a thickness the same size is to scale the length by the cube root of the capillary number. The two-thirds is a consequence of a third derivative, and the third derivative is a consequence of Laplace pressure being a curvature.

Where the ladder goes next

Below this rung is the capillary length itself, which is the only length in the problem, and the price of a gradient, which is what makes the claim that the flat film is free into a computation rather than an assertion.

Beside it is the drop that is not a tear, which is the same competition between surface tension and something else with the shape as the unknown, and the film that will not be squeezed out, which is a thin film whose thickness is set by a completely different balance and whose behaviour is correspondingly nothing like this one’s.

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.

Capillary lengthCapillary numberCoatingDissipationDissipation functionLubrication filmMatched asymptoticsSurface tensionThin filmViscosity