No Arabic abstract
We have advanced the energy and flux budget (EFB) turbulence closure theory that takes into account a two-way coupling between internal gravity waves (IGW) and the shear-free stably stratified turbulence. This theory is based on the budget equation for the total (kinetic plus potential) energy of IGW, the budget equations for the kinetic and potential energies of fluid turbulence, and turbulent fluxes of potential temperature for waves and fluid flow. The waves emitted at a certain level, propagate upward, and the losses of wave energy cause the production of turbulence energy. We demonstrate that due to the nonlinear effects more intensive waves produce more strong turbulence, and this, in turns, results in strong damping of IGW. As a result, the penetration length of more intensive waves is shorter than that of less intensive IGW. The anisotropy of the turbulence produced by less intensive IGW is stronger than that caused by more intensive waves. The low amplitude IGW produce turbulence consisting up to 90 % of turbulent potential energy. This resembles the properties of the observed high altitude tropospheric strongly anisotropic (nearly two-dimensional) turbulence.
The energy and flux budget (EFB) closure theory for a passive scalar (non-buoyant and non-inertial particles or gaseous admixtures) is developed for stably stratified turbulence. The physical background of the EFB turbulence closures is based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, as well as the turbulent flux of particles. The EFB turbulence closure is designed for stratified geophysical flows from neutral to very stable stratification and it implies that turbulence is maintained by the velocity shear at any stratification. In a steady-state, expressions for the turbulent flux of passive scalar and the anisotropic non-symmetric turbulent diffusion tensor are derived, and universal flux Richardson number dependencies of the components of this tensor are obtained. The diagonal component in the vertical direction of the turbulent diffusion tensor is suppressed by strong stratification, while the diagonal components in the horizontal directions are not suppressed, and they are dominant in comparison with the other components of turbulent diffusion tensor. This implies that any initially created strongly inhomogeneous particle cloud is evolved into a thin pancake in horizontal plane with very slow increase of its thickness in the vertical direction. The turbulent Schmidt number increases linearly with the gradient Richardson number. Considering the applications of these results to the atmospheric boundary-layer turbulence, the theoretical relationships are derived which allow to determine the turbulent diffusion tensor as a function of the vertical coordinate measured in the units of the local Obukhov length scale. The obtained relations are potentially useful in modelling applications of particle dispersion in the atmospheric boundary-layer turbulence and free atmosphere turbulence.
In this paper we advance physical background of the energy- and flux-budget turbulence closure based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time-scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth rotation. In accordance to modern experimental evidence, the closure implies maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: strong turbulence at Ri << 1 typical of boundary-layer flows and characterised by the practically constant turbulent Prandtl number; and weak turbulence at Ri > 1 typical of the free atmosphere or deep ocean, where the turbulent Prandtl number asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy-viscosity and eddy-conductivity, and general non-gradient model based on prognostic equations for all basic parameters of turbulence including turbulent fluxes.
We find and investigate via numerical simulations self-sustained two-dimensional turbulence in a magnetohydrodynamic flow with a maximally simple configuration: plane, noninflectional (with a constant shear of velocity) and threaded by a parallel uniform background magnetic field. This flow is spectrally stable, so the turbulence is subcritical by nature and hence it can be energetically supported just by transient growth mechanism due to shear flow nonnormality. This mechanism appears to be essentially anisotropic in spectral (wavenumber) plane and operates mainly for spatial Fourier harmonics with streamwise wavenumbers less than a ratio of flow shear to the Alfv{e}n speed, $k_y < S/u_A$ (i.e., the Alfv{e}n frequency is lower than the shear rate). We focused on the analysis of the character of nonlinear processes and underlying self-sustaining scheme of the turbulence, i.e., on the interplay between linear transient growth and nonlinear processes, in spectral plane. Our study, being concerned with a new type of the energy-injecting process for turbulence -- the transient growth, represents an alternative to the main trends of MHD turbulence research. We find similarity of the nonlinear dynamics to the related dynamics in hydrodynamic flows -- to the emph{bypass} concept of subcritical turbulence. The essence of the analyzed nonlinear MHD processes appears to be a transverse redistribution of kinetic and magnetic spectral energies in wavenumber plane [as occurs in the related hydrodynamic flow, see Horton et al., Phys. Rev. E {bf 81}, 066304 (2010)] and differs fundamentally from the existing concepts of (anisotropic direct and inverse) cascade processes in MHD shear flows.
Stratified turbulence shows scale- and direction-dependent anisotropy and the coexistence of weak turbulence of internal gravity waves and strong turbulence of eddies. Straightforward application of standard analyses developed in isotropic turbulence sometimes masks important aspects of the anisotropic turbulence. To capture detailed structures of the energy distribution in the wave-number space, it is indispensable to examine the energy distribution with non-integrated spectra by fixing the codimensional wave-number component or in the two-dimensional domain spanned by both the horizontal and vertical wave numbers. Indices which separate the range of the anisotropic weak-wave turbulence in the wave-number space are proposed based on the decomposed energies. In addition, the dominance of the waves in the range is also verified by the small frequency deviation from the linear dispersion relation. In the wave-dominant range, the linear wave periods given by the linear dispersion relation are smaller than approximately one third of the eddy-turnover time. The linear wave periods reflect the anisotropy of the system, while the isotropic Brunt-Vaisala period is used to evaluate the Ozmidov wave number, which is necessarily isotropic. It is found that the time scales in consideration of the anisotropy of the flow field must be appropriately selected to obtain the critical wave number separating the weak-wave turbulence.
In this study, the stability dependence of turbulent Prandtl number ($Pr_t$) is quantified via a novel and simple analytical approach. Based on the variance and flux budget equations, a hybrid length scale formulation is first proposed and its functional relationships to well-known length scales are established. Next, the ratios of these length scales are utilized to derive an explicit relationship between $Pr_t$ and gradient Richardson number. In addition, theoretical predictions are made for several key turbulence variables (e.g., dissipation rates, normalized fluxes). The results from our proposed approach are compared against other competing formulations as well as published datasets. Overall, the agreement between the different approaches is rather good despite their different theoretical foundations and assumptions.