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.

Condensation of water vapor on active cloud condensation nuclei produces micron-size water droplets. To form rain, they must grow rapidly into at least 50-100 $\mu$m droplets. Observations show that this process takes only 15-20 minutes. The unexplained physical mechanism of such fast growth, is crucial for understanding and modeling of rain, and known as "condensation-coalescence bottleneck in rain formation". We show that the recently discovered phenomenon of the tangling clustering instability of small droplets in temperature-stratified turbulence (Phys. Fluids 25, 085104, 2013) results in the formation of droplet clusters with drastically increased droplet number densities. The mechanism of the tangling clustering instability is much more effective than the previously considered by us the inertial clustering instability caused by the centrifugal effect of turbulent vortices. This is the reason of strong enhancement of the collision-coalescence rate inside the clusters. The mean-field theory of the droplet growth developed in this study can be useful for explanation of the observed fast growth of cloud droplets in warm clouds from the initial 1 micron size droplets to 40-50 micron size droplets within 15-20 minutes.

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.

The paper analyzes the phenomenon of turbulent thermal diffusion in the Earth atmosphere, its relation to the turbulent diffusion and its potential impact on aerosol distribution. This phenomenon was predicted theoretically more than 10 years ago and detected recently in the laboratory experiments. This effect causes a non-diffusive flux of aerosols in the direction of the heat flux and results in formation of long-living aerosol layers in the vicinity of temperature inversions. We demonstrated that the theory of turbulent thermal diffusion explains the GOMOS aerosol observations near the tropopause (i.e., the observed shape of aerosol vertical profiles with elevated concentrations located almost symmetrically with respect to temperature profile). In combination with the derived expression for the dependence of the turbulent thermal diffusion ratio on the turbulent diffusion, these measurements yield an independent method for determining the coefficient of turbulent diffusion at the tropopause. We evaluated the impact of turbulent thermal diffusion to the lower-troposphere vertical profiles of aerosol concentration by means of numerical dispersion modelling, and found a regular upward forcing of aerosols with coarse particles affected stronger than fine aerosols.

We advance our prior energy- and flux-budget turbulence closure model (Zilitinkevich et al., 2007, 2008) for the stably stratified atmospheric flows and extend it accounting for additional vertical flux of momentum and additional productions of turbulent kinetic energy, turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). Main effects of IGW are following: the maximal value of the flux Richardson number (universal constant 0.2-0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. In the heterogeneous stratification, when IGW propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. IGW also reduce anisotropy of turbulence and increase the share of TPE in the turbulent total energy. Depending on the direction (downward or upward), IGW either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.

Traditionally, turbulence energetics is characterized by turbulent kinetic energy (TKE) and modelled using solely the TKE budget equation. In stable stratification, TKE is generated by the velocity shear and expended through viscous dissipation and work against buoyancy forces. The effect of stratification is characterized by the ratio of the buoyancy gradient to squared shear, called Richardson number, Ri. It is widely believed that at Ri exceeding a critical value, Ric, local shear cannot maintain turbulence, and the flow becomes laminar. We revise this concept by extending the energy analysis to turbulent potential and total energies (TPE and TTE = TKE + TPE), consider their budget equations, and conclude that TTE is a conservative parameter maintained by shear in any stratification. Hence there is no "energetics Ric", in contrast to the hydrodynamic-instability threshold, Ric-instability, whose typical values vary from 0.25 to 1. We demonstrate that this interval, 0.25<Ri<1, separates two different turbulent regimes: strong mixing and weak mixing rather than the turbulent and the laminar regimes, as the classical concept states. This explains persistent occurrence of turbulence in the free atmosphere and deep ocean at Ri>>1, clarify principal difference between turbulent boundary layers and free flows, and provide basis for improving operational turbulence closure models.

We propose a new turbulence closure model based on the budget equations for the key second moments: turbulent kinetic and potential energies: TKE and TPE (comprising the turbulent total energy: TTE = TKE + TPE) and vertical turbulent fluxes of momentum and buoyancy (proportional to potential temperature). Besides the concept of TTE, we take into account the non-gradient correction to the traditional buoyancy flux formulation. The proposed model grants the existence of turbulence at any gradient Richardson number, Ri. Instead of its critical value separating - as usually assumed - the turbulent and the laminar regimes, it reveals a transition interval, 0.1< Ri <1, which separates two regimes of essentially different nature but both turbulent: strong turbulence at Ri<<1; and weak turbulence, capable of transporting momentum but much less efficient in transporting heat, at Ri>1. Predictions from this model are consistent with available data from atmospheric and lab experiments, direct numerical simulation (DNS) and large-eddy simulation (LES).

A new mean-field theory of turbulent convection is developed. This theory predicts the convective wind instability in a shear-free turbulent convection which causes formation of large-scale semi-organized fluid motions in the form of cells or rolls. Spatial characteristics of these motions, such as the minimum size of the growing perturbations and the size of perturbations with the maximum growth rate, are determined. This study predicts also the existence of the convective shear instability in a sheared turbulent convection which results in generation of convective shear waves with a nonzero hydrodynamic helicity. Increase of shear promotes excitation of the convective shear instability. Applications of the obtained results to the atmospheric turbulent convection and the laboratory experiments on turbulent convection are discussed. This theory can be applied also for the describing a mesogranular turbulent convection in astrophysics.