Monin–Obukhov similarity theory

The Obukhov length ${\displaystyle L}$ 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,

${\displaystyle L=-{\dfrac {u_{*}^{3}}{\kappa {\dfrac {g}{T}}{\dfrac {Q}{\rho c_{p}}}}}}$

where ${\displaystyle \kappa \approx 0.40}$ is the von Kármán constant, ${\displaystyle u_{*}}$ friction velocity, ${\displaystyle Q}$ turbulent heat flux, and ${\displaystyle c_{p}}$ heat capacity.[4] Virtual potential temperature ${\displaystyle \theta _{v}}$ is often used instead of temperature ${\displaystyle T}$ to correct for the effects of pressure and water vapor. ${\displaystyle Q}$ can be written as vertical eddy flux,

${\displaystyle Q=\rho c_{p}{\overline {w'\theta _{v}'}}}$

with ${\displaystyle w'}$ and ${\displaystyle \theta _{v}'}$ perturbations of vertical velocity and virtual potential temperature, respectively. Therefore, the Obukhov length can also be defined as,[6]

${\displaystyle L=-{\dfrac {u_{*}^{3}}{\kappa {\dfrac {g}{\overline {\theta _{v}}}}{\overline {w'\theta _{v}'}}}}}$

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

M–O similarity theory parameterizes fluxes in the surface layer as a function of the dimensionless length parameter ${\displaystyle \zeta =z/L}$. From the Buckingham Pi theorem of dimensional analysis, two dimensionless group can be formed from the basic parameter set ${\displaystyle \{u_{*},g/{\overline {\theta _{v}}},\partial {\overline {u}}/\partial z,z,{\overline {w'\theta _{v}'}}\}}$,

${\displaystyle {\dfrac {\kappa z}{u_{*}}}{\dfrac {\partial {\overline {u}}}{\partial z}}}$, and

${\displaystyle \zeta ={\dfrac {z}{L}}}$

From there, a function ${\displaystyle \varphi _{M}(\zeta )}$ can be determined to empirically describe the relationship between the two dimensionless quantities, called a universal function. Similarly, ${\displaystyle \varphi _{H}(\zeta )}$ can be defined for the dimensionless group of mean temperature profile. Mean wind and temperature profiles therefore satisfy the following relations,[1][5]

${\displaystyle {\dfrac {\partial {\overline {u}}}{\partial z}}={\dfrac {u_{*}}{\kappa z}}\varphi _{M}(\zeta )}$

${\displaystyle {\dfrac {\partial {\overline {\theta _{v}}}}{\partial z}}={\dfrac {\theta _{*}}{\kappa z}}\varphi _{H}(\zeta )}$

where ${\displaystyle \theta _{*}=-{\dfrac {\overline {w'\theta _{v}'}}{u_{*}}}}$ is the characteristic dynamical temperature, ${\displaystyle \varphi _{M}}$ and ${\displaystyle \varphi _{H}}$ are the universal functions of momentum and heat. The eddy diffusivity coefficients for momentum and heat fluxes are defined as follows,

${\displaystyle K_{M}=\kappa z{\dfrac {u_{*}}{\varphi _{M}(\zeta )}},\ K_{H}=\kappa z{\dfrac {u_{*}}{\varphi _{H}(\zeta )}}}$

${\displaystyle K_{M}}$ and ${\displaystyle K_{H}}$ can be related with the turbulent Prandtl number ${\displaystyle Pr_{t}}$,

${\displaystyle {\dfrac {K_{H}}{K_{M}}}={\dfrac {1}{Pr_{t}}}>1}$

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 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 ${\displaystyle L{\Big (}{\dfrac {\partial Ri}{\partial z}}{\Big )}_{z=0}=1}$, where ${\displaystyle Ri}$ is the Richardson number, the following condition must be satisfied by the universal function chosen,[1]

${\displaystyle \varphi (0)=1}$

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

${\displaystyle \varphi _{M}(\zeta )=1+\beta \zeta }$

where ${\displaystyle \beta \approx 6}$.[5] However this is only applicable when ${\displaystyle 0<\zeta <1}$. For conditions where ${\displaystyle \zeta <0}$, the relation is,

${\displaystyle \varphi _{M}^{4}-\gamma \zeta \varphi _{M}^{3}=1}$

where ${\displaystyle \gamma }$ is a coefficient to be determined from experimental data. This equation can be further approximated by ${\displaystyle \varphi _{M}=(1+\gamma \zeta )^{-1/4}}$ when ${\displaystyle -2<\zeta <0}$.

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]

${\displaystyle \varphi _{M}(\zeta )=(1-15\zeta )^{-1/4}\quad -2<\zeta <0}$

${\displaystyle \varphi _{M}(\zeta )=1+4.7\zeta \quad 0<\zeta <1}$

${\displaystyle \varphi _{H}(\zeta )=0.74(1-9\zeta )^{-1/2}\quad -2<\zeta <0}$

${\displaystyle \varphi _{H}(\zeta )=0.74+4.7\zeta \quad 0<\zeta <1}$

Other methods which determine the universal functions using the relation between ${\displaystyle \zeta }$ and ${\displaystyle Ri}$ 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]

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

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]

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 ${\displaystyle {\dfrac {gQ}{T}}}$. 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 ${\displaystyle \epsilon _{0}}$.[12] This scheme is able to explain anomalies of M–O similarity theory, but involves non-locality to modeling and experiments.

• Obukhov length
• Surface layer
• Logarithmic wind profile
• Mixing length model