the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# On dissipation timescales of the basic second-order moments: the effect on the energy and flux budget (EFB) turbulence closure for stably stratified turbulence

### Evgeny Kadantsev

### Evgeny Mortikov

### Andrey Glazunov

### Nathan Kleeorin

### Igor Rogachevskii

The dissipation rates of the basic second-order moments are the key parameters playing a vital role in turbulence modelling and controlling turbulence energetics and spectra and turbulent fluxes of momentum and heat. In this paper, we use the results of direct numerical simulations (DNSs) to evaluate dissipation rates of the basic second-order moments and revise the energy and flux budget (EFB) turbulence closure theory for stably stratified turbulence. We delve into the theoretical implications of this approach and substantiate our closure hypotheses through DNS data. We also show why the concept of down-gradient turbulent transport becomes incomplete when applied to the vertical turbulent flux of potential temperature under stable stratification. We reveal essential feedback between the turbulent kinetic energy (TKE), the vertical turbulent flux of buoyancy, and the turbulent potential energy (TPE), which is responsible for maintaining shear-produced stably stratified turbulence for any Richardson number.

- Article
(2014 KB) - Full-text XML
- BibTeX
- EndNote

Turbulence and associated turbulent transport have been studied theoretically, experimentally, observationally, and numerically throughout several decades (see books by Batchelor, 1953; Monin and Yaglom, 1971, 2013; Tennekes and Lumley, 1972; Frisch, 1995; Pope, 2000; Davidson, 2013; Rogachevskii, 2021, and references therein), but some important questions remain. This is particularly true in applications to atmospheric physics and geophysics, where Reynolds and Péclet numbers are extremely large, so the governing equations are strongly nonlinear. The classical Kolmogorov theory (Kolmogorov, 1941a, b, 1942, 1991) has been formulated for neutrally stratified homogeneous and isotropic turbulence.

In atmospheric boundary layers, temperature stratification causes turbulence to become anisotropic and inhomogeneous, making some assumptions underlying Kolmogorov's theory questionable. Numerous alternative turbulence closure theories (see reviews by Weng and Taylor, 2003; Umlauf and Burchard, 2005; Mahrt, 2014) have been formulated using the budget equations for not only turbulent kinetic energy (TKE), but also turbulent potential energy (TPE) (see, for example, Holloway, 1986; Ostrovsky and Troitskaya, 1987; Dalaudier and Sidi, 1987; Hunt et al., 1988; Canuto and Minotti, 1993; Schumann and Gerz, 1995; Hanazaki and Hunt, 1996; Keller and van Atta, 2000; Canuto et al., 2001; Stretch et al., 2001; Cheng et al., 2002; Hanazaki and Hunt, 2004; Rehmann and Hwang, 2005; Umlauf, 2005). The budget equations for all three energies, TKE, TPE, and total turbulent energy (TTE), were considered by Canuto and Minotti (1993), Elperin et al. (2002, 2006), Zilitinkevich et al. (2007), and Canuto et al. (2008).

The energy and flux budget (EFB) turbulence closure theory, which is based on the budget equations for the densities of TKE, TPE, and turbulent fluxes of momentum and heat, has been developed for the stably stratified atmospheric flows (Zilitinkevich et al., 2007, 2008, 2009, 2013, 2019; Kleeorin et al., 2019), surface layers in atmospheric convective turbulence (Rogachevskii et al., 2022), and core of the convective boundary layer (Rogachevskii and Kleeorin, 2024), as well as for passive scalar transport (Kleeorin et al., 2021). The EFB closure theory has shown that strong atmospheric stably stratified turbulence is maintained by large-scale shear (mean wind) for any stratification, and the “critical” Richardson number, considered for many years to be a threshold between the turbulent and laminar states of the flow, actually separates two turbulent regimes: the strong turbulence typical of atmospheric boundary layers and the weak three-dimensional turbulence typical of the free atmosphere and characterised by a strong decrease in the turbulent heat transfer in comparison to the momentum transfer.

Some other turbulent closure models (Mauritsen et al., 2007; Canuto et al., 2008; Sukoriansky and Galperin, 2008; Li et al., 2016) do not employ the concept of the critical Richardson number, so shear-generated turbulent mixing may persist for any stratification. In particular, Mauritsen et al. (2007) have developed a turbulent closure based on the budget equation for TTE (instead of TKE) and different observational findings to take into account the mean flow stability. They used this turbulent closure model to study the turbulent transfer of heat and momentum under very stable stratification. In their model, although the turbulent heat flux tends toward zero beyond a certain stability limit, the turbulent stress stays finite. However, the model by Mauritsen et al. (2007) does not use the budget equation for TPE and the vertical turbulent heat flux.

L'vov et al. (2008) performed detailed analyses of the budget equations for the Reynolds stresses in the turbulent boundary layer (relevant to the strong turbulence regime), explicitly taking into consideration the dissipative effect in the horizontal turbulent heat flux budget equation in contrast to the EFB “effective dissipation approximation”, adopted in the EFB turbulent closure model. However, the theory by L'vov et al. (2008) still contains the critical gradient Richardson number for the existence of the shear-produced turbulence.

Sukoriansky and Galperin (2008) apply a quasi-normal-scale elimination theory that is similar to the renormalization group analysis. Sukoriansky and Galperin (2008) do not use the budget equations for TKE, TPE, and TTE in their analysis. This theory correctly describes the dependence of the turbulent Prandtl number versus the gradient Richardson number and does not employ the concept of the critical gradient Richardson number for the existence of turbulence. However, this approach does not have detailed Richardson number dependences of the other non-dimensional parameters, such as the ratio between TPE and TTE, dimensionless turbulent flux of momentum, or dimensionless vertical turbulent flux of potential temperature. Their background non-stratified shear-produced turbulence is assumed to be isotropic and homogeneous. Canuto et al. (2008) have generalised their original model (see Cheng et al., 2002), introducing the new parameterisation for the buoyancy timescale to accommodate the existence of stably stratified shear-produced turbulence at arbitrary Richardson numbers.

Li et al. (2016) have developed the co-spectral budget (CSB) closure approach, which is formulated in the Fourier space and integrated across all turbulent scales to obtain turbulent characteristics in physical space. The CSB model allows turbulence to exist at any gradient Richardson number. However, the CSB model yields different predictions for the vertical anisotropy versus the Richardson number compared to the EFB theory. All state-of-the-art turbulent closures follow the so-called Kolmogorov hypothesis: all dissipation timescales of turbulent second-order moments are assumed to be proportional to each other, which, at first glance, looks reasonable but is, in fact, hypothetical for stably stratified turbulence.

The present study aims to demonstrate the dependence of dissipation timescales of basic second-order moments on stability through direct numerical simulation (DNS) experiments. The obtained numerical results allow us to modify the EFB turbulence closure theory to account for that dependency. It is worth noting that the DNSs presented here are limited to bulk Richardson numbers (based on the wall velocity and temperature differences and channel height) up to *Ri*_{b}=0.11 and Reynolds numbers (based on the wall velocity difference and channel height; see Sect. 3) up to *Re*=120 000.

This paper is organised as follows. In Sect. 2, we formulate basic budget equations and main assumptions in the framework of the EFB turbulence closure theory. Section 3 describes the setup for DNSs of stably stratified turbulent plane Couette flow to determine the vertical profiles of the dissipation timescales of turbulent second-order moments. In Sect. 4, we formulate the modified EFB turbulence closure theory while considering the dependencies of the dissipation timescales of basic second-order moments on the gradient Richardson number obtained from DNSs. There, we also perform the validation of the modified EFB turbulence closure model, which yields vertical profiles of the basic turbulence parameters (including the turbulent Prandtl number, the ratio of TPE to TKE, and the normalised turbulent heat flux) using the data from the DNS. Finally, in Sect. 5, we discuss the obtained results and draw the conclusions.

We consider plane-parallel, stably stratified dry-air flow and employ the familiar budget equations underlying turbulence closure theory (e.g. Kleeorin et al., 2021; Zilitinkevich et al., 2013, 2019; Kaimal and Finnigan, 1994; Canuto et al., 2008) for the Reynolds stress, *τ*_{ij}=〈*u*_{i}*u*_{j}〉; the turbulent flux of potential temperature, *F*_{i}=〈*θ**u*_{i}〉; and the intensity of potential temperature fluctuations, ${E}_{\mathit{\theta}}=\langle {\mathit{\theta}}^{\mathrm{2}}\rangle /\mathrm{2}$:

Here, *x*_{1}=*x* and *x*_{2}=*y* are horizontal coordinates; *x*_{3}=*z* is the vertical coordinate; *t* is the time; $\mathbf{U}=\left({U}_{\mathrm{1}},{U}_{\mathrm{2}},{U}_{\mathrm{3}}\right)=\left(U,V,W\right)$ is the mean flow velocity; $\mathbf{u}=\left({u}_{\mathrm{1}},{u}_{\mathrm{2}},{u}_{\mathrm{3}}\right)=\left(u,v,w\right)$ is the velocity fluctuations; $\mathrm{\Theta}=T{\left({P}_{\mathrm{0}}/P\right)}^{\mathrm{1}-\mathrm{1}/\mathit{\gamma}}$ is the mean potential temperature (expressed through absolute temperature, *T*, and pressure, *P*); *T*_{0}, *P*_{0}, and *ρ*_{0} are reference values of temperature, pressure, and density, respectively; $\mathit{\gamma}={c}_{p}/{c}_{v}=\mathrm{1.41}$ is the ratio of specific heats; *θ* and *p* are fluctuations of potential temperature and pressure; $\mathrm{D}/\mathrm{D}t=\partial /\partial t+{U}_{k}\partial /\partial {x}_{k}$ is the advective derivative; angle brackets denote the averaging; $\mathit{\beta}=g/{T}_{\mathrm{0}}$ is the buoyancy parameter; *g* is the acceleration due to gravity; *δ*_{ij} is the unit tensor (*δ*_{ij}=1 for *i*=*j* and *δ*_{ij}=0 for *i*≠*j*); and ${\mathrm{\Phi}}_{ij\mathrm{3}}^{\left(\mathit{\tau}\right)}$, ${\mathrm{\Phi}}_{i}^{\left(F\right)}$, and Φ^{(θ)} are the third-order moments, which describe turbulent transport of the second-order moments under consideration.

where *Q*_{ij} is the correlation between fluctuations of pressure and strain-rate tensor, which controls the interactions between the Reynolds stress components as follows:

Here, ${\mathit{\epsilon}}_{ij}^{\left(\mathit{\tau}\right)}$, ${\mathit{\epsilon}}_{i}^{\left(F\right)}$, and *ε*_{θ} are the dissipation rates of the second-order moments:

where *ν* is kinematic viscosity and *κ* is thermal conductivity.

The budget of TKE components, ${E}_{i}=\u2329{u}_{i}^{\mathrm{2}}\u232a/\mathrm{2}$ ($i=\mathrm{1},\mathrm{2},\mathrm{3}$), is determined by Eq. (1) for *i*=*j*, which yields the familiar TKE budget equation:

where ${E}_{K}=\sum {E}_{i}$ is TKE and ${\mathit{\epsilon}}_{K}=\sum {\mathit{\epsilon}}_{ii}^{\left(\mathit{\tau}\right)}/\mathrm{2}$ is the TKE dissipation rate. The sum of the term *Q*_{ii} (the trace of the tensor *Q*_{ij}) is equal to zero because of the incompressibility, constraint on the flow velocity field, $\partial {u}_{i}/\partial {x}_{i}=\mathrm{0}$; that is, *Q*_{ij} only redistributes energy between TKE components.

Likewise, *ε*_{θ} is the dissipation rate of the intensity of potential temperature fluctuations, *E*_{θ}, and ${\mathit{\epsilon}}_{i}^{\left(F\right)}$ are the dissipation rates of the three components of the turbulent flux of potential temperature, *F*_{i}.

Following Kolmogorov (1941b, 1942), the dissipation rates, *ε*_{K} and *ε*_{θ}, are taken to be proportional to the dissipating quantities divided by corresponding timescales:

where *t*_{K} is the TKE dissipation timescale and *t*_{θ} is the dissipation timescale of *E*_{θ}. Here, the formulation of the dissipation rates is not hypothetical: it merely expresses one unknown (dissipation rate) through another (dissipation timescale).

In this study, we consider the EFB theory in its simplest, algebraic form, neglecting non-steady terms in all budget equations and neglecting divergence of the fluxes of TKE and TPE (determined by third-order moments). This approach is reasonable because, for example, the characteristic times of variations of the second moments are much larger than the turbulent timescales for large Reynolds and Péclet numbers. We also assume that the terms related to the divergence of the fluxes of TKE and TPE for stably stratified turbulence are much smaller than the rates of production and dissipation in budget Eqs. (3) and (11). In this case, the TKE budget equation, Eq. (11), and the budget equation for *E*_{θ}, Eq. (3), become

The intensity of the potential temperature fluctuations, *E*_{θ}, determines TPE as follows:

so that Eq. (14) becomes

where ${\mathit{\epsilon}}_{P}={E}_{P}/{t}_{\mathit{\theta}}$ is the TPE dissipation time.

The first term on the right-hand side of Eq. (13), $-\mathit{\tau}\partial U/\partial z$, is the rate of the TKE production, while the second term, *β**F*_{z}, is the buoyancy, which, in stably stratified flow, causes decay of TKE; that is, it results in the conversion of TKE into TPE. The ratio of these terms is the flux Richardson number:

and this dimensionless parameter characterises the effect of stratification on turbulence.

Taking into account Eq. (17), the steady-state versions of TKE and TPE budget equations, Eqs. (13) and (14), can be rewritten as

Thus, the ratio of TPE to TKE is

Zilitinkevich et al. (2013) suggested the following relation, linking *R**i*_{f} with another stratification parameter, $z/L$:

where $L=-{\mathit{\tau}}^{\mathrm{3}/\mathrm{2}}/\mathit{\beta}{F}_{z}$ is the Obukhov length scale, *k*=0.4 is the von Kármán constant, and *R*_{∞}=0.2 is the maximum value of the flux Richardson number.

On the right-hand side of Eq. (20), there is an unknown ratio between two dissipation timescales, ${t}_{\mathit{\theta}}/{t}_{K}$. The Kolmogorov hypothesis suggests that it is a universal constant. We do not imply that this assumption applies but instead investigate a possible stability dependency of dissipation timescale ratios and improve the EFB turbulence closure model accounting for it. To this end, we perform DNSs of stably stratified turbulent plane Couette flow (see Sect. 3) to measure the dissipation timescales of basic second-order moments and validate the modified EFB turbulence closure model (see Sect. 4).

For our study, we conducted a series of direct numerical simulations of stably stratified turbulent plane Couette flow. This flow occurs between two parallel plates that move relative to each other, producing shear and turbulence, with the plates at different temperatures, thus creating stable stratification. In Couette flow, the total (turbulent plus molecular) vertical fluxes of momentum and potential temperature remain constant, independent of distance from the walls, which, in particular, ensures a very certain fixed value of the Obukhov length scale. Figure 1 illustrates the profiles of mean flow velocity and mean potential temperature. We recall that all our derivations are relevant to the well-developed turbulence regime where molecular transports are negligible compared to turbulent transports, so turbulent fluxes practically coincide with total fluxes. This is the case in our DNS, except for the narrow near-wall viscous turbulent-flow transition layers. Data from these layers, obviously irrelevant to the turbulence regime we consider, are shown by grey points in the figures and are ignored in fitting procedures. In further analysis, we primarily utilise $z/L$ as a stratification parameter instead of *Ri* or *Ri*_{f} because it offers a better dynamic range in our experiments. While *Ri* remains practically constant in each DNS run and *Ri*_{f} is limited in its growth, the parameter $z/L$ is determined by the distance from the walls, thus varying significantly in every DNS run.

Numerical simulation of stably stratified turbulent Couette flow was performed using the unified DNS, LES, and Reynolds-averaged Navier–Stokes equations (RANS) code developed at the Moscow State University (MSU) and the Institute of Numerical Mathematics (INM) of the Russian Academy of Science (see Mortikov, 2016; Mortikov et al., 2019; Bhattacharjee et al., 2022; Debolskiy et al., 2023; Gladskikh et al., 2023; Zasko et al., 2023). The code is designed for high-resolution simulations on modern-day high-performance computing (HPC) systems. The DNS part of the code solves the finite-difference approximation of the incompressible Navier–Stokes system of equations under the Boussinesq approximation. Conservative schemes on the staggered grid (Morinishi et al., 1998; Vasilyev, 2000) of fourth-order accuracy are used in the horizontal direction, while, in the vertical direction, the spatial approximation is restricted to second-order accuracy with near-wall grid resolution refinement sufficient to resolve the near-wall viscous region. The time step used in the simulations was determined by Courant–Friedrichs–Lewy (CFL) restrictions, with CFL maintained at approximately 0.1 in all runs. This corresponds to a value of ${u}_{\ast}^{\mathrm{2}}\mathrm{\Delta}t/\mathit{\nu}$ on the order of 0.01. The projection method (Brown et al., 2001) is used for the time advancement of momentum equations coupled with the incompressibility condition, while the multigrid method is applied to solve the Poisson equation to ensure that the velocity is divergence-free at each time step. For the Couette flow, periodic boundary conditions are used in the horizontal directions, and no-slip and no-penetration conditions are set on the channel walls for the velocity. The stable stratification is maintained by prescribed Dirichlet boundary conditions on the potential temperature. In all experiments, the value of molecular Prandtl number (ratio of kinematic viscosity and thermal diffusivity of the fluid) was fixed at 0.7 based on its typical value for air (Monin and Yaglom, 1971). The simulations were performed for a wide range of the Reynolds number, *Re*, defined by the wall velocity difference, channel height, and kinematic viscosity, that is from 40 000 up to 120 000 (see Table 1). All experiments were carried out using the resources of MSU and Center for Scientific Computing (CSC) HPC facilities. For the maximum *Re* values achieved, the numerical grid consisted of more than 2×10^{8} cells and the calculations used about 10 000 CPU cores.

For each Reynolds number, we conducted a series of experiments. Beginning with neutral conditions (no imposed gradient of the mean potential temperature), we incrementally increased the bulk Richardson number, which characterises the stable stratification, in each successive experiment. By gradually increasing the stability in each experiment, we were able to cover a wide range of *Ri* values, extending from neutral to stably stratified states. In each run, the turbulent flow was allowed sufficient time to develop and reach statistical steady-state conditions, which required a spin-up period of at least 15 $H/{u}_{\ast}$ periods. This ensured that parameters such as the total momentum flux remained constant and that the TKE balance was in a steady state. The fully developed steady state was used as the initial conditions for the higher *Ri* or *Re* experiment setups. Additionally, all terms in the second-order moment budget equations (Eqs. 1–3) were evaluated consistently using the finite-difference approximation used, resulting in a negligible residual. This approach enabled us to comprehensively study the characteristics of shear-produced stably stratified turbulence, explicitly resolving all dissipation timescales of turbulent second-order moments.

In the steady-state, Eq. (1) for the vertical component of the turbulent flux of momentum, *τ*, becomes

Following Zilitinkevich et al. (2007, 2013, 2019) we define the sum of all terms in square brackets on the right-hand side of Eq. (22) as the “effective dissipation”:

Thus, Eq. (22) becomes

yielding the well-known down-gradient formulation of the vertical turbulent flux of momentum:

where ${A}_{z}\equiv {E}_{z}/{E}_{K}$ is the vertical share of TKE (the vertical anisotropy parameter). By substituting Eq. (25) into Eq. (18), we obtain

In Eq. (26), all the variables are exactly resolved numerically in DNSs, making a detailed investigation on ${t}_{\mathit{\tau}}/{t}_{K}$ possible. Figure 2 demonstrates that the dissipation timescale ratio, ${t}_{\mathit{\tau}}/{t}_{K}$, is a function of the stratification parameter $z/L$ rather than a constant. We propose to approximate this function with a ratio of two first-order polynomials:

Here, the dimensionless empirical constants are obtained from the best fit of Eq. (27) to DNS bin-averaged data: ${C}_{\mathrm{1}}^{\mathit{\tau}K}=\mathrm{0.08}$, ${C}_{\mathrm{2}}^{\mathit{\tau}K}=\mathrm{0.4}$, and ${C}_{\mathrm{3}}^{\mathit{\tau}K}=\mathrm{2}$. The fitting is done using the rational regression model of Curve Fitting Toolbox version 3.5.13 (R2021a) (2022). The ratio of two first-order polynomials is chosen as a simpler fitting function that could provide monotonicity, reasonable smoothness, and clear asymptotes. The only three adjustable parameters of this approximation correspond to the function value at $z/L=\mathrm{0}$, the $z/L\to \mathrm{\infty}$ limit, and the transition between them.

Proceeding to the vertical flux of potential temperature, *F*_{z}, we derive its steady-state budget equation from Eq. (2):

DNS modelling has shown that $\frac{\partial}{\partial z}{\mathrm{\Phi}}_{z}^{\left(F\right)}$ term is of the same order of magnitude as *ε*_{F} and of the same sign, so we introduce the “effective dissipation rate”, ${\mathit{\epsilon}}_{F}^{\left(\mathrm{eff}\right)}$:

Consequently, Eq. (28) reduces to

Traditionally, the pressure term was either assumed to be negligible or declared to be proportional to the *β*〈*θ*^{2}〉 term (see Zilitinkevich et al., 2007, 2013). However, our DNS data have shown that it is neither negligible nor proportional to any other term in the budget equation, Eq. (30). Instead, we found that it is well approximated by a linear combination of the production and transport terms of Eq. (30) (see Fig. 3):

The dimensionless constants *C*_{θ}=0.82 and ${C}_{\mathrm{\nabla}}=-\mathrm{0.80}$ are obtained from the best fit of Eq. (31) to DNS data.

By substituting Eq. (31) into Eq. (30), we rewrite the budget equation as

Substituting Eq. (15) for 〈*θ*^{2}〉 into Eq. (32) allows us to express *F*_{z} by a familiar temperature gradient expression:

Substituting Eq. (33) into Eq. (14) gives

Next, the turbulent Prandtl number, defined as ${\mathit{Pr}}_{T}={K}_{M}/{K}_{H}$, is given by

Equations (34) and (35) provide us with two constrains on the function in the square brackets. First, the left-hand side of Eq. (34) is non-negative by definition, implying that the same requirement applies for the right-hand side of the equation. Second, the turbulent Prandtl number grows with the increase of the gradient Richardson number – that is, ${\left({\mathit{Pr}}_{T}\right|}_{(z/L\to \mathrm{\infty})}\to \mathit{Ri}/{R}_{\mathrm{\infty}}$ – requiring the function in the square brackets to approach zero under extreme stratification. This leads us to the following approximation (see Fig. 4):

This function monotonically decreases from 1 to $\mathrm{0}<z/L<\mathrm{\infty}$, satisfying our requirements with *C*_{Pr}=0.65. The observed spread of data points might be explained by the simulation time being insufficient to reach a fully statistical steady state for this specific ratio. Although the fully developed steady state was achieved (verified using the standard criterion of stabilised TKE, which showed no significant fluctuations over time), the parameter involving ratios of temperature fluctuations, *θ*, might require additional time to stabilise. We believe that increasing the experiment time would decrease the spread, but we leave the validation of this hypothesis for future studies.

This leads us to a similar approximation of ${t}_{\mathit{\tau}}/{t}_{F}$ (see Fig. 5):

Now, to complete the closure, we need to determine one more dimensionless ratio, ${t}_{\mathit{\theta}}/{t}_{K}$. It is explicitly required for the ratio of turbulent energies, ${E}_{P}/{E}_{K}$, and consequently for *A*_{z} through Eqs. (20) and (36). We approximate it once again with the ratio of two first-order polynomials:

With Eq. (38), our turbulence closure is now complete, allowing us to proceed with the verification process using quantities not utilised in the fitting procedures. Figure 7 provides empirical evidence that supports the stability dependencies given by Eqs. (20), (26), (27), and (34)–(38). Table 2 summarises the proposed approximations and provides a summary of the resulting turbulent closure.

For practical reasons, most operational numerical weather prediction and climate models parameterise these dimensionless ratios as functions of the gradient Richardson number rather than $z/L$. This preference arises from the fact that the gradient Richardson number is defined solely by mean quantities – namely, buoyancy and shear productions – which in practice imposes fewer computational restrictions on the model's time step. Since *Ri*=*Pr*_{T}*Ri*_{f} and both *Pr*_{T} and *Ri*_{f} are defined as functions of $z/L$ by Eqs. (35) and (21), respectively, we can derive an expression for the gradient Richardson number, *Ri*, as the function of $z/L$, as shown in Fig. 8:

For many years, our understanding of dissipation rates for turbulent second-order moments has been hindered by a lack of direct observations in fully controlled conditions, particularly in a strongly stable stratification. To address this limitation, we conducted topical DNS experiments of stably stratified Couette flows. The main finding of this study is that, contrary to the traditional assumption of them being proportional to a single universal dissipation timescale, the ratios of the dissipation timescales of the basic second-order moments depend on the temperature stratification (e.g. characterised by the gradient Richardson number).

This finding laid the foundation for empirically approximating these ratios with simple universal functions of stability parameters, valid for a wide range of stratifications. Consequently, this allowed us to refine the EFB turbulent closure by accounting for dissipation timescales that are intrinsic to the basic second-order moments. As a result, the revised formulations for eddy viscosity and eddy conductivity reveal greater physical consistency in stratified conditions, thereby enhancing the representation of turbulence in numerical weather prediction and climate modelling.

We have also observed that the dimensionless parameters involving *θ* fluctuations demonstrate a wider spread of values within and across the DNS experiments, making it more challenging to approximate them with stability functions. This suggests that the stabilisation time for these parameters may be significantly longer than for TKE components.

It is important to note that our DNS experiments were limited to gradient Richardson numbers up to *Ri*=0.17. Any data reliably indicating different asymptotic values of the timescale dimensionless ratios or demonstrating their different dependency on the temperature stratification would impose the need to readjust the proposed parameterisation.

We deliberately avoided discussing intermittency issues; for that, one needs to determine higher-order two-point (or multi-point) moments. Intermittency is important for small-scale effects, and intermittency implies that higher-order moments of velocity and temperature fields have non-Gaussian statistics. In this study, we focused on larger scales, determining one-point second-order correlation functions that barely touch one-point third-order correlation functions only when it is necessary. However, addressing this topic would be crucial for advancing numerical simulations towards higher stratifications and warrants detailed investigation.

With these considerations in mind, we believe that the most challenging step will be to explicitly explore the transitional region between traditional weakly stratified turbulence and extremely stable stratification, where the behaviour of the turbulent Prandtl number shifts from nearly constant to a linear function with respect to the gradient Richardson number. Investigating this phenomenon would require unprecedented computational resources for DNS or specialised in situ or laboratory experiments.

The DNS code is available on GitLab at http://tesla.parallel.ru/emortikov/nse-couette-dns (Mortikov et al., 2024). The access can be granted by Evgeny Mortikov (mortikov@srcc.msu.ru). The datasets generated and analysed during the current study are available at https://doi.org/10.23728/b2share.7a1d875b872748c7bf566ece352c0a10 (Kadantsev and Mortikov, 2024).

EK conceptualised the paper, performed data analysis, wrote the initial text, and prepared the figures. EM contributed to the conceptualisation of the study, developed the DNS code, and performed the numerical simulations. AG contributed to the conceptualisation of the study and code development. NK and IR contributed to the conceptualisation of the study and assisted with literature overview and manuscript editing.

The contact author has declared that none of the authors has any competing interests.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

This paper was not only inspired by but also conducted under the supervision of the esteemed researcher Sergej Zilitinkevich, who unfortunately is no longer with us. We wish to express our profound gratitude to Sergej for the incredible honour of collaborating with him and for the immense inspiration he generously bestowed upon us. The authors wish to acknowledge CSC – IT Center for Science, Finland, and the shared research facilities at MSU for computational resources.

This research has been supported by the EU Horizon 2020 (grant no. 101036245), the Research Council of Finland (grant nos. 334798 and 337552), the FSTP project no. 124042700008-6 titled “Research in geophysical boundary layers and the development of new modelling approaches for Earth system models” within the programme “Improvement of the global world-level Earth system model for research purposes and scenarios forecasting of climate change”, the RSF grant (agreement no. 21-71-30003; development of the DNS model), and the Ministry of Education and Science of the Russian Federation as part of the programme of the Moscow Center of Fundamental and Applied Mathematics (agreement no. 075-15-2022-284; DNS of stably stratified Couette flow).

Open-access funding was provided by the Helsinki University Library.

This paper was edited by Christian Franzke and reviewed by three anonymous referees.

Batchelor, G. K.: The Theory of Homogeneous Turbulence, Cambridge, Cambridge University Press, Cambridge, 1953.

Bhattacharjee, S., Mortikov, E. V., Debolskiy, A. V., Kadantsev, E., Pandit., R., Vesala, T., and Sahoo, G.: Direct Numerical Simulation of a Turbulent Channel Flow with Forchheimer Drag, Bound.-Lay. Meteorol., 185, 259–276, https://doi.org/10.1007/s10546-022-00731-8, 2022.

Brown, D. L., Cortez, R., and Minion, M. L.: Accurate projection methods for the incompressible Navier–Stokes equations, J. Comp. Phys., 168, 464–499, 2001.

Canuto, V. and Minotti, F.: Stratified turbulence in the atmosphere and oceans: A new subgrid model, J. Atmos. Sci., 50, 1925–1935, https://doi.org/10.1175/1520-0469(1993)050<1925:STITAA>2.0.CO;2, 1993.

Canuto, V., Howard, A., Cheng, Y., and Dubovikov, M.: Ocean turbulence, part I: One-point closure model—Momentum and heat vertical diffusivities, J. Phys. Oceanogr., 31, 1413–1426, https://doi.org/10.1175/1520-0485(2001)031<1413:OTPIOP>2.0.CO;2, 2001.

Canuto, V., Cheng, Y., Howard, A. M., and Esau, I.: Stably stratified flows: a model with no Ri(cr), J. Atmos. Sci., 65, 2437–2447, 2008.

Cheng, Y., Canuto, V., and Howard, A. M.: An improved model for the turbulent PBL, J. Atmos. Sci., 59, 1550–1565, 2002.

Curve Fitting Toolbox version: 3.5.13 (R2021a), Natick, The MathWorks Inc., Massachusetts, https://www.mathworks.com/products/curvefitting.html (last access: 15 September 2024), 2022.

Dalaudier, F. and Sidi, C.: Evidence and Interpretation of a Spectral Gap in the Turbulent Atmospheric Temperature Spectra, J. Atmos. Sci., 44, 3121–3126, https://doi.org/10.1175/1520-0469(1987)044<3121:EAIOAS>2.0.CO;2, 1987.

Davidson, P. A.: Turbulence in Rotating, Stratified and Electrically Conducting Fluids, Cambridge University Press, Cambridge, https://doi.org/10.1017/CBO9781139208673, 2013.

Debolskiy, A. V., Mortikov, E. V., Glazunov, A. V., and Lüpkes, C.: Evaluation of Surface Layer Stability Functions and Their Extension to First Order Turbulent Closures for Weakly and Strongly Stratified Stable Boundary Layer, Bound.-Lay. Meteorol., 187, 73–93, https://doi.org/10.1007/s10546-023-00784-3, 2023.

Elperin, T., Kleeorin, N., Rogachevskii, I., and Zilitinkevich, S.: Formation of large-scale semi-organized structures in turbulent convection, Phys. Rev. E, 66, 066305, https://doi.org/10.1103/PhysRevE.66.066305, 2002.

Elperin, T., Kleeorin, N., Rogachevskii, I., and Zilitinkevich, S.: Tangling turbulence and semi-organized structures in convective boundary layers, Bound.-Lay. Meteorol., 119, 449, https://doi.org/10.1007/s10546-005-9041-5, 2006.

Frisch, U.: Turbulence: the Legacy of A. N. Kolmogorov, Cambridge University Press, Cambridge, https://doi.org/10.1017/CBO9781139170666, 1995.

Gladskikh, D., Ostrovsky, L., Troitskaya, Yu., Soustova, I., and Mortikov, E.: Turbulent Transport in a Stratified Shear Flow, J. Mar. Sci. Eng., 11, 136, https://doi.org/10.3390/jmse11010136, 2023.

Hanazaki, H. and Hunt, J.: Linear processes in unsteady stably stratified turbulence, J. Fluid Mech., 318, 303–337, https://doi.org/10.1017/S0022112096007136, 1996.

Hanazaki, H. and Hunt, J.: Structure of unsteady stably stratified turbulence with mean shear, J. Fluid Mech., 507, 1–42, https://doi.org/10.1017/S0022112004007888, 2004.

Holloway, G.: Estimation of oceanic eddy transports from satellite altimetry, Nature, 323, 243–244, https://doi.org/10.1038/323243a0, 1986.

Hunt, J., Wray, A., and Moin, P.: Eddies, Stream, and Convergence Zones in Turbulent Flows, Proceeding of the Summer Program in Center for Turbulence Research, 193–208, 1988.

Kadantsev, E. and Mortikov, E.: Direct Numerical Simulations of stably stratified turbulent plane Couette flow, B2Share [data set], https://doi.org/10.23728/B2SHARE.7A1D875B872748C7BF566 ECE352C0A10, 2024.

Kaimal, J. C. and Finnigan, J. J.: Atmospheric boundary layer flows, Oxford University Press, New York, 289 pp., 1994.

Keller, K. and van Atta, C.: An experimental investigation of the vertical temperature structure of homogeneous stratified shear turbulence, J. Fluid Mech., 425, 1–29, https://doi.org/10.1017/S0022112000002111, 2000.

Kleeorin, N., Rogachevskii, I., Soustova, I. A., Troitskaya, Y. I., Ermakova, O. S., and Zilitinkevich S.: Internal gravity waves in the energy and flux budget turbulence closure theory for shear-free stably stratified flows, Phys. Rev. E, 99, 063106, https://doi.org/10.1103/PhysRevE.99.063106, 2019.

Kleeorin, N., Rogachevskii, I., and Zilitinkevich, S.: Energy and flux budget closure theory for passive scalar in stably stratified turbulence, Phys. Fluids, 33, 076601, https://doi.org/10.1063/5.0052786, 2021.

Kolmogorov, A. N.: Dissipation of energy in the locally isotropic turbulence, Dokl. Akad. Nauk SSSR A, 32, 16, 1941a.

Kolmogorov, A. N.: Energy dissipation in locally isotropic turbulence, Dokl. Akad. Nauk. SSSR A, 32, 19, 1941b.

Kolmogorov, A. N.: The equations of turbulent motion in an incompressible fluid, Izvestia Akad. Sci., USSR; Phys., 6, 56, 1942.

Kolmogorov, A. N.: The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers, P. Roy. Soc. Lond. A, 434, 9, https://www.jstor.org/stable/51980 (last access: 18 September 2024), 1991.

L'vov, V. S., Procaccia, I., and Rudenko, O.: Turbulent fluxes in stably stratified boundary layers, Phys. Scr., 132, 1–15, 2008.

Li, D., Katul, G., and Zilitinkevich, S.: Closure Schemes for Stably Stratified Atmospheric Flows without Turbulence Cutoff, J. Atmos. Sci., 73, 4817–4832, https://doi.org/10.1175/JAS-D-16-0101.1, 2016.

Mahrt, L.: Stably stratified atmospheric boundary layers, Annu. Rev. Fluid Mech., 46, 23, https://doi.org/10.1146/annurev-fluid-010313-141354, 2014.

Mauritsen, T., Svensson, G., Zilitinkevich, S., Esau, I., Enger, L., and Grisogono, B.: A total turbulent energy closure model for neutrally and stably stratified atmospheric boundary layers, J. Atmos. Sci., 64, 4117–4130, 2007.

Monin, A. S. and Yaglom, A. M.: Statistical Fluid Mechanics, vol. 1, MIT Press, Cambridge, ISBN-10 0262130629, ISBN-13 978-0262130622, 1971.

Monin, A. S. and Yaglom, A. M.: Statistical Fluid Mechanics, vol. 2, Courier Corporation, ISBN-10 0486318141, ISBN-13 978-0486318141, 2013.

Morinishi, Y., Lund, T. S., Vasilyev, O. V., and Moin, P.: Fully conservative higher order finite difference schemes for incompressible flows, J. Comp. Phys., 143, 90–124, 1998.

Mortikov, E. V.: Numerical simulation of the motion of an ice keel in a stratified flow, Izv. – Atmos. Ocean. Phys., 52, 108–115, https://doi.org/10.1134/S0001433816010072, 2016.

Mortikov, E. V., Glazunov, A. V., and Lykosov V. N.: Numerical study of plane Couette flow: turbulence statistics and the structure of pressure-strain correlations, Russ. J. Numer. Anal. Math. Model., 34, 119–132, https://doi.org/10.1515/rnam-2019-0010, 2019.

Mortikov, E., Glazunov, A., Debolskiy, A., and Gashuk, E.: nse-couette-dns, GitLab [code], http://tesla.parallel.ru/emortikov/nse-couette-dns, last access: 15 September 2024.

Ostrovsky, L, and Troitskaya, Yu.: A model of turbulent transfer and dynamics of turbulence in a stratified shear flow, Izvestiya AN SSSR FAO, 23, 1031–1040, 1987.

Pope, S. B.: Turbulent Flows, Cambridge University Press, Cambridge, https://doi.org/10.1017/CBO9780511840531, 2000.

Rehmann, C. R. and Hwang, J. H.: Small-Scale Structure of Strongly Stratified Turbulence, J. Phys. Oceanogr., 35, 151–164, 2005.

Rogachevskii, I.: Introduction to Turbulent Transport of Particles, Temperature and Magnetic Fields, Cambridge University Press, Cambridge, https://doi.org/10.1017/9781009000918, 2021.

Rogachevskii, I. and Kleeorin, N.: Semi-organised structures and turbulence in the atmospheric convection, Phys. Fluids, 36, 026610, https://doi.org/10.1063/5.0188732, 2024.

Rogachevskii, I., Kleeorin, N., and Zilitinkevich, S.: The energy- and flux budget theory for surface layers in atmospheric convective turbulence, Phys. Fluids, 34, 116602, https://doi.org/10.1063/5.0123401, 2022.

Schumann, U. and Gerz, T.: Turbulent mixing in stably stratified shear flows, J. Appl. Meteorol., 34, 33–48, 1995.

Stretch, D. D., Rottman, J. W., Nomura, K. K., and Venayagamoorthy, S. K.: Transient mixing events in stably stratified turbulence, in: 14th Australasian Fluid Mechanics Conference, 10–14 December 2001, Adelaide, Australia, 2001.

Sukoriansky, S. and Galperin, B.: Anisotropic turbulence and internal waves in stably stratified flows (QNSE theory), Phys. Scr., 132, 1–8, 2008.

Tennekes, H. and Lumley, J. L.: A First Course in Turbulence, MIT Press, Cambridge, ISBN-10 0262200198, ISBN-13 978-0262200196, 1972.

Umlauf, L.: Modelling the effects of horizontal and vertical shear in stratified turbulent flows, Deep-Sea Res. Pt. II, 52, 1181–1201, 2005.

Umlauf, L. and Burchard, H.: Second-order turbulence closure models for geophysical boundary layers. A review of recent work, Cont. Shelf Res., 25, 795, https://doi.org/10.1016/j.csr.2004.08.004, 2005.

Vasilyev, O. V.: High order finite difference schemes on non-uniform meshes with good conservation properties, J. Comp. Phys., 157, 746–761, 2000.

Weng, W. and Taylor, P.: On modelling the one-dimensional Atmospheric Boundary Layer, Bound.-Lay. Meteorol., 107, 371–400, 2003.

Zasko, G., Glazunov, A., Mortikov, E., Nechepurenko, Yu., and Perezhogin, P.: Optimal Energy Growth in Stably Stratified Turbulent Couette Flow, Bound.-Lay. Meteorol., 187, 395–421, https://doi.org/10.1007/s10546-022-00744-3, 2023.

Zilitinkevich, S., Elperin, T., Kleeorin, N., and Rogachevskii, I.: Energy- and flux budget (EFB) turbulence closure model for stably stratified flows. I: Steady-state, homogeneous regimes, Bound.-Lay. Meteorol., 125, 167, https://doi.org/10.1007/s10546-007-9189-2, 2007.

Zilitinkevich, S., Elperin, T., Kleeorin, N., Rogachevskii, I., Esau, I., Mauritsen, T., and Miles, M. W.: Turbulence energetics in stably stratified geophysical flows: Strong and weak mixing regimes, Q. J. Roy. Meteor. Soc., 134, 793, https://doi.org/10.1002/qj.264, 2008.

Zilitinkevich, S., Elperin, T., Kleeorin, N., L'vov, V., and Rogachevskii, I.: Energy-and flux-budget turbulence closure model for stably stratified flows. II: The role of internal gravity waves, Bound.-Lay. Meteorol., 133, 139, https://doi.org/10.1007/s10546-009-9424-0, 2009.

Zilitinkevich, S., Elperin, T., Kleeorin, N., Rogachevskii, I., and Esau, I.: A hierarchy of energy- and flux budget (EFB) turbulence closure models for stably stratified geophysical flows, Bound.-Lay. Meteorol., 146, 341, https://doi.org/10.1007/s10546-012-9768-8, 2013.

Zilitinkevich, S., Druzhinin, O., Glazunov, A., Kadantsev, E., Mortikov, E., Repina, I., and Troitskaya, Y.: Dissipation rate of turbulent kinetic energy in stably stratified sheared flows, Atmos. Chem. Phys., 19, 2489–2496, https://doi.org/10.5194/acp-19-2489-2019, 2019.