Reynolds number

In fluid dynamics, the Reynolds number (Re) is a dimensionless quantity that helps predict fluid flow patterns in different situations by measuring the ratio between inertial and viscous forces.^[2] At low Reynolds numbers, flows tend to be dominated by laminar (sheet-like) flow, while at high Reynolds numbers, flows tend to be turbulent. The 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 transition from laminar (left) to turbulent (right) flow of water from a tap occurs as the Reynolds number increases.

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 behavior on a larger scale, such as in local or global air or water movement, and thereby the associated meteorological and climatological effects.

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 who popularized its use in 1883.^[5][6] (cf. this list)

The Reynolds number is the ratio of inertial forces to viscous forces within a fluid that is subjected to relative internal movement due to different fluid velocities. A region where these forces change behavior is known as a boundary layer, such as the bounding surface in 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]

where:

• ρ is the density of the fluid (SI units: kg/m³)
• u is the flow speed (m/s)
• L is a characteristic length (m)
• μ is the dynamic viscosity of the fluid (Pa·s or N·s/m² or kg/(m·s))
• ν is the kinematic viscosity of the fluid (m²/s).

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.

If we know that the relevant physical quantities in a physical system are only ${\displaystyle \rho ,u,L,\mu }$, then the Reynolds number is essentially fixed by the Buckingham π theorem.

In detail, since there are 4 quantities ${\displaystyle \rho ,u,L,\mu }$, but they have only 3 dimensions (length, time, mass), we can consider ${\displaystyle \rho ^{x_{1}}u^{x_{2}}L^{x_{3}}\mu ^{x_{4}}}$, where ${\displaystyle x_{1},...,x_{4}}$ are real numbers. Setting the three dimensions of ${\displaystyle \rho ^{x_{1}}u^{x_{2}}L^{x_{3}}\mu ^{x_{4}}}$ to zero, we obtain 3 independent linear constraints, so the solution space has 1 dimension, and it is spanned by the vector ${\displaystyle (1,1,1,-1)}$.

Thus, any dimensionless quantity constructed out of ${\displaystyle \rho ,u,L,\mu }$ is a function of ${\displaystyle \rho uL\mu ^{-1}}$, the Reynolds number.

Alternatively, we can take the incompressible Navier–Stokes equations (convective form):$${\displaystyle {\frac {\partial \mathbf {u} }{\partial t}}+(\mathbf {u} \cdot \nabla )\mathbf {u} -\nu \,\nabla ^{2}\mathbf {u} =-{\frac {1}{\rho }}\nabla p+\mathbf {g} }$$Remove the gravity term ${\displaystyle {\mathbf {g}}}$, then the left side consists of inertial force ${\displaystyle {\frac {\partial \mathbf {u} }{\partial t}}+(\mathbf {u} \cdot \nabla )\mathbf {u} }$, and viscous force ${\displaystyle \nu \,\nabla ^{2}\mathbf {u} }$. Their ratio has the order of ${\displaystyle {\frac {(\mathbf {u} \cdot \nabla )\mathbf {u} }{\nu \,\nabla ^{2}\mathbf {u} }}\sim {\frac {u^{2}/L}{\nu u/L^{2}}}={\frac {uL}{\nu }}}$, the Reynolds number. This argument is written out in detail on the Scallop theorem page.

Alternative 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/m³).

If we now set

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

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.

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. 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 throughout 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.

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

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

where

• D_H is the hydraulic diameter of the pipe (the inside diameter if the pipe is circular) (m),
• Q is the volumetric flow rate (m³/s),
• A is the pipe's cross-sectional area (A = ⁠πD²/4⁠) (m²),
• u is the mean velocity of the fluid (m/s),
• μ (mu) is the dynamic viscosity of the fluid (Pa·s = N·s/m² = kg/(m·s)),
• ν (nu) is the kinematic viscosity (ν = ⁠μ/ρ⁠) (m²/s),
• ρ (rho) is the density of the fluid (kg/m³),
• W is the mass flowrate of the fluid (kg/s).

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, D_H, 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

• D_o is the inside diameter of the outer pipe,
• D_i 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]

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 Re_x ≈ 5×10⁵,^[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 Re_D < 2300 and turbulent flow occurs when Re_D > 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 Re_D > 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]

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 the distance between the plates.^[16] This is consistent with the annular duct and rectangular duct cases above, taken to a limiting aspect ratio.

For calculating the 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 a quarter of 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,^[17] while others use the hydraulic radius as the characteristic length-scale with consequently different values of Re for transition and turbulent flow.

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).^[18] Fluid dynamicists define the chord Reynolds number 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⁻⁵ m²/s for the atmosphere at sea level.^[19] 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 span-wise stations on a wing, where chord is still used.^[20]

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

In viscous fluids

This section does not cite any sources. (February 2014)

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.^[21]

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.^[22] 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 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.^[23] 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.^[24]

At higher Reynolds numbers the drag on a sphere depends on surface roughness. Thus, for example, adding dimples on the surface of a golf ball causes the boundary layer on the upstream side of the ball to transition from laminar to turbulent. The turbulent boundary layer is able to remain attached to the surface of the ball much longer than a laminar boundary and so creates a narrower low-pressure wake and hence less pressure drag. The reduction in pressure drag causes the ball to travel farther.^[25]

Pressure drops^[28] 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.

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 the 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.^[29][30] 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^[31][32]

• Dictyostelium amoebae: ~ 1 × 10^−6[33]
• Bacterium ~ 1 × 10⁻⁴
• Ciliate ~ 1 × 10⁻¹
• Smallest fish ~ 1
• Blood flow in brain ~ 1 × 10²
• Blood flow in aorta ~ 1 × 10³
• Onset of turbulent flow ~ 2.3 × 10³ to 5.0 × 10⁴ for pipe flow to 10⁶ for boundary layers
• Typical pitch in Major League Baseball ~ 2 × 10⁵
• Person swimming ~ 4 × 10⁶
• Fastest fish ~ 1 × 10⁸
• Blue whale ~ 4 × 10⁸
• A large ship (Queen Elizabeth 2) ~ 5 × 10⁹
• Atmospheric tropical cyclone ~ 1 x 10¹²

This section does not cite any sources. (January 2019)

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 smokestack, 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.^[34]

This section does not cite any sources. (December 2014)

Poiseuille's law on blood circulation in the body is dependent on laminar flow.^[35] 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.^[36]

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

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.

• Reynolds transport theorem – 3D generalization of the Leibniz integral rule
• Drag coefficient – Dimensionless parameter to quantify fluid resistance
• Deposition (geology) – Geological process in which sediments, soil and rocks are added to a landform or landmass
• Kelvin–Helmholtz instability