Rayleigh number

In fluid mechanics, the Rayleigh number (Ra, after Lord Rayleigh ^[1] ) for a fluid is a dimensionless number associated with buoyancy-driven flow, also known as free (or natural) convection. ^{[2] [3] [4]} It characterises the fluid's flow regime: ^[5] a value in a certain lower range denotes laminar flow; a value in a higher range, turbulent flow. Below a certain critical value, there is no fluid motion and heat transfer is by conduction rather than convection. For most engineering purposes, the Rayleigh number is large, somewhere around 10⁶ to 10⁸.

The Rayleigh number is defined as the product of the Grashof number (Gr), which describes the relationship between buoyancy and viscosity within a fluid, and the Prandtl number (Pr), which describes the relationship between momentum diffusivity and thermal diffusivity: Ra = Gr × Pr. ^{[4] [3]} Hence it may also be viewed as the ratio of buoyancy and viscosity forces multiplied by the ratio of momentum and thermal diffusivities: Ra = B/μ × ν/α. It is closely related to the Nusselt number (Nu). ^[5]

Derivation

The Rayleigh number describes the behaviour of fluids (such as water or air) when the mass density of the fluid is non-uniform. The mass density differences are usually caused by temperature differences. Typically a fluid expands and becomes less dense as it is heated. Gravity causes denser parts of the fluid to sink, which is called convection. Lord Rayleigh studied ^[2] the case of Rayleigh-Bénard convection. ^[6] When the Rayleigh number, Ra, is below a critical value for a fluid, there is no flow and heat transfer is purely by conduction; when it exceeds that value, heat is transferred by natural convection. ^[3]

When the mass density difference is caused by temperature difference, Ra is, by definition, the ratio of the time scale for diffusive thermal transport to the time scale for convective thermal transport at speed ${\displaystyle u}$: ^[4]

$${\displaystyle \mathrm {Ra} ={\frac {\text{time scale for thermal transport via diffusion}}{{\text{time scale for thermal transport via convection at speed}}~u}}.}$$

This means the Rayleigh number is a type ^[4] of Péclet number. For a volume of fluid of size ${\displaystyle l}$ in all three dimensions^{[ clarification needed ]} and mass density difference ${\displaystyle \Delta \rho }$, the force due to gravity is of the order ${\displaystyle \Delta \rho l^{3}g}$, where ${\displaystyle g}$ is acceleration due to gravity. From the Stokes equation, when the volume of fluid is sinking, viscous drag is of the order ${\displaystyle \eta lu}$, where ${\displaystyle \eta }$ is the dynamic viscosity of the fluid. When these two forces are equated, the speed ${\displaystyle u\sim \Delta \rho l^{2}g/\eta }$. Thus the time scale for transport via flow is ${\displaystyle l/u\sim \eta /\Delta \rho lg}$. The time scale for thermal diffusion across a distance ${\displaystyle l}$ is ${\displaystyle l^{2}/\alpha }$, where ${\displaystyle \alpha }$ is the thermal diffusivity. Thus the Rayleigh number Ra is

$${\displaystyle \mathrm {Ra} ={\frac {l^{2}/\alpha }{\eta /\Delta \rho lg}}={\frac {\Delta \rho l^{3}g}{\eta \alpha }}={\frac {\rho \beta \Delta Tl^{3}g}{\eta \alpha }}}$$

where we approximated the density difference ${\displaystyle \Delta \rho =\rho \beta \Delta T}$ for a fluid of average mass density ${\displaystyle \rho }$, thermal expansion coefficient ${\displaystyle \beta }$ and a temperature difference ${\displaystyle \Delta T}$ across distance ${\displaystyle l}$.

The Rayleigh number can be written as the product of the Grashof number and the Prandtl number: ^{[4] [3]} $${\displaystyle \mathrm {Ra} =\mathrm {Gr} \mathrm {Pr} .}$$

Classical definition

For free convection near a vertical wall, the Rayleigh number is defined as:

$${\displaystyle \mathrm {Ra} _{x}={\frac {g\beta }{\nu \alpha }}(T_{s}-T_{\infty })x^{3}=\mathrm {Gr} _{x}\mathrm {Pr} }$$

where:

• x is the characteristic length
• Ra_x is the Rayleigh number for characteristic length x
• g is acceleration due to gravity
• β is the thermal expansion coefficient (equals to 1/T, for ideal gases, where T is absolute temperature).
• ${\displaystyle \nu }$ is the kinematic viscosity
• α is the thermal diffusivity
• T_s is the surface temperature
• T_∞ is the quiescent temperature (fluid temperature far from the surface of the object)
• Gr_x is the Grashof number for characteristic length x
• Pr is the Prandtl number

In the above, the fluid properties Pr, ν, α and β are evaluated at the film temperature, which is defined as:

$${\displaystyle T_{f}={\frac {T_{s}+T_{\infty }}{2}}.}$$

For a uniform wall heating flux, the modified Rayleigh number is defined as:

$${\displaystyle \mathrm {Ra} _{x}^{*}={\frac {g\beta q''_{o}}{\nu \alpha k}}x^{4}}$$

where:

• q″_o is the uniform surface heat flux
• k is the thermal conductivity. ^[7]

Other applications

Solidifying alloys

The Rayleigh number can also be used as a criterion to predict convectional instabilities, such as A-segregates, in the mushy zone of a solidifying alloy. The mushy zone Rayleigh number is defined as:

$${\displaystyle \mathrm {Ra} ={\frac {{\frac {\Delta \rho }{\rho _{0}}}g{\bar {K}}L}{\alpha \nu }}={\frac {{\frac {\Delta \rho }{\rho _{0}}}g{\bar {K}}}{R\nu }}}$$

where:

• K is the mean permeability (of the initial portion of the mush)
• L is the characteristic length scale
• α is the thermal diffusivity
• ν is the kinematic viscosity
• R is the solidification or isotherm speed. ^[8]

A-segregates are predicted to form when the Rayleigh number exceeds a certain critical value. This critical value is independent of the composition of the alloy, and this is the main advantage of the Rayleigh number criterion over other criteria for prediction of convectional instabilities, such as Suzuki criterion.

Torabi Rad et al. showed that for steel alloys the critical Rayleigh number is 17. ^[8] Pickering et al. explored Torabi Rad's criterion, and further verified its effectiveness. Critical Rayleigh numbers for lead–tin and nickel-based super-alloys were also developed. ^[9]

Porous media

The Rayleigh number above is for convection in a bulk fluid such as air or water, but convection can also occur when the fluid is inside and fills a porous medium, such as porous rock saturated with water. ^[10] Then the Rayleigh number, sometimes called the Rayleigh-Darcy number, is different. In a bulk fluid, i.e., not in a porous medium, from the Stokes equation, the falling speed of a domain of size ${\displaystyle l}$ of liquid ${\displaystyle u\sim \Delta \rho l^{2}g/\eta }$. In porous medium, this expression is replaced by that from Darcy's law ${\displaystyle u\sim \Delta \rho kg/\eta }$, with ${\displaystyle k}$ the permeability of the porous medium. The Rayleigh or Rayleigh-Darcy number is then

$${\displaystyle \mathrm {Ra} ={\frac {\rho \beta \Delta Tklg}{\eta \alpha }}}$$

This also applies to A-segregates, in the mushy zone of a solidifying alloy. ^[8]

Geophysical applications

In geophysics, the Rayleigh number is of fundamental importance: it indicates the presence and strength of convection within a fluid body such as the Earth's mantle. The mantle is a solid that behaves as a fluid over geological time scales. The Rayleigh number for the Earth's mantle due to internal heating alone, Ra_H, is given by:

$${\displaystyle \mathrm {Ra} _{H}={\frac {g\rho _{0}^{2}\beta HD^{5}}{\eta \alpha k}}}$$

where:

• H is the rate of radiogenic heat production per unit mass
• η is the dynamic viscosity
• k is the thermal conductivity
• D is the depth of the mantle. ^[11]

A Rayleigh number for bottom heating of the mantle from the core, Ra_T, can also be defined as:

$${\displaystyle \mathrm {Ra} _{T}={\frac {\rho _{0}^{2}g\beta \Delta T_{\text{sa}}D^{3}C_{P}}{\eta k}}}$$

where:

• ΔT_sa is the superadiabatic temperature difference (the superadiabatic temperature difference is the actual temperature difference minus the temperature difference in a fluid whose entropy gradient is zero, but has the same profile of the other variables appearing in the equation of state) between the reference mantle temperature and the core–mantle boundary
• C_P is the specific heat capacity at constant pressure. ^[11]

High values for the Earth's mantle indicates that convection within the Earth is vigorous and time-varying, and that convection is responsible for almost all the heat transported from the deep interior to the surface.

See also

• Grashof number
• Prandtl number
• Reynolds number
• Péclet number
• Nusselt number
• Rayleigh–Bénard convection

Notes

1. Chandrasekhar, S. (1961). Hydrodynamic and Hydromagnetic Stability . London: Oxford University Press. p.  10. ISBN   978-0-19-851237-0.
2. Baron Rayleigh (1916). "On convection currents in a horizontal layer of fluid, when the higher temperature is on the under side". London Edinburgh Dublin Phil. Mag. J. Sci. 32 (192): 529–546. doi:10.1080/14786441608635602.
3. Çengel, Yunus; Turner, Robert; Cimbala, John (2017). Fundamentals of thermal-fluid sciences (Fifth ed.). New York, NY. ISBN   9780078027680. OCLC   929985323.
4. Squires, Todd M.; Quake, Stephen R. (2005-10-06). "Microfluidics: Fluid physics at the nanoliter scale" (PDF). Reviews of Modern Physics. 77 (3): 977–1026. Bibcode:2005RvMP...77..977S. doi:10.1103/RevModPhys.77.977.
5. Çengel, Yunus A. (2002). Heat and Mass Transfer (Second ed.). McGraw-Hill. p. 466.
6. Ahlers, Guenter; Grossmann, Siegfried; Lohse, Detlef (2009-04-22). "Heat transfer and large scale dynamics in turbulent Rayleigh-Bénard convection". Reviews of Modern Physics. 81 (2): 503–537. arXiv: 0811.0471 . Bibcode:2009RvMP...81..503A. doi:10.1103/RevModPhys.81.503. S2CID   7566961.
7. M. Favre-Marinet and S. Tardu, Convective Heat Transfer, ISTE, Ltd, London, 2009
8. Torabi Rad, M.; Kotas, P.; Beckermann, C. (2013). "Rayleigh number criterion for formation of A-Segregates in steel castings and ingots". Metall. Mater. Trans. A. 44A (9): 4266–4281. Bibcode:2013MMTA...44.4266R. doi:10.1007/s11661-013-1761-4. S2CID   137652216.
9. Pickering, E.J.; Al-Bermani, S.; Talamantes-Silva, J. (2014). "Application of criterion for A-segregation in steel ingots". Materials Science and Technology. 31 (11): 1313. Bibcode:2015MatST..31.1313P. doi:10.1179/1743284714Y.0000000692. S2CID   137549220.
10. Lister, John R.; Neufeld, Jerome A.; Hewitt, Duncan R. (2014). "High Rayleigh number convection in a three-dimensional porous medium". Journal of Fluid Mechanics. 748: 879–895. arXiv: 0811.0471 . Bibcode:2014JFM...748..879H. doi:10.1017/jfm.2014.216. ISSN   1469-7645. S2CID   43758157.
11. Bunge, Hans-Peter; Richards, Mark A.; Baumgardner, John R. (1997). "A sensitivity study of three-dimensional spherical mantle convection at 10⁸ Rayleigh number: Effects of depth-dependent viscosity, heating mode, and endothermic phase change". Journal of Geophysical Research . 102 (B6): 11991–12007. Bibcode:1997JGR...10211991B. doi: 10.1029/96JB03806 .

References

• Turcotte, D.; Schubert, G. (2002). Geodynamics (2nd ed.). New York: Cambridge University Press. ISBN   978-0-521-66186-7.

External links

• Rayleigh number calculator