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.
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.
Timeline
| 1842 | Born 23 August in Belfast, son of the Rev. Osborne Reynolds, headmaster of the Belfast Collegiate School. The family soon moves to Dedham, Essex. |
| 1861 | Apprenticed to the engineer Edward Hayes at Stony Stratford. |
| 1867 | Seventh wrangler in the Mathematical Tripos; fellow of Queens’ College, Cambridge. A year as a civil engineer with Lawson and Mansergh in London. |
| 1868 | Appointed to the new chair of engineering at Owens College, Manchester, aged 25. Marries Charlotte Jemima Chadwick. |
| 1869 | Elected to the Manchester Literary and Philosophical Society. His first wife dies after the birth of their son. |
| 1874 | Paper on the heating surface of steam boilers: the Reynolds analogy. Secretary of the Lit & Phil, 1874–1884. |
| 1877 | Fellow of the Royal Society. |
| 1879 | Long memoir on thermal transpiration and the dynamical theory of gases in the Philosophical Transactions. |
| 1883 | March: submits the paper on direct and sinuous flow in pipes, the origin of the Reynolds number. |
| 1884 | Honorary LLD, Glasgow. Candidate for the Cavendish chair at Cambridge, which goes to J. J. Thomson. |
| 1885 | Names and explains dilatancy in granular media. |
| 1886 | On the theory of lubrication, explaining Beauchamp Tower’s bearing experiments. |
| 1888 | Royal Medal of the Royal Society. President of the Lit & Phil. |
| 1892 | Publishes his memoir of James Prescott Joule. |
| 1895 | Paper on the dynamical theory of incompressible viscous fluids: Reynolds averaging, the Reynolds stresses and the energy method. |
| 1897 | Bakerian Lecture of the Royal Society. |
| 1900–03 | Papers on Mechanical and Physical Subjects, three volumes. |
| 1902–03 | Rede Lecture at Cambridge; The Sub-Mechanics of the Universe. Dalton Medal, 1903. |
| 1905 | Retires because of failing health. |
| 1912 | Dies of influenza on 21 February at Watchet, Somerset; buried at St Decuman’s church. |
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?
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.
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.
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.
- 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\).
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
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
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.
- 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}. \]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.
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
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
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
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.
- 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.
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
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
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.
- 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}. \]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
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
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
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.
- 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.
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
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
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.
- 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\).
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
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
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.
- 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.
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
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
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
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.
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
- J. J. Thomson 1870-1876
- Thomas Stanton 1888-1896
Influenced by
- George Gabriel Stokes
- James Prescott Joule
- Joseph Boussinesq
Influenced
- Arnold Sommerfeld
- Ludwig Prandtl
- William McFadden Orr
Colleagues and collaborators
- James Clerk Maxwell 1879
- James Prescott Joule
- John William Strutt, Lord Rayleigh 1863-1867
- William Crookes 1879
- Horace Lamb 1885-1905
Also mentioned in this life
- Henry Cavendish
- James Chadwick
- John Dalton
- Leonhard Euler
- Joseph-Louis Lagrange
- Gottfried Wilhelm Leibniz
- Claude-Louis Navier
- Henri Poincaré
Fields
- fluid mechanics 6 lives
- engineering 6 lives
- turbulence 5 lives
Sources and further reading
Biography
- J. J. O’Connor and E. F. Robertson, Osborne Reynolds, MacTutor History of Mathematics.
- Osborne Reynolds, Wikipedia.
- J. Allen, Reynolds, Osborne, Dictionary of Scientific Biography; online at Encyclopedia.com.
- H. Lamb, Osborne Reynolds (obituary), Proc. Roy. Soc. A 88 (1913).
- A. H. Gibson, Osborne Reynolds and His Work in Hydraulics and Hydrodynamics, Longmans, 1946.
- D. M. McDowell and J. D. Jackson (eds.), Osborne Reynolds and Engineering Science Today, Manchester University Press, 1970.
- J. J. Thomson, Recollections and Reflections, Bell, 1936.
- APS News, March 1883: Osborne Reynolds and the onset of turbulence.
Physics
- O. Reynolds, On the extent and action of the heating surface of steam boilers, Proc. Manchester Lit. Phil. Soc. 14 (1874–75).
- O. Reynolds, On certain dimensional properties of matter in the gaseous state, Phil. Trans. R. Soc. 170 (1879).
- 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.
- O. Reynolds, On the dilatancy of media composed of rigid particles in contact, Phil. Mag. 20 (1885).
- O. Reynolds, On the theory of lubrication and its application to Mr. Beauchamp Tower’s experiments, Phil. Trans. R. Soc. 177 (1886) 157–234.
- 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.
- O. Reynolds, Papers on Mechanical and Physical Subjects, 3 vols., Cambridge University Press, 1900–1903.
- N. Rott, Note on the history of the Reynolds number, Annu. Rev. Fluid Mech. 22 (1990) 1–11.
- K. Avila et al., The onset of turbulence in pipe flow, Science 333 (2011) 192–196.
- S. B. Pope, Turbulent Flows, Cambridge University Press, 2000.
- D. Dowson, History of Tribology, Longman, 1979.
- Reynolds number and Reynolds transport theorem, Wikipedia.