courseshub.world · Lives in Mathematics and Physics

Osborne Reynolds

Osborne Reynolds, FRS

A course in two strands: the life of the clergyman’s son who served an engineering apprenticeship before Cambridge and became Manchester’s first professor of engineering, and the physics of heat transfer, turbulence, lubrication and granular matter that he founded.

Born 23 August 1842, Belfast Died 21 February 1912, Watchet, Somerset Professor of Engineering Owens College, Manchester, 1868–1905 Level advanced undergraduate → graduate
The two manners of motion (1883). Top: a dye filament in slow and fast pipe flow. Bottom: friction factor against Re; the laminar law 64/Re gives way near Re ≈ 2000 to the turbulent Blasius law.
Lesson 00

How to use this course

Every engineering student meets Osborne Reynolds through one dimensionless number. That number is the smallest part of his legacy. Reynolds explained why hot fluids cool faster when they flow faster, why a film of oil can carry a railway axle, why wet sand dries around a footprint, and how the chaotic eddies of a turbulent flow exert a stress on the mean motion. He did this as a professor of engineering, with a workshop training and a Cambridge mathematical education, at a time when engineering was barely an academic subject.

Prerequisites for Part II
  • Vector calculus and ordinary differential equations.
  • The Navier–Stokes equations for an incompressible fluid, at the level of the Stokes course.
  • Elementary thermodynamics and kinetic theory for Lesson 04.
  • Basic probability (means and variances) for Lesson 09.

Each technical lesson begins with a paragraph for the general reader, then gives precise statements, derivations and exercises with folded solutions. The order roughly follows Reynolds’s own: heat transfer and gases (1874–1879), pipe flow and similarity (1883–1887), lubrication (1886), granular matter (1885), and the averaged equations and energy method of 1895.

Reference

Timeline

1842Born 23 August in Belfast, son of the Rev. Osborne Reynolds, headmaster of the Belfast Collegiate School. The family soon moves to Dedham, Essex.
1861Apprenticed to the engineer Edward Hayes at Stony Stratford.
1867Seventh wrangler in the Mathematical Tripos; fellow of Queens’ College, Cambridge. A year as a civil engineer with Lawson and Mansergh in London.
1868Appointed to the new chair of engineering at Owens College, Manchester, aged 25. Marries Charlotte Jemima Chadwick.
1869Elected to the Manchester Literary and Philosophical Society. His first wife dies after the birth of their son.
1874Paper on the heating surface of steam boilers: the Reynolds analogy. Secretary of the Lit & Phil, 1874–1884.
1877Fellow of the Royal Society.
1879Long memoir on thermal transpiration and the dynamical theory of gases in the Philosophical Transactions.
1883March: submits the paper on direct and sinuous flow in pipes, the origin of the Reynolds number.
1884Honorary LLD, Glasgow. Candidate for the Cavendish chair at Cambridge, which goes to J. J. Thomson.
1885Names and explains dilatancy in granular media.
1886On the theory of lubrication, explaining Beauchamp Tower’s bearing experiments.
1888Royal Medal of the Royal Society. President of the Lit & Phil.
1892Publishes his memoir of James Prescott Joule.
1895Paper on the dynamical theory of incompressible viscous fluids: Reynolds averaging, the Reynolds stresses and the energy method.
1897Bakerian Lecture of the Royal Society.
1900–03Papers on Mechanical and Physical Subjects, three volumes.
1902–03Rede Lecture at Cambridge; The Sub-Mechanics of the Universe. Dalton Medal, 1903.
1905Retires because of failing health.
1912Dies of influenza on 21 February at Watchet, Somerset; buried at St Decuman’s church.
Part I · Formation
Lesson 01

Belfast, Dedham and the workshop, 1842–1863

Osborne Reynolds was born in Belfast on 23 August 1842. His father, also Osborne, was an Anglican clergyman and a Cambridge graduate of 1837 who had been a fellow of Queens’ College. At the time of his son’s birth he was headmaster of the Belfast Collegiate School. Three generations of the family had been rectors of Debach-with-Boulge in Suffolk. Soon after Osborne’s birth the family moved to Dedham in Essex, where the father became headmaster of the grammar school. A younger brother, Edward, became a priest.

The boy was taught first at Dedham and then privately, and above all by his father, a capable mathematician who also took out patents on agricultural machinery. Reynolds later wrote that he had enjoyed his father’s constant guidance in mechanics through his boyhood. The combination of mathematics and machines stayed with him for life.

The apprenticeship

Instead of going directly to university, Reynolds in 1861 became an apprentice in the workshop of Edward Hayes, a mechanical engineer at Stony Stratford in Buckinghamshire who built engines and small coastal steamers. There he learned how machines were made and how they failed. He came to believe that the practical questions he met there could only be answered by mathematics, and that conviction sent him to Cambridge.

DiscussionReynolds reversed the usual order of training, putting the workshop before the university. What does a practical apprenticeship teach that a mathematical education does not, and how might this order have shaped the problems he later chose?

Lesson 02

Cambridge and the London sewers, 1863–1868

At Queens’ College, Cambridge, his father’s college, Reynolds read mathematics. The Mathematical Tripos of the 1860s was a hard competitive examination in mathematics and mathematical physics, taught through private coaches. John William Strutt, later Lord Rayleigh, was a year ahead of him, and they attended some classes together. In 1867 Reynolds graduated as seventh wrangler, that is, seventh in the first class, and was elected to a fellowship at Queens’.

He did not stay to teach. He spent a year as a practising civil engineer with the London firm of Lawson and Mansergh, working on the sewage system for Croydon. Moving large volumes of water through pipes and channels was the bread and butter of Victorian civil engineering. The resistance of pipes, the subject of his most famous paper, was in front of him from the start.

Engineering in the universities

In the 1860s engineering was learned by apprenticeship. Glasgow had had a chair of engineering since 1840, held from 1855 by W. J. M. Rankine, and University College London offered engineering courses. In England outside London there was nothing. When the industrialists of Manchester endowed a chair at Owens College in 1868, they created one of the first university posts in the country with the title Professor of Engineering.

Lesson 03

The Manchester chair, 1868–1880

The post at Owens College was advertised in 1868 at £500 a year. Reynolds applied, writing that he had an irresistible liking for mechanics and the physical laws on which it rests. He was appointed at 25 and stayed for 37 years. Owens College became a constituent of the new federal Victoria University in 1880 and later the University of Manchester. Reynolds built its engineering school and, later, the Whitworth engineering laboratories.

In 1868 he married Charlotte Jemima Chadwick. She died in 1869 after the birth of their son Osborne, who died in childhood in 1879. In the early 1880s Reynolds married Annie Charlotte Wilkinson; they had three sons and a daughter.

From electricity to fluids

His first research, from 1868 to about 1873, was on electricity and magnetism and on electrical theories of the Sun and comets. Around 1873 he turned to mechanics, and for more than two decades fluid motion was his main subject. He published in the Memoirs and Proceedings of the Manchester Literary and Philosophical Society, which he joined in 1869 and served as secretary from 1874 to 1884, and increasingly in the journals of the Royal Society, which elected him a fellow in 1877.

His teaching was demanding. Students found his lectures hard to follow, because he worked out ideas at the board rather than delivering polished notes. J. J. Thomson, who studied at Owens College in the 1870s before going to Cambridge, later recalled him as a teacher in his memoirs. Horace Lamb, professor of mathematics at Owens from 1885 and Reynolds’s obituarist, remarked that he had no taste for advertisement.

Part II · The physics
Lesson 04

Boilers, gases and the Reynolds analogy

For the general reader. Blowing on soup cools it faster. Reynolds explained why: the same swirling motions that drag a fast fluid against a wall also carry heat between the fluid and the wall. So if you know the friction, you know the heat transfer. This idea, from a paper on steam boilers, is still the first estimate engineers make.

Goals
  • Derive the Reynolds analogy between wall friction and heat transfer.
  • Define the Stanton number and skin-friction coefficient and see where the analogy fails.
  • Understand thermal transpiration, the effect Reynolds studied in his 1879 memoir on gases.

The analogy (1874)

Reynolds’s 1874 paper on the heating surface of steam boilers argued that heat passes between a flowing fluid and a wall in two ways: by molecular conduction, and by the mixing motions of the fluid itself, which carry hot fluid to the wall and cold fluid away. The mixing part grows with speed. The same mixing transports momentum toward the wall and appears as friction.

Model both transports in a turbulent layer next to a wall by a common eddy diffusivity \(\varepsilon\). With \(y\) the distance from the wall, \(u(y)\) the mean velocity and \(T(y)\) the mean temperature,

\[ \tau=\rho\,\varepsilon\,\frac{du}{dy},\qquad q=-\rho c_p\,\varepsilon\,\frac{dT}{dy}. \]

If \(\tau\) and \(q\) vary across the layer in the same proportion, then \(q/\tau=-c_p\,dT/du\) is constant, and integrating from the wall (\(u=0,\ T=T_w\)) to the outer flow (\(u=U,\ T=T_\infty\)) gives \(q_w/\tau_w=c_p(T_w-T_\infty)/U\).

Reynolds analogy

Define the heat-transfer coefficient \(h=q_w/(T_w-T_\infty)\), the Stanton number \(\mathrm{St}=h/(\rho c_p U)\), and the skin-friction coefficient \(C_f=\tau_w/(\tfrac12\rho U^2)\). If momentum and heat are transported by the same mechanism, with equal diffusivities, then

\[ \mathrm{St}=\tfrac12 C_f. \]

The hypothesis of equal diffusivities holds near a wall only when the molecular diffusivities are also equal, that is, when the Prandtl number \(\mathrm{Pr}=\nu/\kappa\) is close to 1, as for gases. For other fluids Thomas Chilton and Allan Colburn proposed in 1933 the empirical correction \(\mathrm{St}\,\mathrm{Pr}^{2/3}=\tfrac12C_f\), valid roughly for \(0.6<\mathrm{Pr}<60\). The analogy also fails when pressure gradients change the friction without changing the heat flux in the same way, for example in separated flow.

Thermal transpiration (1879)

The radiometer of William Crookes (1873), a vane that spins in a partial vacuum when lit, provoked a long argument about its mechanism. Reynolds traced the force to gas flowing along temperature gradients at surfaces, and in 1879 published a memoir on thermal transpiration: when two vessels at temperatures \(T_1\) and \(T_2\) are joined by a porous plug whose pores are narrow compared with the mean free path, gas creeps from the cold side to the hot side until a pressure difference builds up. James Clerk Maxwell refereed the memoir and wrote his own paper on rarefied gases in response. In the limit of very narrow pores, a steady state requires equal molecular fluxes through the pores in both directions, and the flux of molecules through a small hole is proportional to \(p/\sqrt{T}\). Hence

\[ \frac{p_1}{p_2}=\sqrt{\frac{T_1}{T_2}}\qquad\text{(free-molecular limit).} \]

For wide pores the pressures equalize. Reynolds measured the transition between these regimes as the ratio of pore size to mean free path changed, an early systematic study of what is now called rarefied gas dynamics.

Exercise 4.1Air flows at \(U=10\ \mathrm{m\,s^{-1}}\) over a flat plate, with \(\rho=1.2\ \mathrm{kg\,m^{-3}}\), \(c_p=1005\ \mathrm{J\,kg^{-1}K^{-1}}\). The measured skin-friction coefficient is \(C_f=0.004\). Estimate the heat-transfer coefficient \(h\) by the Reynolds analogy, and the heat flux when the plate is 30 K hotter than the air.

Show solution
\(\mathrm{St}=0.002\), so \(h=\mathrm{St}\,\rho c_pU=0.002\times1.2\times1005\times10\approx24\ \mathrm{W\,m^{-2}K^{-1}}\). The flux is \(q\approx24\times30\approx720\ \mathrm{W\,m^{-2}}\). Since \(\mathrm{Pr}\approx0.71\) for air, the Chilton–Colburn form would raise this by \(\mathrm{Pr}^{-2/3}\approx1.26\), to about 30 \(\mathrm{W\,m^{-2}K^{-1}}\).

Exercise 4.2Derive the flux of molecules through a small hole of area \(A\) in a wall, \(\Phi=\tfrac14 n\bar c A\), where \(\bar c=\sqrt{8k_BT/\pi m}\) is the mean speed, and deduce the free-molecular transpiration law.

Show solution
Molecules with speed \(c\) moving at angle \(\theta\) to the normal cross area \(A\) at rate \(n f(c)\,c\cos\theta\,A\) per unit of the speed distribution and solid angle fraction \(\sin\theta\,d\theta\,d\phi/4\pi\). Integrating \(\cos\theta\sin\theta\) over the forward hemisphere gives \(\pi\), so \(\Phi=nA\bar c\cdot\pi/4\pi=\tfrac14n\bar cA\). With \(n=p/k_BT\), \(\Phi\propto p/\sqrt{T}\). Equal fluxes in both directions give \(p_1/\sqrt{T_1}=p_2/\sqrt{T_2}\).
Lesson 05

The two manners of motion

For the general reader. Turn on a tap gently and the water runs in a smooth glassy stream. Open it further and the stream becomes ragged. In 1883 Reynolds showed with a thread of dye in a glass pipe that water flows in two different ways, smoothly or chaotically, and that one number built from the speed, the pipe width and the viscosity decides which.

Goals
  • Describe Reynolds’s dye experiment and its findings.
  • Derive the Reynolds number from the Navier–Stokes equations.
  • Compare the laminar and turbulent laws of pipe resistance.

The experiment

The paper, submitted to the Royal Society in March 1883 and published in the Philosophical Transactions, has a long title: an experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels. Reynolds mounted glass tubes with flared, bell-mouthed entrances inside a tank of still water and drew water through them with a valve at the far end. A fine jet of coloured water was released at the centre of the entrance.

At low speeds the dye formed a straight thread along the whole tube. As the speed increased, at a definite point the thread suddenly mixed into the water and filled the tube with colour. Lit by an electric spark, the coloured cloud resolved into distinct eddies. Reynolds called the two regimes direct and sinuous; we say laminar and turbulent. He found that the critical speed was inversely proportional to the tube diameter and proportional to the kinematic viscosity, which he varied by changing the water temperature. The critical combination was therefore a pure number.

He found that the value at which turbulence appeared depended on how calm the water was in the tank, and could be raised far above its usual value with great care. On the other hand, a turbulent flow slowed down gradually always became laminar below a lower value, close to 2000 in modern units. Precise measurements in 2011 placed the threshold for sustained pipe turbulence at about 2040.

Dynamical similarity

Write the incompressible Navier–Stokes equations in terms of a velocity scale \(U\), a length \(L\), time \(L/U\) and pressure \(\rho U^2\). With \(\mathbf u=U\mathbf u^*\), \(\mathbf x=L\mathbf x^*\), and so on,

\[ \frac{\partial\mathbf u^*}{\partial t^*}+(\mathbf u^*\cdot\nabla^*)\mathbf u^*=-\nabla^*p^*+\frac{1}{\mathrm{Re}}\nabla^{*2}\mathbf u^*,\qquad \nabla^*\!\cdot\mathbf u^*=0,\qquad \mathrm{Re}=\frac{\rho UL}{\mu}=\frac{UL}{\nu}. \]
Principle of dynamical similarity

Two incompressible flows in geometrically similar domains, with similar boundary conditions, satisfy the same dimensionless equations if and only if they have the same Reynolds number. Their velocity fields are then the same up to scaling, and every dimensionless output, such as a friction factor, is a function of \(\mathrm{Re}\) alone.

George Gabriel Stokes had identified this combination in his 1851 paper on pendulums. Reynolds showed experimentally that it controls the most important qualitative change a flow can undergo. Arnold Sommerfeld attached Reynolds’s name to it in 1908.

Laws of resistance

For a pipe of diameter \(D\), mean speed \(V\) and length \(\ell\), the pressure drop is written \(\Delta p=f\,\frac{\ell}{D}\,\frac12\rho V^2\), with \(f\) the Darcy friction factor. In the second half of his paper Reynolds measured the resistance of long pipes and found that it rises in proportion to the speed below the critical value, as the Poiseuille law requires, and roughly as a power of the speed close to 2 above it.

Two regimes of pipe resistance

Laminar (Hagen–Poiseuille): \(f=64/\mathrm{Re}\), with \(\mathrm{Re}=VD/\nu\).
Turbulent, smooth pipe (Paul Blasius, 1913, empirical for \(4\times10^3\lesssim\mathrm{Re}\lesssim10^5\)): \(f\approx0.316\,\mathrm{Re}^{-1/4}\), so \(\Delta p\propto V^{7/4}\).

Exercise 5.1From the Poiseuille profile \(w(r)=\frac{G}{4\mu}(R^2-r^2)\), show that \(f=64/\mathrm{Re}\).

Show solution
The flux is \(Q=\pi GR^4/8\mu\), so the mean speed is \(V=Q/\pi R^2=GR^2/8\mu=GD^2/32\mu\). Hence \(G=\Delta p/\ell=32\mu V/D^2\). Then \(f=\frac{\Delta p}{\ell}\frac{D}{\frac12\rho V^2}=\frac{32\mu V}{D^2}\cdot\frac{2D}{\rho V^2}=\frac{64\mu}{\rho VD}=\frac{64}{\mathrm{Re}}\).

Exercise 5.2Water (\(\nu=1.0\times10^{-6}\ \mathrm{m^2s^{-1}}\)) flows through a domestic pipe of diameter 15 mm. Find the mean speed at which \(\mathrm{Re}=2000\). Is a tap delivering 6 litres per minute through this pipe laminar?

Show solution
\(V=\mathrm{Re}\,\nu/D=2000\times10^{-6}/0.015\approx0.13\ \mathrm{m\,s^{-1}}\). A flow of 6 L/min is \(10^{-4}\ \mathrm{m^3s^{-1}}\); the cross-section is \(1.77\times10^{-4}\ \mathrm{m^2}\), so \(V\approx0.57\ \mathrm{m\,s^{-1}}\) and \(\mathrm{Re}\approx8500\): turbulent.

Exercise 5.3At what Reynolds number do the laminar and Blasius laws give the same friction factor? Why is that value not the transition point?

Show solution
Solve \(64/\mathrm{Re}=0.316\,\mathrm{Re}^{-1/4}\): \(\mathrm{Re}^{3/4}=202.5\), so \(\mathrm{Re}\approx1200\). Transition is not decided by which law gives less resistance. Laminar flow is a solution at every \(\mathrm{Re}\) and is linearly stable in a pipe; turbulence needs finite disturbances and sustains itself only above about 2000. Between the curves’ intersection and transition, laminar flow simply persists.
Lesson 06

Dynamical similarity and models

For the general reader. Engineers test small model ships, aircraft and harbours before building the real thing. A model only tells the truth if the right dimensionless numbers match. Reynolds was among the first to use scale models systematically in research, building small working models of tidal estuaries to study how sand banks form.

Goals
  • State the Buckingham Π theorem and apply it to pipe flow.
  • Understand why Froude and Reynolds scaling conflict in ship testing.
  • See how Reynolds used models of estuaries.
Π theorem (Buckingham, 1914)

If a physical relation \(F(q_1,\dots,q_n)=0\) among \(n\) dimensional quantities is independent of the choice of units, and the dimension matrix of the \(q_i\) has rank \(r\), then it is equivalent to a relation \(\Phi(\Pi_1,\dots,\Pi_{n-r})=0\) among \(n-r\) independent dimensionless products of the \(q_i\).

Proof idea. Dimensions form a vector space over \(\mathbb Q\) with basis M, L, T. Monomials \(\prod q_i^{a_i}\) are dimensionless when \((a_i)\) lies in the kernel of the \(r\times n\) dimension matrix, a space of dimension \(n-r\). Rescaling the base units acts on the remaining \(r\) quantities freely, so a unit-independent relation can depend only on the invariants.

Froude against Reynolds

William Froude’s towing-tank experiments of the 1870s showed that wave-making resistance of ships scales with the Froude number \(\mathrm{Fr}=U/\sqrt{gL}\). Frictional resistance, which Reynolds’s work showed depends on \(\mathrm{Re}\), scales differently. Both cannot be matched in water, so model tests split the resistance into a wave part scaled by Froude and a friction part estimated from flat-plate data as a function of \(\mathrm{Re}\). Reynolds wrote on ship design and screw propulsion for the Institution of Naval Architects. His studies of the racing of screw propellers, when the blades lose their grip on the water, touched on what is now called cavitation, the formation of vapour cavities where the local pressure falls too low.

Estuary models

In the late 1880s Reynolds built small tidal models, with sand beds and a mechanism to raise and lower the sea level periodically, to reproduce the banks and channels of estuaries. He reported them in his 1887 address to the engineering section of the British Association. Such models cannot match \(\mathrm{Re}\), yet he found that they reproduced the main features of real estuaries. Hydraulic laboratories later used the same approach for rivers and harbours throughout the twentieth century.

Exercise 6.1The pressure drop per unit length \(G\) in a pipe depends on \(\rho,\mu,V,D\) and the wall roughness height \(k\). How many dimensionless groups are there? Write the relation in standard form.

Show solution
Six quantities, three base dimensions of full rank, so three groups: \(\frac{GD}{\rho V^2}=\phi\big(\frac{\rho VD}{\mu},\frac kD\big)\), that is, \(f=f(\mathrm{Re},k/D)\). This is the content of the Moody diagram.

Exercise 6.2A ship of length \(L_p\) is tested with a model of length \(L_m=L_p/25\). Find the model speed required for Froude similarity, and the kinematic viscosity a model fluid would need to match the Reynolds number as well.

Show solution
Equal \(\mathrm{Fr}\) gives \(U_m=U_p\sqrt{L_m/L_p}=U_p/5\). Equal \(\mathrm{Re}\) then needs \(\nu_m=\nu_p\frac{U_mL_m}{U_pL_p}=\nu_p(L_m/L_p)^{3/2}=\nu_p/125\). No practical liquid is 125 times less viscous than water, so both cannot be matched.
Lesson 07

The theory of lubrication

For the general reader. A spinning shaft in an oiled bearing does not touch the metal around it. It floats on a thin film of oil whose pressure can exceed many atmospheres, generated by the motion itself. Reynolds explained this in 1886 with an equation that is still the basis of bearing design and of the whole science of tribology.

Goals
  • Derive the Reynolds equation of thin-film lubrication.
  • Solve the inclined slider bearing and compute its load capacity.
  • Understand why a converging gap is needed.

Tower’s experiments

In 1883 Beauchamp Tower, working for the Institution of Mechanical Engineers, tested a railway axle bearing running in an oil bath. When he drilled a hole in the bearing cap for an oiler, oil was forced out of it; a pressure gauge on the hole showed a pressure far above the average load per unit area. The oil was supporting the load hydrodynamically. Reynolds’s paper On the theory of lubrication and its application to Mr. Beauchamp Tower’s experiments, in the Philosophical Transactions of 1886, explained this and also reported his measurements of the viscosity of olive oil as a function of temperature.

Thin-film equations

Let the film occupy \(0<z<h(x)\), with \(h\) small compared with the length \(L\) of the bearing, the lower surface sliding with speed \(U\) in the \(x\) direction and the upper surface fixed. If \(\mathrm{Re}\,(h/L)^2\ll1\), inertia is negligible and the Navier–Stokes equations reduce to

\[ \frac{\partial p}{\partial x}=\mu\frac{\partial^2u}{\partial z^2},\qquad \frac{\partial p}{\partial z}=0. \]

So \(p=p(x)\), and with \(u(0)=U,\ u(h)=0\),

\[ u=U\Big(1-\frac zh\Big)-\frac{1}{2\mu}\frac{dp}{dx}\,z(h-z),\qquad q=\int_0^hu\,dz=\frac{Uh}{2}-\frac{h^3}{12\mu}\frac{dp}{dx}. \]
Reynolds equation (1886), one-dimensional form

For a steady incompressible film the flux \(q\) is independent of \(x\), and differentiating gives

\[ \frac{d}{dx}\Big(h^3\frac{dp}{dx}\Big)=6\mu U\frac{dh}{dx}. \]

In two dimensions, with surface velocity \((U,V)\) and a gap that may change in time, \(\nabla\!\cdot\big(h^3\nabla p\big)=6\mu\big(U\partial_xh+V\partial_yh\big)+12\mu\,\partial_th\).

The slider bearing

Take \(h\) decreasing linearly from \(h_1\) at \(x=0\) to \(h_2\) at \(x=L\), with \(p=0\) at both ends. Integrating once, \(dp/dx=6\mu U(h-h^*)/h^3\), where \(h^*\) is the gap at the pressure maximum. The boundary conditions give \(h^*=2h_1h_2/(h_1+h_2)\) and

\[ p(x)=\frac{6\mu UL}{h_1^2-h_2^2}\,\frac{(h_1-h)(h-h_2)}{h^2}. \]

With \(k=h_1/h_2\), the load per unit width is

\[ W=\int_0^Lp\,dx=\frac{6\mu UL^2}{(k-1)^2h_2^2}\Big(\ln k-\frac{2(k-1)}{k+1}\Big). \]

The load grows as \(1/h_2^2\): a heavier load squeezes the film thinner, which raises the pressure until it balances. This self-regulation is why hydrodynamic bearings are stable. The bracket is maximal at \(k\approx2.2\), where \(W\approx0.16\,\mu UL^2/h_2^2\).

Exercise 7.1Show that a parallel gap (\(h\) constant) with \(p=0\) at both ends carries no load. What does this say about the role of the wedge?

Show solution
With \(dh/dx=0\), \(h^3p''=0\), so \(p\) is linear; vanishing at both ends forces \(p\equiv0\). Pressure is generated only where the sliding surface drags fluid into a narrowing gap. If the gap widens in the direction of motion the computed pressure is negative, and in practice the film ruptures or cavitates.

Exercise 7.2Verify the formula \(h^*=2h_1h_2/(h_1+h_2)\). Hint: integrate \(dp/dx=6\mu U(h-h^*)/h^3\) with \(dx=-L\,dh/(h_1-h_2)\) between \(h_1\) and \(h_2\) and set the result to zero.

Show solution
\(\int_{h_2}^{h_1}(h^{-2}-h^*h^{-3})\,dh=\big(\tfrac1{h_2}-\tfrac1{h_1}\big)-\tfrac{h^*}{2}\big(\tfrac1{h_2^2}-\tfrac1{h_1^2}\big)=0\). Dividing by \(\tfrac1{h_2}-\tfrac1{h_1}\) gives \(1=\tfrac{h^*}2\big(\tfrac1{h_2}+\tfrac1{h_1}\big)\), so \(h^*=\frac{2h_1h_2}{h_1+h_2}\), the harmonic mean.

Exercise 7.3A pad of length \(L=50\) mm runs on oil with \(\mu=0.03\ \mathrm{Pa\,s}\) at \(U=10\ \mathrm{m\,s^{-1}}\) and \(k=2\). What minimum film thickness \(h_2\) carries a load of \(10^4\ \mathrm{N}\) per metre of width?

Show solution
At \(k=2\), \(6(\ln2-\tfrac23)=0.159\). So \(h_2^2=0.159\,\mu UL^2/W=0.159\times0.03\times10\times0.0025/10^4\approx1.19\times10^{-8}\ \mathrm{m^2}\), giving \(h_2\approx0.11\) mm. Films of tens of micrometres are typical in real bearings.
Lesson 08

Dilatancy of granular media

For the general reader. Step on wet sand at the water’s edge and the sand around your foot turns dry and pale for a moment. Reynolds explained why in 1885. Tightly packed grains cannot be sheared without spreading apart, so the sand expands, sucks water in from the surface and looks dry. He called this property dilatancy.

Goals
  • Understand why shearing a dense packing increases its volume.
  • Compute packing fractions of regular sphere packings.
  • Connect dilatancy to modern soil mechanics.

The idea

In his 1885 paper on the dilatancy of media composed of rigid particles in contact, and in lectures at the Royal Institution in 1885–1886, Reynolds pointed out that a dense packing of rigid grains is locked. To deform it, grains must ride up over one another, and the total volume must increase. A loose packing, by contrast, compacts when sheared. He demonstrated the effect with a rubber bag filled with sand or shot and water: when squeezed, the bag became rigid and drew in water, rather than yielding.

Principle of dilatancy (Reynolds, 1885)

A granular mass of rigid particles packed more densely than a critical density can undergo shear deformation only if its volume increases. A mass looser than critical decreases in volume under shear. When the pore space is filled with liquid and drainage is restricted, dilation produces suction in the liquid and contraction produces excess pore pressure.

The modern form of this idea is the critical state theory of soil mechanics developed at Cambridge in the 1950s and 1960s. Arthur Casagrande introduced the critical void ratio in 1936, the density at which a sand shears without volume change. Dilatancy explains why dense sands are strong and loose saturated sands can liquefy in earthquakes. Reynolds later tried to build a theory of the aether on the same principle (Lesson 11).

A geometric model

Consider layers of equal spheres stacked in two ways. In the simple cubic packing each sphere sits directly on the one below, with packing fraction \(\pi/6\approx0.524\). In the densest packing each sphere sits in a hollow of the layer below, with fraction \(\pi/(3\sqrt2)\approx0.740\). Shearing one dense layer over the next forces the upper spheres to climb out of the hollows, raising the layer.

Exercise 8.1Two square-packed layers of spheres of radius \(a\) sit one on the other. In the dense position each upper sphere rests in the hollow between four lower spheres. Find the height of the upper layer’s centres in the dense position and when the upper layer has been sheared so that each sphere sits directly on top of one below. By what fraction does the volume of a pile of such layers increase?

Show solution
In the hollow, the upper centre lies above the centre of a square of side \(2a\), at horizontal distance \(a\sqrt2\) from four lower centres, at distance \(2a\). Its height is \(\sqrt{4a^2-2a^2}=a\sqrt2\). Directly on top, the height is \(2a\). The layer spacing grows by the factor \(2/\sqrt2=\sqrt2\), so the volume rises by about 41 percent, from packing fraction 0.740 to 0.524.

Exercise 8.2A saturated sand with porosity (pore fraction) 0.36 dilates by 2 percent of its total volume under a footstep. If the water cannot drain in time, what happens to the water pressure and to the strength of the sand? Explain the dry halo around the foot.

Show solution
The pore volume must grow from 0.36 to 0.38 of the original volume, but no water is available, so the pore pressure falls below atmospheric. The effective stress between grains, the total stress minus the pore pressure, rises, and the sand becomes stiffer. The suction draws the surface film of water into the pores, so the surface around the foot looks dry. When the foot is lifted, the sand contracts and water returns.
Lesson 09

Reynolds averaging and the turbulent stresses

For the general reader. A turbulent flow is a tangle of eddies that never repeats. Engineers usually care only about its average behaviour. In 1895 Reynolds showed how to write equations for the average flow. They contain new terms, the Reynolds stresses, which describe how the eddies carry momentum. Finding these terms without solving for every eddy is the closure problem of turbulence, still unsolved in general.

Goals
  • Derive the Reynolds-averaged Navier–Stokes equations.
  • Interpret the Reynolds stress tensor and the closure problem.

The decomposition

In On the dynamical theory of incompressible viscous fluids and the determination of the criterion, read in 1894 and published in the Philosophical Transactions in 1895, Reynolds split each velocity component into a mean and a fluctuation, \(u_i=U_i+u_i'\), with \(\overline{u_i'}=0\). He used local space averages; today one uses ensemble or time averages. The averaging operator is linear, commutes with derivatives and satisfies \(\overline{U_iu_j'}=0\).

Reynolds-averaged equations (1895)

Averaging the incompressible Navier–Stokes equations gives \(\partial_iU_i=0\) and

\[ \rho\Big(\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\Big)=-\frac{\partial P}{\partial x_i}+\frac{\partial}{\partial x_j}\Big(\mu\Big(\frac{\partial U_i}{\partial x_j}+\frac{\partial U_j}{\partial x_i}\Big)-\rho\,\overline{u_i'u_j'}\Big). \]

The symmetric tensor \(-\rho\,\overline{u_i'u_j'}\) is the Reynolds stress.

Derivation. Write the convective term in conservation form \(\partial_j(u_iu_j)\), which is allowed because \(\partial_ju_j=0\). Averaging, \(\overline{u_iu_j}=U_iU_j+\overline{u_i'u_j'}\), because the cross terms average to zero. The extra term is moved to the right-hand side, where it acts like a stress. Its trace is \(-2\rho k\), with \(k=\tfrac12\overline{u_i'u_i'}\) the turbulent kinetic energy.

The closure problem

The averaged equations have more unknowns than equations. Equations for \(\overline{u_i'u_j'}\) can be derived, but they contain triple correlations, and so on without end. Joseph Boussinesq had proposed in 1877 to model turbulent stress by an eddy viscosity, \(-\overline{u_i'u_j'}=\nu_t(\partial_jU_i+\partial_iU_j)-\tfrac23k\delta_{ij}\). Ludwig Prandtl’s mixing length (1925) and the \(k\)–\(\varepsilon\) models of the 1970s are ways of choosing \(\nu_t\). Every Reynolds-averaged (RANS) code used in engineering today starts from Reynolds’s 1895 equation.

Exercise 9.1Steady turbulent flow in a channel \(-h<y<h\) is driven by a constant mean pressure gradient, with mean velocity \(U(y)\) along \(x\). Show that the total shear stress is linear: \(\mu\,dU/dy-\rho\,\overline{u'v'}=-\tau_w\,y/h\), where \(\tau_w\) is the wall stress at \(y=-h\).

Show solution
For a statistically steady flow that is homogeneous in \(x\), the \(x\)-momentum equation reduces to \(0=-dP/dx+\frac{d}{dy}\big(\mu U'-\rho\overline{u'v'}\big)\). Integrating, the total stress \(\tau(y)\) is linear in \(y\). By symmetry \(\tau(0)=0\). At the walls the fluctuations vanish by no slip, so \(\tau(-h)=\mu U'(-h)=\tau_w\). Hence \(\tau(y)=-\tau_wy/h\). Near the wall viscosity carries the stress; in the core the Reynolds stress carries almost all of it.

Exercise 9.2Show that averaging does not commute with products: if \(u=U+u'\) and \(v=V+v'\), then \(\overline{uv}=UV+\overline{u'v'}\). Give an example where \(\overline{u'v'}\ne0\).

Show solution
\(\overline{uv}=UV+U\overline{v'}+V\overline{u'}+\overline{u'v'}=UV+\overline{u'v'}\). Take \(u'=\sin\omega t\), \(v'=-\sin\omega t\), averaged over a period: \(\overline{u'v'}=-\tfrac12\). Correlated fluctuations of this sign, with fast fluid moving toward a wall, are exactly what produce turbulent drag.
Lesson 10

The energy method and the critical number

For the general reader. Why does turbulence appear at a Reynolds number near 2000 and not at 20 or 200 000? Reynolds tried to answer this from the equations themselves by tracking the energy of a disturbance: the mean flow feeds it, viscosity drains it. His bookkeeping gives a speed below which every disturbance must die. The gap between that bound and the observed transition is still a central theme of the subject.

Goals
  • Derive the energy equation for a disturbance.
  • State the energy stability theorem and its variational form.
  • Understand why linear and energy theories both fail to predict pipe transition.

The second half of Reynolds’s 1895 paper sought the critical number theoretically, completing the title’s promise of a criterion.

Reynolds’s energy equation

The second half of the paper sought the critical Reynolds number theoretically. Reynolds derived an equation for the kinetic energy of the fluctuations. In the form refined by William Orr in 1907, for a disturbance \(\mathbf u'\) of a steady flow \(\mathbf U\) in a bounded domain with no-slip walls,

\[ \frac{d}{dt}\int\tfrac12|\mathbf u'|^2\,dV=-\int u_i'u_j'\,\frac{\partial U_i}{\partial x_j}\,dV-\nu\int|\nabla\mathbf u'|^2\,dV. \]

The first term transfers energy from the mean flow to the disturbance; the second is viscous dissipation.

Energy stability theorem (Reynolds 1895, Orr 1907)

Let \(\mathbf U\) be a steady solution of the incompressible Navier–Stokes equations in a bounded domain with no-slip walls, made dimensionless with a velocity scale \(U_0\) and length \(L\), so that \(\mathrm{Re}=U_0L/\nu\). Define

\[ \frac{1}{\mathrm{Re}_E}=\max_{\mathbf v}\ \frac{-\int v_iv_j\,\partial_jU_i\,dV}{\int|\nabla\mathbf v|^2\,dV}, \]

the maximum taken over divergence-free fields vanishing on the walls. If \(\mathrm{Re}<\mathrm{Re}_E\), the energy of every disturbance, of any amplitude, decays at least exponentially.

Proof sketch. In dimensionless form the energy equation reads \(\dot E=-\int u_i'u_j'\partial_jU_i\,dV-\mathrm{Re}^{-1}\int|\nabla\mathbf u'|^2dV\le-\big(\mathrm{Re}^{-1}-\mathrm{Re}_E^{-1}\big)\int|\nabla\mathbf u'|^2dV\). The Poincaré inequality \(\int|\nabla\mathbf u'|^2\ge\lambda_1\int|\mathbf u'|^2\) then gives \(\dot E\le-cE\) with \(c>0\). The maximization is a linear eigenvalue problem, the Euler–Lagrange equation of the ratio.

This gives a rigorous lower bound for the transition Reynolds number. For pipe flow, based on the mean speed and diameter, \(\mathrm{Re}_E\approx81\), far below the observed value of about 2000. Between the two, disturbances can grow for a while before they decay, a phenomenon now called transient growth, and turbulence can sustain itself only above about 2000. Orr corrected and sharpened Reynolds’s analysis in 1907, and the energy method in this form became a standard tool of stability theory.

Linear theory

A second approach linearizes the equations about the base flow and looks for exponentially growing modes. William Orr and Arnold Sommerfeld formulated the resulting Orr–Sommerfeld equation in 1907–1908. For plane Poiseuille flow it predicts instability above \(\mathrm{Re}\approx5772\) (based on centreline speed and half-width), a value computed by Steven Orszag in 1971. For pipe flow linear theory finds no instability at all. Pipe turbulence is therefore a finite-amplitude phenomenon, exactly as Reynolds’s observations of calm and disturbed inlets suggested.

Exercise 10.1Derive the energy equation above for a disturbance of a steady solution \(\mathbf U\) of the Navier–Stokes equations. (Take the dot product of the disturbance equation with \(\mathbf u'\) and integrate.)

Show solution
Subtracting the equations for \(\mathbf U\) and \(\mathbf U+\mathbf u'\): \(\partial_tu_i'+U_j\partial_ju_i'+u_j'\partial_jU_i+u_j'\partial_ju_i'=-\partial_ip'/\rho+\nu\nabla^2u_i'\). Multiply by \(u_i'\) and integrate. Terms of the form \(\int a_j\partial_j(\tfrac12|\mathbf u'|^2)\) with \(\nabla\!\cdot\mathbf a=0\) vanish after integration by parts, as does the pressure term, because \(\mathbf u'=0\) on the walls. The viscous term becomes \(-\nu\int|\nabla\mathbf u'|^2\). What remains is \(-\int u_i'u_j'\partial_jU_i\). The nonlinear term drops out: it redistributes energy among disturbances but never creates it.

Exercise 10.2For plane Couette flow \(U=(Sy,0,0)\), show that the production term is \(-S\int u'v'\,dV\). Which sign of correlation feeds the disturbance?

Show solution
Only \(\partial_yU_x=S\) is non-zero, so \(-\int u_i'u_j'\partial_jU_i=-S\int u'v'\). For \(S>0\) energy grows when \(\overline{u'v'}<0\): fluid moving up (\(v'>0\)) is slower than its surroundings (\(u'<0\)), carrying low momentum upward. This is the same correlation that produces a positive turbulent shear stress \(-\rho\overline{u'v'}\).
Part III · Life and legacy
Lesson 11

Heat, Joule and the sub-mechanics of the universe, 1888–1903

By the late 1880s Reynolds was one of the leading engineering scientists in Britain. The Royal Society awarded him a Royal Medal in 1888, and he became president of the Manchester Literary and Philosophical Society. The new Whitworth laboratories at Owens College gave him space for large experiments.

The largest was a new determination of the mechanical equivalent of heat. Manchester was James Prescott Joule’s city, and after Joule’s death in 1889 Reynolds wrote a long memoir of him for the Lit & Phil, published in 1892. He then used a large steam engine and a hydraulic brake of his own design to measure the work needed to heat water from freezing to boiling point. The result, a mean specific heat of water expressed in mechanical units, was reported to the Royal Society in the late 1890s and is regarded as one of the classical determinations of a physical constant. He gave the Royal Society’s Bakerian Lecture in 1897. He also wrote on the refraction of sound by wind and temperature gradients, on condensation of steam, and on group velocity, where Lamb credited him with showing that the group velocity is also the speed at which wave energy travels.

Collected papers and the granular aether

Between 1900 and 1903 Cambridge University Press published his Papers on Mechanical and Physical Subjects in three volumes, gathering about seventy papers. The third volume was a single work, The Sub-Mechanics of the Universe (1903), developed from his 1902 Rede Lecture at Cambridge. In it Reynolds proposed that space is filled with a dense medium of tiny rigid grains, a dilatant aether like his rubber bag of shot, and tried to derive light, gravitation and electricity from its mechanics. In passing he stated the transport theorem for integrals over moving regions that now carries his name:

\[ \frac{d}{dt}\int_{\Omega(t)}f\,dV=\int_{\Omega(t)}\frac{\partial f}{\partial t}\,dV+\int_{\partial\Omega(t)}f\,\mathbf v_b\cdot\mathbf n\,dA. \]

The book was hard to read and appeared just as relativity and quantum theory were making aether models obsolete. George Bryan reviewed it sympathetically in Nature, but it had little influence.

Exercise 11.1Check the transport theorem in one dimension, where it is the Leibniz rule: \(\frac{d}{dt}\int_{a(t)}^{b(t)}f(x,t)\,dx\). Then use it with \(\mathbf v_b=\mathbf u\) and mass conservation to show \(\frac{d}{dt}\int_{\Omega(t)}\rho\,\phi\,dV=\int_{\Omega(t)}\rho\frac{D\phi}{Dt}\,dV\) for a material region.

Show solution
In one dimension the formula reads \(\int_a^b\partial_tf\,dx+f(b,t)\dot b-f(a,t)\dot a\), the Leibniz rule; the outward normal is \(+1\) at \(b\) and \(-1\) at \(a\). For a material region, apply the theorem to \(f=\rho\phi\) and use the divergence theorem: \(\int[\partial_t(\rho\phi)+\nabla\!\cdot(\rho\phi\mathbf u)]\,dV\). Expanding and using \(\partial_t\rho+\nabla\!\cdot(\rho\mathbf u)=0\) leaves \(\int\rho(\partial_t\phi+\mathbf u\cdot\nabla\phi)\,dV\).
Lesson 12

Decline, death and legacy, 1905–1912

In the early 1900s Reynolds’s health began to fail, and his mental powers declined with it. He retired from his chair in 1905, after 37 years. He received the Dalton Medal of the Manchester Lit & Phil in 1903. He lived his last years in Somerset and died of influenza at Watchet on 21 February 1912, aged 69. He is buried at St Decuman’s church in Watchet. Horace Lamb wrote his obituary for the Royal Society, and Arthur Gibson later wrote a book-length account of his work in hydraulics.

The school

Reynolds trained a generation of engineers at Manchester. Thomas Stanton studied in the Whitworth laboratory under him from 1888 and stayed as demonstrator until 1896. At the National Physical Laboratory, Stanton and J. R. Pannell confirmed in 1914 that pipe friction for air and water collapses onto a single curve of friction factor against Reynolds number, the decisive test of dynamical similarity. J. J. Thomson, who passed through Owens College in the 1870s, went on to the Cavendish chair and the discovery of the electron.

Legacy

  • Fluid mechanics: the Reynolds number organizes all of fluid dynamics, from bacteria swimming at \(\mathrm{Re}\sim10^{-4}\) to atmospheric flows at \(10^{9}\) and more. The question of how laminar pipe flow becomes turbulent was studied intensively for more than a century after 1883.
  • Turbulence: Reynolds averaging and the Reynolds stresses are the starting point of turbulence modelling and computational fluid dynamics.
  • Tribology: the Reynolds equation is the foundation of hydrodynamic bearing design; Duncan Dowson counted him among the founding figures of the subject.
  • Heat transfer and granular physics: the Reynolds analogy and the concept of dilatancy remain standard.
  • Commemoration: the University of Manchester held a centenary symposium in 1968 for the founding of his chair and has held an annual Reynolds Day on fluid mechanics since 2001.

Final projectChoose one of Reynolds’s papers from Part II. Write a two-page account for a strong undergraduate that explains the engineering problem behind it, derives its central equation, works one numerical example, and describes how the result is used today.

Reference

Connections

People linked to Osborne Reynolds across this collection. A name without a link has no page here yet. See Osborne Reynolds in the network of all 337 lives →

Students

Influenced by

Influenced

Colleagues and collaborators

Also mentioned in this life

Fields

Reference

Sources and further reading

Biography

  1. J. J. O’Connor and E. F. Robertson, Osborne Reynolds, MacTutor History of Mathematics.
  2. Osborne Reynolds, Wikipedia.
  3. J. Allen, Reynolds, Osborne, Dictionary of Scientific Biography; online at Encyclopedia.com.
  4. H. Lamb, Osborne Reynolds (obituary), Proc. Roy. Soc. A 88 (1913).
  5. A. H. Gibson, Osborne Reynolds and His Work in Hydraulics and Hydrodynamics, Longmans, 1946.
  6. D. M. McDowell and J. D. Jackson (eds.), Osborne Reynolds and Engineering Science Today, Manchester University Press, 1970.
  7. J. J. Thomson, Recollections and Reflections, Bell, 1936.
  8. APS News, March 1883: Osborne Reynolds and the onset of turbulence.

Physics

  1. O. Reynolds, On the extent and action of the heating surface of steam boilers, Proc. Manchester Lit. Phil. Soc. 14 (1874–75).
  2. O. Reynolds, On certain dimensional properties of matter in the gaseous state, Phil. Trans. R. Soc. 170 (1879).
  3. O. Reynolds, An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, Phil. Trans. R. Soc. 174 (1883) 935–982.
  4. O. Reynolds, On the dilatancy of media composed of rigid particles in contact, Phil. Mag. 20 (1885).
  5. O. Reynolds, On the theory of lubrication and its application to Mr. Beauchamp Tower’s experiments, Phil. Trans. R. Soc. 177 (1886) 157–234.
  6. O. Reynolds, On the dynamical theory of incompressible viscous fluids and the determination of the criterion, Phil. Trans. R. Soc. A 186 (1895) 123–164.
  7. O. Reynolds, Papers on Mechanical and Physical Subjects, 3 vols., Cambridge University Press, 1900–1903.
  8. N. Rott, Note on the history of the Reynolds number, Annu. Rev. Fluid Mech. 22 (1990) 1–11.
  9. K. Avila et al., The onset of turbulence in pipe flow, Science 333 (2011) 192–196.
  10. S. B. Pope, Turbulent Flows, Cambridge University Press, 2000.
  11. D. Dowson, History of Tribology, Longman, 1979.
  12. Reynolds number and Reynolds transport theorem, Wikipedia.