Monin–Obukhov similarity theory

Monin–Obukhov (M–O) similarity theory describes the non-dimensionalized mean flow and mean temperature in the surface layer under non-neutral conditions as a function of the dimensionless height parameter,[1] named after Russian scientists A. S. Monin and A. M. Obukhov. Similarity theory is an empirical method that describes universal relationships between non-dimensionalized variables of fluids based on the Buckingham π theorem. Similarity theory is extensively used in boundary layer meteorology since relations in turbulent processes are not always resolvable from first principles.[2]

An idealized vertical profile of the mean flow for a neutral boundary layer is the logarithmic wind profile derived from Prandtl's mixing length theory,[3] which states that the horizontal component of mean flow is proportional to the logarithm of height. M–O similarity theory further generalizes the mixing length theory in non-neutral conditions by using so-called "universal functions" of dimensionless height to characterize vertical distributions of mean flow and temperature. The Obukhov length (), a characteristic length scale of surface layer turbulence derived by Obukhov in 1946,[4] is used for non-dimensional scaling of the actual height. M–O similarity theory marked a significant landmark of modern micrometeorology, providing a theoretical basis for micrometeorological experiments and measurement techniques.[5]

The Obukhov length

edit

The Obukhov length   is a length parameter for the surface layer in the boundary layer, which characterizes the relative contributions to turbulent kinetic energy from buoyant production and shear production. The Obukhov length was formulated using Richardson's criterion for dynamic stability.[4] It was derived as,

 

where   is the von Kármán constant,   friction velocity,   turbulent heat flux, and   heat capacity.[4] Virtual potential temperature   is often used instead of temperature   to correct for the effects of pressure and water vapor.   can be written as vertical eddy flux,

 

with   and   perturbations of vertical velocity and virtual potential temperature, respectively. Therefore, the Obukhov length can also be defined as,[6]

 

The Obukhov length also acts as a criterion for the static stability of surface layer. When  , the surface layer is statically unstable, and when   the surface layer is statically stable. The absolute magnitude of   indicates the deviation from statically neutral state, with smaller   values corresponding to larger deviations from neutral conditions. When   is small and  , buoyant processes dominate the production of turbulent kinetic energy compared with shear production. By definition, under neutral conditions  . The Obukhov length   is used for non-dimensionalization of height   in similarity theory.

Governing formulae for similarity relations

edit

M–O similarity theory parameterizes fluxes in the surface layer as a function of the dimensionless length parameter  . From the Buckingham Pi theorem of dimensional analysis, two dimensionless group can be formed from the basic parameter set  ,

 , and

 

From there, a function   can be determined to empirically describe the relationship between the two dimensionless quantities, called a universal function. Similarly,   can be defined for the dimensionless group of mean temperature profile. Mean wind and temperature profiles therefore satisfy the following relations,[1][5]

 
 

where   is the characteristic dynamical temperature,   and   are the universal functions of momentum and heat. The eddy diffusivity coefficients for momentum and heat fluxes are defined as follows,

 

  and   can be related with the turbulent Prandtl number  ,

 

In reality, the universal functions need to be determined using experimental data when applying M–O similarity theory. Although the choice of universal functions is not unique, certain functional forms have been proposed and are widely accepted for fitting experimental data.

Universal functions of the Monin–Obukhov similarity theory

edit
 
Universal functions for Monin–Obukhov similarity theory

Several functional forms have been proposed to represent the universal functions of similarity theory. Because the Obukhov length is determined when  , where   is the Richardson number, the following condition must be satisfied by the universal function chosen,[1]

 

A first order approximation of the universal function for momentum flux is,

 

where  .[5] However this is only applicable when  . For conditions where  , the relation is,

 

where   is a coefficient to be determined from experimental data. This equation can be further approximated by   when  .

Based on the results of the 1968 Kansas experiment, the following universal functions are determined for horizontal mean flow and mean virtual potential temperature,[7]

 
 
 
 

Other methods which determine the universal functions using the relation between   and   are also used.[8][9]

For sublayers with significant roughness, e.g. vegetated surfaces or urban areas, the universal functions must be modified to include the effects of surface roughness.[6]

Validations

edit

A myriad of experimental efforts was devoted to the validation of the M–O similarity theory. Field observations and computer simulations have generally demonstrated that the M–O similarity theory is well satisfied.

In field measurements

edit
 
A Kansas wheat field, the flat terrain is needed for the experiment

The 1968 Kansas experiment found great consistency between measurements and predictions from similarity relations for the entire range of stability values.[7] A flat wheat field in Kansas served as the experiment site, with winds measured by anemometers mounted at different heights on a 32 m tower. The temperature profile was also measured in a similar manner. Results from the Kansas field study indicated that the ratio of eddy diffusivities of heat and momentum was approximately 1.35 under neutral conditions. A similar experiment was conducted in a flat field in northwestern Minnesota in 1973. This experiment used both ground and balloon-based observations of the surface layer and further validated the theoretical predictions from similarity.[10]

In large eddy simulations

edit

In addition to field experiments, analysis of M–O similarity theory can be conducted using high-resolution large eddy simulations. The simulation indicates that the temperature field agrees well with M–O similarity. However, the velocity field shows significant anomalies from M–O similarity.[11]

Limitations

edit

M–O similarity theory, albeit successful for surface layers from experimental validations, is essentially a diagnostic empirical theory based upon local first order turbulence closure. Typically, 10%~20% of errors are associated with universal functions. When applied to vegetated areas or complex terrains, it can result in large discrepancies. Because universal functions are often determined under dry conditions, the applicability of M–O similarity theory under moist conditions was not well studied.

The basic parameter set of the M–O similarity theory includes buoyancy production  . It is argued that with such a parameter set, the scaling is applied to the integral features of the flow, whereas an eddy specific similarity relationship prefers the usage of energy dissipation rate  .[12] This scheme is able to explain anomalies of M–O similarity theory, but involves non-locality to modeling and experiments.

See also

edit

References

edit
  1. ^ a b c Monin, A. S.; Obukhov, A. M. (1954). "Basic laws of turbulent mixing in the surface layer of the atmosphere". Tr. Akad. Nauk. SSSR Geophiz. Inst. 24 (151): 163–187.
  2. ^ Stull, Roland (1988). An introduction to boundary layer meteorology. Netherlands: Springer. ISBN 978-94-009-3027-8.
  3. ^ Prandtl, Ludwig (1925). "Bericht über Untersuchungen zur ausgebildeten Turbulenz". Zeitschrift für Angewandte Mathematik und Mechanik. 5 (2): 136–139. Bibcode:1925ZaMM....5..136P. doi:10.1002/zamm.19250050212.
  4. ^ a b c Obukhov, A. M. (1971). "Turbulence in an atmosphere with a non-uniform temperature". Boundary-Layer Meteorology. 2 (1): 7–29. Bibcode:1971BoLMe...2....7O. doi:10.1007/BF00718085. S2CID 121123105.
  5. ^ a b c Foken, T. (2006). "50 Years of the Monin-Obukhov similarity theory". Boundary-Layer Meteorology. 2 (3): 7–29. Bibcode:2006BoLMe.119..431F. doi:10.1007/s10546-006-9048-6. S2CID 122060208.
  6. ^ a b Foken, Thomas (2008). Micrometeorology. Springer-Verlag. pp. 42–49. ISBN 978-3-540-74665-2.
  7. ^ a b Businger, J. A.; J. C. Wyngaard; Y. Izumi; E. F. Bradley (1971). "Flux-profile relationships in the atmospheric surface layer". Journal of the Atmospheric Sciences. 28 (2): 181–189. Bibcode:1971JAtS...28..181B. doi:10.1175/1520-0469(1971)028<0181:FPRITA>2.0.CO;2.
  8. ^ Arya, S. P. (2001). Introduction to Micrometeorology. San Diego: Academic Press.
  9. ^ Högström, U. (1988). "Non-dimensional wind and temperature profiles in the atmospheric surface layer: A re-evaluation". Boundary-Layer Meteorology. 42 (1–2): 55–78. Bibcode:1988BoLMe..42...55H. doi:10.1007/BF00119875. S2CID 117742460.
  10. ^ Kaimal, J. C.; J. C. Wyngaard; D. A. Haugen; O. R. Coté; Y. Izumi; S. J. Caughey; C. J. Readings (1976). "Turbulence Structure in the Convective Boundary Layer". Journal of the Atmospheric Sciences. 33 (11): 2152–2169. Bibcode:1976JAtS...33.2152K. doi:10.1175/1520-0469(1976)033<2152:TSITCB>2.0.CO;2.
  11. ^ Khanna, Samir; Brasseur, James G. (1997). "Analysis of Monin–Obukhov similarity from large-eddy simulation". J. Fluid Mech. 345 (1): 251–286. Bibcode:1997JFM...345..251K. doi:10.1017/S0022112097006277. S2CID 122546650.
  12. ^ McNaughton, Keith (2009). "The rise and fall of Monin-Obukhov theory" (PDF). AsiaFlux Newsletter (30): 1–4.