We develop a theory of three-dimensional slow Rossby waves in rotating spherical density stratified convection. The Rossby waves with frequencies which are much smaller than the rotating frequency, are excited by a non-axisymmetric instability from the equilibrium based on the developed convection. These waves interact with the inertial waves and the density stratified convection. The density stratification is taken into account using the anelastic approximation for very low-Mach-number flows. We study long-term planetary Rossby waves with periods which are larger than two years. We suggest that these waves are related to the Southern Oscillation and El Ni\~{n}o. The El Ni\~{n}o is an irregularly periodical variation in winds and sea surface temperatures over the tropical Pacific Ocean, while the Southern Oscillation is an oscillation in surface air pressure between the tropical eastern and the western Pacific Ocean. The strength of the Southern Oscillation is characterised by the Southern Oscillation Index (SOI). The developed theory is applied for the interpretation of the observed periods of the SOI. This study demonstrates a good agreement between the theoretical predictions and the observations.

Turbulent thermal diffusion is a combined effect of the temperature stratified turbulence and inertia of small particles. It causes the appearance of a non-diffusive turbulent flux of particles in the direction of the turbulent heat flux. This non-diffusive turbulent flux of particles is proportional to the product of the mean particle number density and the effective velocity of inertial particles. The theory of this effect has been previously developed only for small temperature gradients and small Stokes numbers (Phys. Rev. Lett. {\bf 76}, 224, 1996). In this study a generalized theory of turbulent thermal diffusion for arbitrary temperature gradients and Stokes numbers has been developed. The laboratory experiments in the oscillating grid turbulence and in the multi-fan produced turbulence have been performed to validate the theory of turbulent thermal diffusion in strongly stratified turbulent flows. It has been shown that the ratio of the effective velocity of inertial particles to the characteristic vertical turbulent velocity for large Reynolds numbers is less than 1. The effective velocity of inertial particles as well as the effective coefficient of turbulent thermal diffusion increase with Stokes numbers reaching the maximum at small Stokes numbers and decreases for larger Stokes numbers. The effective coefficient of turbulent thermal diffusion also decreases with the mean temperature gradient. It has been demonstrated that the developed theory is in a good agreement with the results of the laboratory experiments.

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.

We study turbulent diffusion of chemically reacting gaseous admixtures in a developed turbulence. In our previous study [Phys. Rev. Lett. {\bf 80}, 69 (1998)] using a path-integral approach for a delta-correlated in time random velocity field, we demonstrated a strong modification of turbulent transport in fluid flows with chemical reactions or phase transitions. In the present study we use the spectral tau approximation, that is valid for large Reynolds and Peclet numbers, and show that turbulent diffusion of the reacting species can be strongly depleted by a large factor that is the ratio of turbulent and chemical times (turbulent Damk\"{o}hler number). We have demonstrated that the derived theoretical dependence of turbulent diffusion coefficient versus the turbulent Damk\"{o}hler number is in a good agreement with that obtained previously in the numerical modelling of a reactive front propagating in a turbulent flow and described by the Kolmogorov-Petrovskii-Piskunov-Fisher equation. We have found that turbulent cross-effects, e.g., turbulent mutual diffusion of gaseous admixtures and turbulent Dufour-effect of the chemically reacting gaseous admixtures, are less sensitive to the values of stoichiometric coefficients. The mechanisms of the turbulent cross-effects are different from the molecular cross effects known in irreversible thermodynamics. In a fully developed turbulence and at large Peclet numbers the turbulent cross-effects are much larger than the molecular ones. The obtained results are applicable also to heterogeneous phase transitions.

We investigated experimentally stably stratified turbulent flows in a lid driven cavity with a non-zero vertical mean temperature gradient in order to identify the parameters governing the mean and turbulent flows and to understand their effects on the momentum and heat transfer. We found that the mean velocity patterns (e.g., the form and the sizes of the large-scale circulations) depend strongly on the degree of the temperature stratification. In the case of strong stable stratification, the strong turbulence region is located in the vicinity of the main large-scale circulation. We detected the large-scale nonlinear oscillations in the case of strong stable stratification which can be interpreted as nonlinear internal gravity waves. The ratio of the main energy-containing frequencies of these waves in velocity and temperature fields in the nonlinear stage is about 2. The amplitude of the waves increases in the region of weak turbulence (near the bottom wall of the cavity), whereby the vertical mean temperature gradient increases.

We study particle clustering in a temperature stratified turbulence with small finite correlation time. It is shown that the temperature stratified turbulence strongly increases the degree of compressibility of particle velocity field. This results in the strong decrease of the threshold for the excitation of the tangling clustering instability even for small particles. The tangling clustering instability in the temperature stratified turbulence is essentially different from the inertial clustering instability that occurs in non-stratified isotropic and homogeneous turbulence. While the inertial clustering instability is caused by the centrifugal effect of the turbulent eddies, the mechanism of the tangling clustering instability is related to the temperature fluctuations generated by the tangling of the mean temperature gradient by the velocity fluctuations. Temperature fluctuations produce pressure fluctuations and cause particle clustering in regions with increased pressure fluctuations. It is shown that the growth rate of the tangling clustering instability is much larger than that of the inertial clustering instability. It is found that depending on the parameters of the turbulence and the mean temperature gradient there is a preferential particle size at which the particle clustering due to the tangling clustering instability is more effective. The particle number density inside the cluster after the saturation of this instability can be in several orders of magnitude larger than the mean particle number density. It is also demonstrated that the evaporation of droplets drastically change the tangling clustering instability, e.g., it increases the instability threshold in the droplet radius. The tangling clustering instability is of a great importance, e.g., in atmospheric turbulence with temperature inversions.

We study experimentally and theoretically transition phenomena caused by the external forcing from Rayleigh-Benard convection with the large-scale circulation (LSC) to the limiting regime of unstably stratified turbulent flow without LSC whereby the temperature field behaves like a passive scalar. In the experiments we use the Rayleigh-B\'enard apparatus with an additional source of turbulence produced by two oscillating grids located nearby the side walls of the chamber. When the frequency of the grid oscillations is larger than 2 Hz, the large-scale circulation (LSC) in turbulent convection is destroyed, and the destruction of the LSC is accompanied by a strong change of the mean temperature distribution. However, in all regimes of the unstably stratified turbulent flow the ratio $\big[(\ell_x \nabla_x T)^2 + (\ell_y \nabla_y T)^2 + (\ell_z \nabla_z T)^2\big] / < \theta^2 >$ varies slightly (even in the range of parameters whereby the behaviour of the temperature field is different from that of the passive scalar). Here $\ell_i$ are the integral scales of turbulence along x, y, z directions, T and \theta are the mean and fluctuating parts of the fluid temperature. At all frequencies of the grid oscillations we have detected the long-term nonlinear oscillations of the mean temperature. The theoretical predictions based on the budget equations for turbulent kinetic energy, turbulent temperature fluctuations and turbulent heat flux, are in agreement with the experimental results.

We have predicted theoretically and detected in laboratory experiments a new type of particle clustering (tangling clustering of inertial particles) in a stably stratified turbulence with imposed mean vertical temperature gradient. In this stratified turbulence a spatial distribution of the mean particle number density is nonuniform due to the phenomenon of turbulent thermal diffusion, that results in formation of a gradient of the mean particle number density, \nabla N, and generation of fluctuations of the particle number density by tangling of the gradient, \nabla N, by velocity fluctuations. The mean temperature gradient, \nabla T, produces the temperature fluctuations by tangling of the gradient, \nabla T, by velocity fluctuations. These fluctuations increase the rate of formation of the particle clusters in small scales. In the laboratory stratified turbulence this tangling clustering is much more effective than a pure inertial clustering that has been observed in isothermal turbulence. In particular, in our experiments in oscillating grid isothermal turbulence in air without imposed mean temperature gradient, the inertial clustering is very weak for solid particles with the diameter 10 microns and Reynolds numbers Re =250. Our theoretical predictions are in a good agreement with the obtained experimental results.

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.

We study an effect of large-scale coherent structures on global properties of turbulent convection in laboratory experiments in air flow in a rectangular chamber with aspect ratios $A \approx 2$ and $A\approx 4$ (with the Rayleigh numbers varying in the range from $5 \times 10^6$ to $10^8$). The large-scale coherent structures comprise the one-cell and two-cell flow patterns. We found that a main contribution to the turbulence kinetic energy production in turbulent convection with large-scale coherent structures is due to the non-uniform large-scale motions. Turbulence in large Rayleigh number convection with coherent structures is produced by shear, rather than by buoyancy. We determined the scalings of global parameters (e.g., the production and dissipation of turbulent kinetic energy, the turbulent velocity and integral turbulent scale, the large-scale shear, etc.) of turbulent convection versus the temperature difference between the bottom and the top walls of the chamber. These scalings are in an agreement with our theoretical predictions. We demonstrated that the degree of inhomogeneity of the turbulent convection with large-scale coherent structures is small.

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 study a large-scale instability in a sheared nonhelical turbulence that causes generation of large-scale vorticity. Three types of the background large-scale flows are considered, i.e., the Couette and Poiseuille flows in a small-scale homogeneous turbulence, and the "log-linear" velocity shear in an inhomogeneous turbulence. It is known that laminar plane Couette flow and antisymmetric mode of laminar plane Poiseuille flow are stable with respect to small perturbations for any Reynolds numbers. We demonstrate that in a small-scale turbulence under certain conditions the large-scale Couette and Poiseuille flows are unstable due to the large-scale instability. This instability causes formation of large-scale vortical structures stretched along the mean sheared velocity. The growth rate of the large-scale instability for the "log-linear" velocity shear is much larger than that for the Couette and Poiseuille background flows. We have found a turbulent analogue of the Tollmien-Schlichting waves in a small-scale sheared turbulence. A mechanism of excitation of turbulent Tollmien-Schlichting waves is associated with a combined effect of the turbulent Reynolds stress-induced generation of perturbations of the mean vorticity and the background sheared motions. These waves can be excited even in a plane Couette flow imposed on a small-scale turbulence when perturbations of mean velocity depend on three spatial coordinates. The energy of these waves is supplied by the small-scale sheared turbulence.

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).

Formation of large-scale coherent structures in a turbulent convection via excitation of large-scale instability is studied. The redistribution of the turbulent heat flux due to non-uniform large-scale motions plays a crucial role in the formation of the coherent large-scale structures in the turbulent convection. The modification of the turbulent heat flux results in strong reduction of the critical Rayleigh number (based on the eddy viscosity and turbulent temperature diffusivity) required for the excitation of the large-scale instability. The mean-field equations which describe the large-scale instability, are solved numerically. We determine the key parameters that affect formation of the large-scale coherent structures in the turbulent convection. In particular, the degree of thermal anisotropy and the lateral background heat flux strongly modify the growth rates of the large-scale instability, the frequencies of the generated convective-shear waves and change the thresholds required for the excitation of the large-scale instability. This study elucidates the origins of the large-scale circulations and rolls in the atmospheric convective boundary layers and the meso-granular structures in the solar convection.

A mechanism of excitation of the large-scale inertial waves in a rotating inhomogeneous turbulence due to an excitation of a large-scale instability is found. This instability is caused by a combined effect of the inhomogeneity of the turbulence and the uniform mean rotation. The source of the large-scale instability is the energy of the small-scale turbulence. We determined the range of parameters at which the large-scale instability occurs, the growth rate of the instability and the frequency of the generated large-scale inertial waves. Geophysical and astrophysical applications of the obtained results are discussed.

We discuss a new phenomenon of turbulent thermal diffusion associated with turbulent transport of aerosols in the atmosphere and in laboratory experiments. The essence of this phenomenon is the appearance of a nondiffusive mean flux of particles in the direction of the mean heat flux, which results in the formation of large-scale inhomogeneities in the spatial distribution of aerosols that accumulate in regions of minimum mean temperature of the surrounding fluid. This effect of turbulent thermal diffusion was detected experimentally. In experiments turbulence was generated by two oscillating grids in two directions of the imposed vertical mean temperature gradient. We used Particle Image Velocimetry to determine the turbulent velocity field, and an Image Processing Technique based on an analysis of the intensity of Mie scattering to determine the spatial distribution of aerosols. Analysis of the intensity of laser light Mie scattering by aerosols showed that aerosols accumulate in the vicinity of the minimum mean temperature due to the effect of turbulent thermal diffusion. Geophysical applications of the obtained results are discussed.

An effect of a mean velocity shear on a turbulence and on the effective force which is determined by the gradient of Reynolds stresses is studied. Generation of a mean vorticity in a homogeneous incompressible turbulent flow with an imposed mean velocity shear due to an excitation of a large-scale instability is found. The instability is caused by a combined effect of the large-scale shear motions (''skew-induced" deflection of equilibrium mean vorticity) and ''Reynolds stress-induced" generation of perturbations of mean vorticity. Spatial characteristics, such as the minimum size of the growing perturbations and the size of perturbations with the maximum growth rate, are determined. This instability and the dynamics of the mean vorticity are associated with the Prandtl's turbulent secondary flows. This instability is similar to the mean-field magnetic dynamo instability. Astrophysical applications of the obtained results are discussed.

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.