Reynolds number

The Reynolds number (Re) is the ratio between the inertial forces in a fluid and the viscous forces. A fluid in motion tends to behave as sheets or layers of infinitely small thicknesses (smaller than the wavelength of light) sliding relative to each other. The viscosity of a fluid is the resistance to flow, which in turn translates as resistance to shear between the layers. The inertia (momentum) of the fluid is a dynamic function of its mass and speed, or resistance to a change in motion, which works to create shear between the fluid layers. This shear is created as a fluid encounters viscous resistance from an outside obstacle, such as the inner wall of a pipe, where flow is slowed by drag against the surface while remaining mostly unimpeded near the center of the pipe. For any given fluid, as flow rate increases, at a certain point the inertial forces begin to overcome the viscous forces; the smoothly-sliding layers of the fluid begin to roll past each other, and rough, turbulent flow results. Reynolds number is used to predict both where and at what scale in the flow this transition will happen for any given fluid in any given flow situation.

Reynolds number is an important dimensionless quantity in fluid mechanics used to help predict flow patterns in different fluid flow-situations. At low Reynolds numbers, flows tend to be dominated by laminar (sheet-like) flow, while at high Reynolds numbers turbulence results from differences in the fluid's speed and direction, which may sometimes intersect or even move counter to the overall direction of the flow (eddy currents). These eddy currents begin to churn the flow, using up energy in the process, which for liquids increases the chances of cavitation. The Reynolds number has wide applications, ranging from liquid flow in a pipe to the passage of air over an aircraft wing. It is used to predict the transition from laminar to turbulent flow, and is used in the scaling of similar but different-sized flow situations, such as between an aircraft model in a wind tunnel and the full size version. The predictions of the onset of turbulence and the ability to calculate scaling effects can be used to help predict fluid behaviour on any scale, ranging from microscopic scales where Reynolds numbers tend to be extremely low and viscosity dominates, to very large scales such as in local or global air or water movement and the associated meteorological and climatological effects, to astronomical scales such as accretion flows around proto-stars, galaxies, and black holes where Reynolds numbers are extremely high and turbulent flow dominates.[2]

The concept was introduced by George Stokes in 1851,[3] but the Reynolds number was named by Arnold Sommerfeld in 1908[4] after Osborne Reynolds (1842–1912), who popularized its use in 1883.[5][6]

Definition

The Reynolds number is the ratio of inertial forces to viscous forces within a fluid which is subjected to relative internal movement due to different fluid velocities, which is known as a boundary layer in the case of a bounding surface such as the interior of a pipe. A similar effect is created by the introduction of a stream of high-velocity fluid into a low-velocity fluid, such as the hot gases emitted from a flame in air. This relative movement generates fluid friction, which is a factor in developing turbulent flow. Counteracting this effect is the viscosity of the fluid, which tends to inhibit turbulence. The Reynolds number quantifies the relative importance of these two types of forces for given flow conditions, and is a guide to when turbulent flow will occur in a particular situation.[7]

This ability to predict the onset of turbulent flow is an important design tool for equipment such as piping systems or aircraft wings, but the Reynolds number is also used in scaling of fluid dynamics problems, and is used to determine dynamic similitude between two different cases of fluid flow, such as between a model aircraft, and its full-size version. Such scaling is not linear and the application of Reynolds numbers to both situations allows scaling factors to be developed.

With respect to laminar and turbulent flow regimes:

• laminar flow occurs at low Reynolds numbers, where viscous forces are dominant, and is characterized by smooth, constant fluid motion;
• turbulent flow occurs at high Reynolds numbers and is dominated by inertial forces, which tend to produce chaotic eddies, vortices and other flow instabilities.[8]

The Reynolds number is defined as[4]

${\displaystyle \mathrm {Re} ={\frac {\rho uL}{\mu }}={\frac {uL}{\nu }}}$

where:

The Reynolds number can be defined for several different situations where a fluid is in relative motion to a surface.[n 1] These definitions generally include the fluid properties of density and viscosity, plus a velocity and a characteristic length or characteristic dimension (L in the above equation). This dimension is a matter of convention – for example radius and diameter are equally valid to describe spheres or circles, but one is chosen by convention. For aircraft or ships, the length or width can be used. For flow in a pipe, or for a sphere moving in a fluid, the internal diameter is generally used today. Other shapes such as rectangular pipes or non-spherical objects have an equivalent diameter defined. For fluids of variable density such as compressible gases or fluids of variable viscosity such as non-Newtonian fluids, special rules apply. The velocity may also be a matter of convention in some circumstances, notably stirred vessels.

In practice, matching the Reynolds number is not on its own sufficient to guarantee similitude. Fluid flow is generally chaotic, and very small changes to shape and surface roughness of bounding surfaces can result in very different flows. Nevertheless, Reynolds numbers are a very important guide and are widely used.

History

Osborne Reynolds famously studied the conditions in which the flow of fluid in pipes transitioned from laminar flow to turbulent flow. In his 1883 paper Reynolds described the transition from laminar to turbulent flow in a classic experiment in which he examined the behaviour of water flow under different flow velocities using a small stream of dyed water introduced into the centre of clear water flow in a larger pipe.

The larger pipe was glass so the behaviour of the layer of the dyed stream could be observed, and at the end of this pipe there was a flow control valve used to vary the water velocity inside the tube. When the velocity was low, the dyed layer remained distinct through the entire length of the large tube. When the velocity was increased, the layer broke up at a given point and diffused throughout the fluid's cross-section. The point at which this happened was the transition point from laminar to turbulent flow.

From these experiments came the dimensionless Reynolds number for dynamic similarity—the ratio of inertial forces to viscous forces. Reynolds also proposed what is now known as the Reynolds-averaging of turbulent flows, where quantities such as velocity are expressed as the sum of mean and fluctuating components. Such averaging allows for 'bulk' description of turbulent flow, for example using the Reynolds-averaged Navier–Stokes equations.

Flow in a pipe

For flow in a pipe or tube, the Reynolds number is generally defined as[9]

${\displaystyle \mathrm {Re} ={\frac {\rho uD_{\text{H}}}{\mu }}={\frac {uD_{\text{H}}}{\nu }}={\frac {QD_{\text{H}}}{\nu A}},}$

where

DH is the hydraulic diameter of the pipe (the inside diameter if the pipe is circular) (m),
Q is the volumetric flow rate (m3/s),
A is the pipe's cross-sectional area (m2),
u is the mean velocity of the fluid (m/s),
μ (mu) is the dynamic viscosity of the fluid (Pa·s = N·s/m2 = kg/(m·s)),
ν (nu) is the kinematic viscosity (ν = μ/ρ) (m2/s),
ρ is the density of the fluid (kg/m3).

For shapes such as squares, rectangular or annular ducts where the height and width are comparable, the characteristic dimension for internal-flow situations is taken to be the hydraulic diameter, DH, defined as

${\displaystyle D_{\text{H}}={\frac {4A}{P}},}$

where A is the cross-sectional area, and P is the wetted perimeter. The wetted perimeter for a channel is the total perimeter of all channel walls that are in contact with the flow.[10] This means that the length of the channel exposed to air is not included in the wetted perimeter.

For a circular pipe, the hydraulic diameter is exactly equal to the inside pipe diameter:

${\displaystyle D_{\text{H}}=D.}$

For an annular duct, such as the outer channel in a tube-in-tube heat exchanger, the hydraulic diameter can be shown algebraically to reduce to

${\displaystyle D_{\text{H,annulus}}=D_{\text{o}}-D_{\text{i}},}$

where

Do is the inside diameter of the outer pipe,
Di is the outside diameter of the inner pipe.

For calculation involving flow in non-circular ducts, the hydraulic diameter can be substituted for the diameter of a circular duct, with reasonable accuracy, if the aspect ratio AR of the duct cross-section remains in the range 1/4 < AR < 4.[11]

Laminar–turbulent transition

In boundary layer flow over a flat plate, experiments confirm that, after a certain length of flow, a laminar boundary layer will become unstable and turbulent. This instability occurs across different scales and with different fluids, usually when Rex5×105,[12] where x is the distance from the leading edge of the flat plate, and the flow velocity is the freestream velocity of the fluid outside the boundary layer.

For flow in a pipe of diameter D, experimental observations show that for "fully developed" flow,[n 2] laminar flow occurs when ReD < 2300 and turbulent flow occurs when ReD > 2900.[13][14] At the lower end of this range, a continuous turbulent-flow will form, but only at a very long distance from the inlet of the pipe. The flow in between will begin to transition from laminar to turbulent and then back to laminar at irregular intervals, called intermittent flow. This is due to the different speeds and conditions of the fluid in different areas of the pipe's cross-section, depending on other factors such as pipe roughness and flow uniformity. Laminar flow tends to dominate in the fast-moving center of the pipe while slower-moving turbulent flow dominates near the wall. As the Reynolds number increases, the continuous turbulent-flow moves closer to the inlet and the intermittency in between increases, until the flow becomes fully turbulent at ReD > 2900.[13] This result is generalized to non-circular channels using the hydraulic diameter, allowing a transition Reynolds number to be calculated for other shapes of channel.[13]

These transition Reynolds numbers are also called critical Reynolds numbers, and were studied by Osborne Reynolds around 1895.[6] The critical Reynolds number is different for every geometry.[15]

Flow in a wide duct

For a fluid moving between two plane parallel surfaces—where the width is much greater than the space between the plates—then the characteristic dimension is equal to twice the distance between the plates. This is consistent with the annular duct and rectangular duct cases above taken to a limiting aspect ratio.

Flow in an open channel

For flow of liquid with a free surface, the hydraulic radius must be determined. This is the cross-sectional area of the channel divided by the wetted perimeter. For a semi-circular channel, it is quarter the diameter (in case of full pipe flow). For a rectangular channel, the hydraulic radius is the cross-sectional area divided by the wetted perimeter. Some texts then use a characteristic dimension that is four times the hydraulic radius, chosen because it gives the same value of Re for the onset of turbulence as in pipe flow,[16] while others use the hydraulic radius as the characteristic length-scale with consequently different values of Re for transition and turbulent flow.

Flow around airfoils

Reynolds numbers are used in airfoil design to (among other things) manage "scale effect" when computing/comparing characteristics (a tiny wing, scaled to be huge, will perform differently).[17] Fluid dynamicists define the chord Reynolds number R like this: R = Vc/ν, where V is the flight speed, c is the chord length, and ν is the kinematic viscosity of the fluid in which the airfoil operates, which is 1.460×10−5 m2/s for the atmosphere at sea level.[18] In some special studies a characteristic length other than chord may be used; rare is the "span Reynolds number", which is not to be confused with spanwise stations on a wing, where chord is still used.[19]

Object in a fluid

The Reynolds number for an object moving in a fluid, called the particle Reynolds number and often denoted Rep, characterizes the nature of the surrounding flow and its fall velocity.

In viscous fluids

Where the viscosity is naturally high, such as polymer solutions and polymer melts, flow is normally laminar. The Reynolds number is very small and Stokes' law can be used to measure the viscosity of the fluid. Spheres are allowed to fall through the fluid and they reach the terminal velocity quickly, from which the viscosity can be determined.

The laminar flow of polymer solutions is exploited by animals such as fish and dolphins, who exude viscous solutions from their skin to aid flow over their bodies while swimming. It has been used in yacht racing by owners who want to gain a speed advantage by pumping a polymer solution such as low molecular weight polyoxyethylene in water, over the wetted surface of the hull.

It is, however, a problem for mixing of polymers, because turbulence is needed to distribute fine filler (for example) through the material. Inventions such as the "cavity transfer mixer" have been developed to produce multiple folds into a moving melt so as to improve mixing efficiency. The device can be fitted onto extruders to aid mixing.

Sphere in a fluid

For a sphere in a fluid, the characteristic length-scale is the diameter of the sphere and the characteristic velocity is that of the sphere relative to the fluid some distance away from the sphere, such that the motion of the sphere does not disturb that reference parcel of fluid. The density and viscosity are those belonging to the fluid.[20] Note that purely laminar flow only exists up to Re = 10 under this definition.

Under the condition of low Re, the relationship between force and speed of motion is given by Stokes' law.[21]

Rectangular object in a fluid

The equation for a rectangular object is identical to that of a sphere, with the object being approximated as an ellipsoid and the axis of length being chosen as the characteristic length scale. Such considerations are important in natural streams, for example, where there are few perfectly spherical grains. For grains in which measurement of each axis is impractical, sieve diameters are used instead as the characteristic particle length-scale. Both approximations alter the values of the critical Reynolds number.

Fall velocity

The particle Reynolds number is important in determining the fall velocity of a particle. When the particle Reynolds number indicates laminar flow, Stokes' law can be used to calculate its fall velocity. When the particle Reynolds number indicates turbulent flow, a turbulent drag law must be constructed to model the appropriate settling velocity.

Packed bed

For fluid flow through a bed, of approximately spherical particles of diameter D in contact, if the voidage is ε and the superficial velocity is vs, the Reynolds number can be defined as[22]

${\displaystyle \mathrm {Re} ={\frac {\rho v_{\text{s}}D}{\mu }},}$

or

${\displaystyle \mathrm {Re} ={\frac {\rho v_{\text{s}}D}{\mu \varepsilon }},}$

or

${\displaystyle \mathrm {Re} ={\frac {\rho v_{\text{s}}D}{\mu (1-\varepsilon )}}.}$

The choice of equation depends on the system involved: the first is successful in correlating the data for various types of packed and fluidized beds, the second Reynolds number suits for the liquid-phase data, while the third was found successful in correlating the fluidized bed data, being first introduced for liquid fluidized bed system.[22]

Laminar conditions apply up to Re = 10, fully turbulent from Re = 2000.[20]

Stirred vessel

In a cylindrical vessel stirred by a central rotating paddle, turbine or propeller, the characteristic dimension is the radius of the agitator r. The velocity V is Nr where N is the rotational speed in rad per second. Then the Reynolds number is:

${\displaystyle \mathrm {Re} ={\frac {\rho Nr^{2}}{\mu }}={\frac {\rho Vr}{\mu }}.}$

The system is fully turbulent for values of Re above 10000.[23]

Pipe friction

Pressure drops[24] seen for fully developed flow of fluids through pipes can be predicted using the Moody diagram which plots the Darcy–Weisbach friction factor f against Reynolds number Re and relative roughness ε/D. The diagram clearly shows the laminar, transition, and turbulent flow regimes as Reynolds number increases. The nature of pipe flow is strongly dependent on whether the flow is laminar or turbulent.

Similarity of flows

In order for two flows to be similar, they must have the same geometry and equal Reynolds and Euler numbers. When comparing fluid behavior at corresponding points in a model and a full-scale flow, the following holds:

{\displaystyle {\begin{aligned}\mathrm {Re} _{\text{m}}&=\mathrm {Re} ,\\\mathrm {Eu} _{\text{m}}&=\mathrm {Eu} ,\end{aligned}}}

where ${\displaystyle \mathrm {Re} _{\text{m}}}$ is the Reynolds number for the model, and ${\displaystyle \mathrm {Re} }$ is full-scale Reynolds number, and similarly for the Euler numbers.

The model numbers and design numbers should be in the same proportion, hence

${\displaystyle {\frac {p_{\text{m}}}{\rho _{\text{m}}v_{\text{m}}^{2}}}={\frac {p}{\rho v^{2}}}.}$

This allows engineers to perform experiments with reduced scale models in water channels or wind tunnels and correlate the data to the actual flows, saving on costs during experimentation and on lab time. Note that true dynamic similitude may require matching other dimensionless numbers as well, such as the Mach number used in compressible flows, or the Froude number that governs open-channel flows. Some flows involve more dimensionless parameters than can be practically satisfied with the available apparatus and fluids, so one is forced to decide which parameters are most important. For experimental flow modeling to be useful, it requires a fair amount of experience and judgment of the engineer.

An example where the mere Reynolds number is not sufficient for similarity of flows (or even the flow regime – laminar or turbulent) are bounded flows, i.e. flows that are restricted by walls or other boundaries. A classical example of this is the Taylor–Couette flow, where the dimensionless ratio of radii of bounding cylinders is also important, and many technical applications where these distinctions play an important role.[25][26] Principles of these restrictions were developed by Maurice Marie Alfred Couette and Geoffrey Ingram Taylor and developed further by Floris Takens and David Ruelle.

Typical values of Reynolds number[27][28]

Smallest scales of turbulent motion

In a turbulent flow, there is a range of scales of the time-varying fluid motion. The size of the largest scales of fluid motion (sometimes called eddies) are set by the overall geometry of the flow. For instance, in an industrial smoke stack, the largest scales of fluid motion are as big as the diameter of the stack itself. The size of the smallest scales is set by the Reynolds number. As the Reynolds number increases, smaller and smaller scales of the flow are visible. In a smoke stack, the smoke may appear to have many very small velocity perturbations or eddies, in addition to large bulky eddies. In this sense, the Reynolds number is an indicator of the range of scales in the flow. The higher the Reynolds number, the greater the range of scales. The largest eddies will always be the same size; the smallest eddies are determined by the Reynolds number.

What is the explanation for this phenomenon? A large Reynolds number indicates that viscous forces are not important at large scales of the flow. With a strong predominance of inertial forces over viscous forces, the largest scales of fluid motion are undamped—there is not enough viscosity to dissipate their motions. The kinetic energy must "cascade" from these large scales to progressively smaller scales until a level is reached for which the scale is small enough for viscosity to become important (that is, viscous forces become of the order of inertial ones). It is at these small scales where the dissipation of energy by viscous action finally takes place. The Reynolds number indicates at what scale this viscous dissipation occurs.

In physiology

Poiseuille's law on blood circulation in the body is dependent on laminar flow. In turbulent flow the flow rate is proportional to the square root of the pressure gradient, as opposed to its direct proportionality to pressure gradient in laminar flow.

Using the definition of the Reynolds number we can see that a large diameter with rapid flow, where the density of the blood is high, tends towards turbulence. Rapid changes in vessel diameter may lead to turbulent flow, for instance when a narrower vessel widens to a larger one. Furthermore, a bulge of atheroma may be the cause of turbulent flow, where audible turbulence may be detected with a stethoscope.

Complex systems

Reynolds number interpretation has been extended into the area of arbitrary complex systems. Such as financial flows,[29] nonlinear networks,[30] etc. In the latter case an artificial viscosity is reduced to nonlinear mechanism of energy distribution in complex network media. Reynolds number then represents a basic control parameter which expresses a balance between injected and dissipated energy flows for open boundary system. It has been shown [30] that Reynolds critical regime separates two types of phase space motion: accelerator (attractor) and decelerator. High Reynolds number leads to a chaotic regime transition only in frame of strange attractor model.

Derivation

The Reynolds number can be obtained when one uses the nondimensional form of the incompressible Navier–Stokes equations for a newtonian fluid expressed in terms of the Lagrangian derivative:

${\displaystyle \rho {\frac {D\mathbf {v} }{Dt}}=-\nabla p+\mu \nabla ^{2}\mathbf {v} +\rho \mathbf {f} .}$

Each term in the above equation has the units of a "body force" (force per unit volume) with the same dimensions of a density times an acceleration. Each term is thus dependent on the exact measurements of a flow. When one renders the equation nondimensional, that is when we multiply it by a factor with inverse units of the base equation, we obtain a form that does not depend directly on the physical sizes. One possible way to obtain a nondimensional equation is to multiply the whole equation by the factor

${\displaystyle {\frac {L}{\rho V^{2}}},}$

where

V is the mean velocity, v or v, relative to the fluid (m/s),
L is the characteristic length (m),
ρ is the fluid density (kg/m3).

If we now set

${\displaystyle \mathbf {v} '={\frac {\mathbf {v} }{V}},\quad p'=p{\frac {1}{\rho V^{2}}},\quad \mathbf {f} '=\mathbf {f} {\frac {L}{V^{2}}},\quad {\frac {\partial }{\partial t'}}={\frac {L}{V}}{\frac {\partial }{\partial t}},\quad \nabla '=L\nabla ,}$

we can rewrite the Navier–Stokes equation without dimensions:

${\displaystyle {\frac {D\mathbf {v} '}{Dt'}}=-\nabla 'p'+{\frac {\mu }{\rho LV}}\nabla '^{2}\mathbf {v} '+\mathbf {f} ',}$

where the term μ/ρLV = 1/Re.

Finally, dropping the primes for ease of reading:

${\displaystyle {\frac {D\mathbf {v} }{Dt}}=-\nabla p+{\frac {1}{\mathrm {Re} }}\nabla ^{2}\mathbf {v} +\mathbf {f} .}$

This is why mathematically all Newtonian, incompressible flows with the same Reynolds number are comparable. Notice also that in the above equation, the viscous terms vanish for Re → ∞. Thus flows with high Reynolds numbers are approximately inviscid in the free stream.

Relationship to other dimensionless parameters

There are many dimensionless numbers in fluid mechanics. The Reynolds number measures the ratio of advection and diffusion effects on structures in the velocity field, and is therefore closely related to Péclet numbers, which measure the ratio of these effects on other fields carried by the flow, for example temperature and magnetic fields. Replacement of the kinematic viscosity ν = μ/ρ in Re by the thermal or magnetic diffusivity results in respectively the thermal Péclet number and the magnetic Reynolds number. These are therefore related to Re by products with ratios of diffusivities, namely the Prandtl number and magnetic Prandtl number.

References

Footnotes

1. The definition of the Reynolds number is not to be confused with the Reynolds equation or lubrication equation.
2. Full development of the flow occurs as the flow enters the pipe, the boundary layer thickens and then stabilizes after several diameters distance into the pipe.

Citations

1. Tansley & Marshall 2001, pp. 3274–3283.
3. Stokes 1851, pp. 8–106.
4. Sommerfeld 1908, pp. 116–124.
5. Reynolds 1883, pp. 935–982.
6. Rott 1990, pp. 1–11.
7. Hall, Nancy (5 May 2015). "Boundary Layer". Glenn Research Center. Retrieved 17 September 2019.
8. "Reynolds Number". Engineeringtoolbox.com. 2003.
9. Fox, McDonald & Pritchard 2004, p. 348.
10. Schlichting & Gersten 2017, pp. 416–419.
11. Holman 2002, p. 207.
12. Potter, Ramadan & Wiggert 2012, p. 105.
13. Lissaman 1983, pp. 223–239.
14. "International Standard Atmosphere". eng.cam.ac.uk. Retrieved 17 September 2019.
15. Ehrenstein & Eloy 2013, pp. 321-346.
16. Rhodes 1989, p. 29.
17. Dusenbery 2009, p. 49.
18. Dwivedi 1977, pp. 157-165.
19. "Major Head Loss - Friction Loss". Nuclear Power. Retrieved 17 September 2019.
20. "Laminar, transitional and turbulent flow". rheologic.net. Retrieved 17 September 2019.
21. Manneville & Pomeau 2009, p. 2072.
22. Patel, Rodi & Scheuerer 1985, pp. 1308-1319.
23. Dusenbery 2009, p. 136.
24. Los 2006, p. 369.
25. Kamenshchikov 2013, pp. 63-71.
26. Benson, Tom (12 June 2014). "Reynolds Number". Glenn Research Center. Retrieved 17 September 2019.

Sources

• Batchelor, G. K. (1967). An Introduction to Fluid Dynamics. Cambridge University Press. pp. 211–215.
• Fouz, Infaz "Fluid Mechanics," Mechanical Engineering Dept., University of Oxford, 2001, p. 96
• Hughes, Roger "Civil Engineering Hydraulics," Civil and Environmental Dept., University of Melbourne 1997, pp. 107–152
• Jermy M., "Fluid Mechanics A Course Reader," Mechanical Engineering Dept., University of Canterbury, 2005, pp. d5.10.
• Purcell, E. M. "Life at Low Reynolds Number", American Journal of Physics vol 45, pp. 3–11 (1977)
• Truskey, G. A., Yuan, F, Katz, D. F. (2004). Transport Phenomena in Biological Systems Prentice Hall, pp. 7. ISBN 0-13-042204-5. ISBN 978-0-13-042204-0.
• Zagarola, M. V. and Smits, A. J., "Experiments in High Reynolds Number Turbulent Pipe Flow." AIAA paper #96-0654, 34th AIAA Aerospace Sciences Meeting, Reno, Nevada, January 15–18, 1996.