Richardson number

From HandWiki

The Richardson number (Ri) is named after Lewis Fry Richardson (1881–1953).[1] It is the dimensionless number that expresses the ratio of the buoyancy term to the flow shear term:[2]

[math]\displaystyle{ \mathrm{Ri} = \frac{\text{buoyancy term}}{\text{flow shear term}} = \frac{g}{\rho} \frac{\partial \rho/\partial z}{(\partial u / \partial z)^2} }[/math]

where [math]\displaystyle{ g }[/math] is gravity, [math]\displaystyle{ \rho }[/math] is density, [math]\displaystyle{ u }[/math] is a representative flow speed, and [math]\displaystyle{ z }[/math] is depth.

The Richardson number, or one of several variants, is of practical importance in weather forecasting and in investigating density and turbidity currents in oceans, lakes, and reservoirs.

When considering flows in which density differences are small (the Boussinesq approximation), it is common to use the reduced gravity g' and the relevant parameter is the densimetric Richardson number[further explanation needed]

[math]\displaystyle{ \mathrm{Ri} = \frac{-\partial g'/\partial z}{(\partial u / \partial z)^2} }[/math]

which is used frequently when considering atmospheric or oceanic flows[citation needed].

If the Richardson number is much less than unity, buoyancy is unimportant in the flow. If it is much greater than unity, buoyancy is dominant (in the sense that there is insufficient kinetic energy to homogenize the fluids).

If the Richardson number is of order unity, then the flow is likely to be buoyancy-driven: the energy of the flow derives from the potential energy in the system originally.

Aviation

In aviation, the Richardson number is used as a rough measure of expected air turbulence. A lower value indicates a higher degree of turbulence. Values in the range 10 to 0.1 are typical[citation needed], with values below unity indicating significant turbulence.

Thermal convection

In thermal convection problems, Richardson number represents the importance of natural convection relative to the forced convection. The Richardson number in this context is defined as

[math]\displaystyle{ \mathrm{Ri} = \frac{g \beta (T_\text{hot} - T_\text{ref})L}{V^2} }[/math]

where g is the gravitational acceleration, [math]\displaystyle{ \beta }[/math] is the thermal expansion coefficient, Thot is the hot wall temperature, Tref is the reference temperature, L is the characteristic length, and V is the characteristic velocity.

The Richardson number can also be expressed by using a combination of the Grashof number and Reynolds number,

[math]\displaystyle{ \mathrm{Ri} = \frac{\mathrm{Gr}}{\mathrm{Re}^2}. }[/math]

Typically, the natural convection is negligible when Ri < 0.1, forced convection is negligible when Ri > 10, and neither is negligible when 0.1 < Ri < 10. It may be noted that usually the forced convection is large relative to natural convection except in the case of extremely low forced flow velocities. However, buoyancy often plays a significant role in defining the laminar–turbulent transition of a mixed convection flow.[3] In the design of water filled thermal energy storage tanks, the Richardson number can be useful.[4]

Meteorology

In atmospheric science, several different expressions for the Richardson number are commonly used: the flux Richardson number (which is fundamental), the gradient Richardson number, and the bulk Richardson number.

  • The flux Richardson number [math]\displaystyle{ Ri_{f} }[/math] is the ratio of buoyant production of turbulence kinetic energy to the suppression of turbulence by shear.[5] Mathematically, this is:
[math]\displaystyle{ Ri_{f} = \frac{(g/T_{v}) \overline{w' \theta'}}{\overline{u' w'} \frac{\partial u_{z}}{\partial z} + \overline{v' w'} \frac{\partial v_{z}}{\partial z}} }[/math],

where [math]\displaystyle{ T_{v} }[/math] is the virtual temperature, [math]\displaystyle{ \theta_{v} }[/math] is the virtual potential temperature, [math]\displaystyle{ z }[/math] is the altitude, [math]\displaystyle{ u }[/math] is the [math]\displaystyle{ x }[/math] component of the wind, [math]\displaystyle{ v }[/math] is the [math]\displaystyle{ y }[/math] component of the wind, and [math]\displaystyle{ w }[/math] is the [math]\displaystyle{ z }[/math] (vertical) component of the wind. A prime (e.g. [math]\displaystyle{ w' }[/math]) denotes a deviation of the respective field from its Reynolds average.


  • The gradient Richardson number [math]\displaystyle{ Ri_{g} }[/math] is arrived at by approximating the flux Richardson number using "K-theory". This results in:[6]
[math]\displaystyle{ Ri_{g} = \frac{(g/T_{v}) \frac{\partial \theta_{v}}{\partial z}}{(\frac{\partial u}{\partial z})^2 + (\frac{\partial v}{\partial z})^{2}} }[/math].
  • The bulk Richardson number [math]\displaystyle{ Ri_{b} }[/math] results from making a finite difference approximation to the derivatives in the expression for the gradient Richardson number, giving:[7]
[math]\displaystyle{ Ri_{b} = \frac{(g/T_{v0})\Delta \theta_{v} \Delta z}{(\Delta u)^2 + (\Delta v)^2} }[/math].

Here, for any variable [math]\displaystyle{ f }[/math], [math]\displaystyle{ \Delta f := f_{z1} - f_{z0} }[/math], i.e. the difference between [math]\displaystyle{ f }[/math] at altitude [math]\displaystyle{ z1 }[/math] and altitude [math]\displaystyle{ z0 }[/math]. If the lower reference level is taken to be [math]\displaystyle{ z0=0 }[/math], then [math]\displaystyle{ u_{z0}=v_{z0}=0 }[/math] (due to the no-slip boundary condition), so the expression simplifies to:

[math]\displaystyle{ Ri_{b} = \frac{(g/\theta_{v0})(\theta_{vz1} - \theta_{v0})z}{(u_{z1})^2 + (v_{z1})^{2}} }[/math].

Oceanography

In oceanography, the Richardson number has a more general form[citation needed] which takes stratification into account. It is a measure of relative importance of mechanical and density effects in the water column, as described by the Taylor–Goldstein equation, used to model Kelvin–Helmholtz instability which is driven by sheared flows.

[math]\displaystyle{ \mathrm{Ri} = \frac{N^2}{(\mathrm{d}u/\mathrm{d}z)^2} }[/math]

where N is the Brunt–Väisälä frequency and u the wind speed.

The Richardson number defined above is always considered positive. A negative value of (i.e. complex N) indicates unstable density gradients with active convective overturning. Under such circumstances the magnitude of negative Ri is not generally of interest. It can be shown that Ri < 1/4 is a necessary condition for velocity shear to overcome the tendency of a stratified fluid to remain stratified, and some mixing (turbulence) will generally occur. When Ri is large, turbulent mixing across the stratification is generally suppressed.[8]

References

  1. Hunt, J.C.R. (1998). "Lewis Fry Richardson and His Contributions to Mathematics, Meteorology, and Models of Conflict" (in en). Annual Review of Fluid Mechanics 30 (1): xiii–xxxvi. doi:10.1146/annurev.fluid.30.1.0. ISSN 0066-4189. Bibcode1998AnRFM..30D..13H. https://www.annualreviews.org/doi/full/10.1146/annurev.fluid.30.1.0. 
  2. Encyclopædia Britannica: Richardson number
  3. Garbrecht, Oliver (August 23, 2017). "Large eddy simulation of three-dimensional mixed convection on a vertical plate". RWTH Aachen University. http://publications.rwth-aachen.de/record/718097/files/718097.pdf. 
  4. Robert Huhn Beitrag zur thermodynamischen Analyse und Bewertung von Wasserwärmespeichern in Energieumwandlungsketten, ISBN:978-3-940046-32-1, Andreas Oberhammer Europas größter Fernwärmespeicher in Kombination mit dem optimalen Ladebetrieb eines Gas- und Dampfturbinenkraftwerkes (Vortrag 2007)
  5. "Flux Richardson number". AMS Glossary. American Meteorological Society. https://glossary.ametsoc.org/wiki/Flux_richardson_number. 
  6. "Gradient richardson number". AMS Glossary. American Meteorological Society. https://glossary.ametsoc.org/wiki/Gradient_richardson_number. 
  7. "Bulk richardson number". AMS Glossary. American Meteorological Society. https://glossary.ametsoc.org/wiki/Bulk_richardson_number. 
  8. A good reference on this subject is Turner, J. S. (1973). Buoyancy Effects in Fluids. Cambridge University Press. ISBN 978-0-521-08623-3. https://archive.org/details/buoyancyeffectsi0000turn.