The tilt angle, current helicity and twist of solar magnetic fields can be observed in solar active regions. We carried out estimates of these parameters by two ways. Firstly, we consider the model of turbulent convective cells (super-granules) which have a loop floating structure towards the surface of the Sun. Their helical properties are attained during the rising process in the rotating stratified convective zone. The other estimate is obtained from a simple mean-field dynamo model that accounts magnetic helicity conservation. The both values are shown to be capable to give important contributions to the observable tilt, helicity and twist.

A theory of the mean tilt of sunspot bipolar regions (the angle between a line connecting the leading and following sunspots and the solar equator) is developed. A mechanism of formation of the mean tilt is related to the effect of Coriolis force on meso-scale motions of super-granular convection and large-scale meridional circulation. The balance between the Coriolis force and the Lorentz force (the magnetic tension) determines an additional contribution caused by the large-scale magnetic field to the mean tilt of the sunspot bipolar regions at low latitudes. The latitudinal dependence of the solar differential rotation affects the mean tilt which can explain deviations from the Joy's law for the sunspot bipolar regions at high latitudes. The obtained theoretical results and performed numerical simulations based on the nonlinear mean-field dynamo theory which takes into account conservation of the total magnetic helicity and the budget equation for the evolution of the Wolf number density, are in agreement with observational data of the mean tilt of sunspot bipolar regions over individual solar cycles 15 - 24.

We find an instability resulting in generation of large-scale vorticity in a fast rotating small-scale turbulence or turbulent convection with inhomogeneous fluid density along the rotational axis in anelastic approximation. The large-scale instability causes excitation of two modes: (i) the mode with dominant vertical vorticity and with the mean velocity being independent of the vertical coordinate; (ii) the mode with dominant horizontal vorticity and with the mean momentum being independent of the vertical coordinate. The mode with the dominant vertical vorticity can be excited in a fast rotating density stratified hydrodynamic turbulence or turbulent convection. For this mode, the mean entropy is depleted inside the cyclonic vortices, while it is enhanced inside the anti-cyclonic vortices. The mode with the dominant horizontal vorticity can be excited only in a fast rotating density stratified turbulent convection. The developed theory may be relevant for explanation of an origin of large spots observed as immense storms in great planets, e.g., the Great Red Spot in Jupiter and large spots in Saturn. It may be also useful for explanation of an origin of high-latitude spots in rapidly rotating late-type stars.

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.

We present a physical mechanism of formation of tilt angles of sunspots due to the process of formation of active regions below the solar photosphere. The contribution of Coriolis force factors on large-scale flows of super-granular convection in turbulent media has been investigated in details. On the basis of earlier works by Kleeorin et al. (2016) and Safiullin et al. (2018) we give physical estimates of orders of magnitude of the effect and estimate the tilt angles near the solar equator, in the "Royal" zone of solar activity. The above model is based on the balance of the small-scale and large scale magnetic helicities and describes in details the sunspot formation process over the last five solar cycles (since 1964). We adopt this model for a wider class of manifestations of solar activity. We present latitudinal dependence of the mean tilt on these five solar cycles and time-latitudinal diagrams over a limited range of latitudes and phases of the solar cycle.

We develop a mean-field theory of compressibility effects in turbulent magnetohydrodynamics and passive scalar transport using the quasi-linear approximation and the spectral $\tau$-approach. We find that compressibility decreases the $\alpha$ effect and the turbulent magnetic diffusivity both at small and large magnetic Reynolds numbers, Rm. Similarly, compressibility decreases the turbulent diffusivity for passive scalars both at small and large P\'eclet numbers, Pe. On the other hand, compressibility does not affect the effective pumping velocity of the magnetic field for large Rm, but it decreases it for small Rm. Density stratification causes turbulent pumping of passive scalars, but it is found to become weaker with increasing compressibility. No such pumping effect exists for magnetic fields. However, compressibility results in a new passive scalar pumping effect from regions of low to high turbulent intensity both for small and large P\'eclet numbers. It can be interpreted as compressible turbophoresis of noninertial particles and gaseous admixtures, while the classical turbophoresis effect exists only for inertial particles and causes them to be pumped to regions with lower turbulent intensity.

An asymmetry in the number density of left- and right-handed fermions is known to give rise to a new term in the induction equation that can result in a dynamo instability. At high temperatures, when a chiral asymmetry can survive for long enough, this chiral dynamo instability can amplify magnetic fields efficiently, which in turn drive turbulence via the Lorentz force. While it has been demonstrated in numerical simulations that this chiral magnetically driven turbulence exists and strongly affects the dynamics of the magnetic field, the details of this process remain unclear. The goal of this paper is to analyse the energetics of chiral magnetically driven turbulence and its effect on the generation and dynamics of the magnetic field using direct numerical simulations. We study these effects for different initial conditions, including a variation of the initial chiral chemical potential and the magnetic Prandtl number, Pm. In particular, we determine the ratio of kinetic to magnetic energy, Upsilon, in chiral magnetically driven turbulence. Within the parameter space explored in this study, Upsilon reaches a value of approximately 0.064-0.074 -independently of the initial chiral asymmetry and for Pm=1. Our simulations suggest, that Upsilon decreases as a power law when increasing Pm by decreasing the viscosity. While the exact scaling depends on the details of the fitting criteria and the Reynolds number regime, an approximate result of Upsilon(Pm)=0.1 Pm^{-0.4} is reported. Using the findings from our numerical simulations, we analyse the energetics of chiral magnetically driven turbulence in the early Universe.

The formation of sunspots and starspots is not yet fully understood and is therefore one of the major open problems in solar and stellar physics. Magnetic flux concentrations can be produced by the negative effective magnetic pressure instability (NEMPI). This instability is strongly suppressed by rotation. However, the presence of an outer coronal envelope was previously found to strengthen the flux concentrations and make them more prominent. It also allows for the formation of bipolar regions (BRs). We want to know whether the presence of an outer coronal envelope also changes the excitation conditions and the rotational dependence of NEMPI. We use direct numerical simulations and mean-field simulations. We adopt a simple two-layer model of turbulence that mimics the jump between the convective turbulent and coronal layers below and above the surface of a star, respectively. The computational domain is Cartesian and located at a certain latitude of a rotating sphere. We investigate the effects of rotation on NEMPI by changing the Coriolis number, the latitude, and the box resolution. Rotation has a strong impact on the process of BR formation. Even rather slow rotation is found to suppress their formation. However, increasing the imposed magnetic field strength also makes the structures stronger and alleviates the rotational suppression somewhat. The presence of a coronal layer itself does not significantly alleviate the effects of rotational suppression.

We apply a nonlinear mean-field dynamo model which includes a budget equation for the dynamics of Wolf numbers to predict solar activity. This dynamo model takes into account the algebraic and dynamic nonlinearities of the alpha effect, where the equation for the dynamic nonlinearity is derived from the conservation law for the magnetic helicity. The budget equation for the evolution of the Wolf number is based on a formation mechanism of sunspots related to the negative effective magnetic pressure instability. This instability redistributes the magnetic flux produced by the mean-field dynamo. To predict solar activity on the time scale of one month we use a method based on a combination of the numerical solution of the nonlinear mean-field dynamo equations and the artificial neural network. A comparison of the results of the prediction of the solar activity with the observed Wolf numbers demonstrates a good agreement between the forecast and observations.

Using direct numerical simulations (DNS), we study laminar and turbulent dynamos in chiral magnetohydrodynamics (MHD) with an extended set of equations that accounts for an additional contribution to the electric current due to the chiral magnetic effect (CME). This quantum phenomenon originates from an asymmetry between left- and right-handed relativistic fermions in the presence of a magnetic field and gives rise to a chiral dynamo. We show that the magnetic field evolution proceeds in three stages: (1) a small-scale chiral dynamo instability; (2) production of chiral magnetically driven turbulence and excitation of a large-scale dynamo instability due to a new chiral effect (alpha_mu effect); and (3) saturation of magnetic helicity and magnetic field growth controlled by a conservation law for the total chirality. The $\alpha_\mu$ effect becomes dominant at large fluid and magnetic Reynolds numbers and is not related to kinetic helicity. The growth rate of the large-scale magnetic field and its characteristic scale measured in the numerical simulations agree well with theoretical predictions based on mean-field theory. The previously discussed two-stage chiral magnetic scenario did not include stage (2) during which the characteristic scale of magnetic field variations can increase by many orders of magnitude. Based on the findings from numerical simulations, the relevance of the CME and the chiral effects revealed in the relativistic plasma of the early universe and of proto-neutron stars are discussed.

We discuss a mean-field theory of generation of large-scale vorticity in a rotating density stratified developed turbulence with inhomogeneous kinetic helicity. We show that the large-scale nonuniform flow is produced due to ether a combined action of a density stratified rotating turbulence and uniform kinetic helicity or a combined effect of a rotating incompressible turbulence and inhomogeneous kinetic helicity. These effects result in the formation of a large-scale shear, and in turn its interaction with the small-scale turbulence causes an excitation of the large-scale instability (known as a vorticity dynamo) due to a combined effect of the large-scale shear and Reynolds stress-induced generation of the mean vorticity. The latter is due to the effect of large-scale shear on the Reynolds stress. A fast rotation suppresses this large-scale instability.

A mean-field theory of differential rotation in a density stratified turbulent convection has been developed. This theory is based on a combined effect of the turbulent heat flux and anisotropy of turbulent convection on the Reynolds stress. A coupled system of dynamical budget equations consisting in the equations for the Reynolds stress, the entropy fluctuations and the turbulent heat flux has been solved. To close the system of these equations, the spectral tau approach which is valid for large Reynolds and Peclet numbers, has been applied. The adopted model of the background turbulent convection takes into account an increase of the turbulence anisotropy and a decrease of the turbulent correlation time with the rotation rate. This theory yields the radial profile of the differential rotation which is in agreement with that for the solar differential rotation.

It is known that unconfined dust explosions consist of a relatively weak primary (turbulent) deflagrations followed by a devastating secondary explosion. The secondary explosion may propagate with a speed of up to 1000 m/s producing overpressures of over 8-10 atm. Since detonation is the only established theory that allows a rapid burning producing a high pressure that can be sustained in open areas, the generally accepted view was that the mechanism explaining the high rate of combustion in dust explosions is deflagration to detonation transition. In the present work we propose a theoretical substantiation of the alternative propagation mechanism explaining origin of the secondary explosion producing the high speeds of combustion and high overpressures in unconfined dust explosions. We show that clustering of dust particles in a turbulent flow gives rise to a significant increase of the thermal radiation absorption length ahead of the advancing flame front. This effect ensures that clusters of dust particles are exposed to and heated by the radiation from hot combustion products of large gaseous explosions sufficiently long time to become multi-point ignition kernels in a large volume ahead of the advancing flame front. The ignition times of fuel-air mixture by the radiatively heated clusters of particles is considerably reduced compared to the ignition time by the isolated particle. The radiation-induced multi-point ignitions of a large volume of fuel-air ahead of the primary flame efficiently increase the total flame area, giving rise to the secondary explosion, which results in high rates of combustion and overpressures required to account for the observed level of overpressures and damages in unconfined dust explosions, such as e.g. the 2005 Buncefield explosion and several vapor cloud explosions of severity similar to that of the Buncefield incident.

The magnetohydrodynamic (MHD) description of plasmas with relativistic particles necessarily includes an additional new field, the chiral chemical potential associated with the axial charge (i.e., the number difference between right- and left-handed relativistic fermions). This chiral chemical potential gives rise to a contribution to the electric current density of the plasma (\emph{chiral magnetic effect}). We present a self-consistent treatment of the \emph{chiral MHD equations}, which include the back-reaction of the magnetic field on a chiral chemical potential and its interaction with the plasma velocity field. A number of novel phenomena are exhibited. First, we show that the chiral magnetic effect decreases the frequency of the Alfv\'{e}n wave for incompressible flows, increases the frequencies of the Alfv\'{e}n wave and of the fast magnetosonic wave for compressible flows, and decreases the frequency of the slow magnetosonic wave. Second, we show that, in addition to the well-known laminar chiral dynamo effect, which is not related to fluid motions, there is a dynamo caused by the joint action of velocity shear and chiral magnetic effect. In the presence of turbulence with vanishing mean kinetic helicity, the derived mean-field chiral MHD equations describe turbulent large-scale dynamos caused by the chiral alpha effect, which is dominant for large fluid and magnetic Reynolds numbers. The chiral alpha effect is due to an interaction of the chiral magnetic effect and fluctuations of the small-scale current produced by tangling magnetic fluctuations (which are generated by tangling of the large-scale magnetic field by sheared velocity fluctuations). These dynamo effects may have interesting consequences in the dynamics of the early universe, neutron stars, and the quark--gluon plasma.

The presence of asymmetry between fermions of opposite handedness in plasmas of relativistic particles can lead to exponential growth of a helical magnetic field via a small-scale chiral dynamo instability known as the chiral magnetic effect. Here, we show, using dimensional arguments and numerical simulations, that this process produces through the Lorentz force chiral magnetically driven turbulence. A k^{-2} magnetic energy spectrum emerges via inverse transfer over a certain range of wavenumbers k. The total chirality (magnetic helicity plus normalized chiral chemical potential) is conserved in this system. Therefore, as the helical magnetic field grows, most of the total chirality gets transferred into magnetic helicity until the chiral magnetic effect terminates. Quantitative results for height, slope, and extent of the spectrum are obtained. Consequences of this effect for cosmic magnetic fields are discussed.

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.

It is known that unconfined dust explosions typically starts off with a relatively weak primary flame followed by a severe secondary explosion. We show that clustering of dust particles in a temperature stratified turbulent flow ahead of the primary flame may give rise to a significant increase in the radiation penetration length. These particle clusters, even far ahead of the flame, are sufficiently exposed and heated by the radiation from the flame to become ignition kernels capable to ignite a large volume of fuel-air mixtures. This efficiently increases the total flame surface area and the effective combustion speed, defined as the rate of reactant consumption of a given volume. We show that this mechanism explains the high rate of combustion and overpressures required to account for the observed level of damage in unconfined dust explosions, e.g., at the 2005 Buncefield vapor-cloud explosion. The effect of the strong increase of radiation transparency due to turbulent clustering of particles goes beyond the state-of-the-art of the application to dust explosions and has many implications in atmospheric physics and astrophysics.

Sunspots are of basic interest in the study of the Sun. Their relevance ranges from them being an activity indicator of magnetic fields to being the place where coronal mass ejections and flares erupt. They are therefore also an important ingredient of space weather. Their formation, however, is still an unresolved problem in solar physics. Observations utilize just 2D surface information near the spot, but it is debatable how to infer deep structures and properties from local helioseismology. For a long time, it was believed that flux tubes rising from the bottom of the convection zone are the origin of the bipolar sunspot structure seen on the solar surface. However, this theory has been challenged, in particular recently by new surface observation, helioseismic inversions, and numerical models of convective dynamos. In this article we discuss another theoretical approach to the formation of sunspots: the negative effective magnetic pressure instability. This is a large-scale instability, in which the total (kinetic plus magnetic) turbulent pressure can be suppressed in the presence of a weak large-scale magnetic field, leading to a converging downflow, which eventually concentrates the magnetic field within it. Numerical simulations of forced stratified turbulence have been able to produce strong super-equipartition flux concentrations, similar to sunspots at the solar surface. In this framework, sunspots would only form close to the surface due to the instability constraints on stratification and rotation. Additionally, we present some ideas from local helioseismology, where we plan to use the Hankel analysis to study the pre-emergence phase of a sunspot and to constrain its deep structure and formation mechanism.

Recent direct numerical simulations (DNS) of large-scale turbulent dynamos in strongly stratified layers have resulted in surprisingly sharp bipolar structures at the surface. Here we present new DNS of helically and non-helically forced turbulence with and without rotation and compare with corresponding mean-field simulations (MFS) to show that these structures are a generic outcome of a broader class of dynamos in density-stratified layers. The MFS agree qualitatively with the DNS, but the period of oscillations tends to be longer in the DNS. In both DNS and MFS, the sharp structures are produced by converging flows at the surface and might be driven in nonlinear stage of evolution by the Lorentz force associated with the large-scale dynamo-driven magnetic field if the dynamo number is at least 2.5 times supercritical.

In the presence of strong density stratification, hydromagnetic turbulence attains qualitatively new properties: the formation of magnetic flux concentrations. We review here the theoretical foundations of this mechanism in terms of what is now called the negative effective magnetic pressure instability. We also present direct numerical simulations of forced turbulence in strongly stratified layers and discuss the qualitative and quantitative similarities with corresponding mean-field simulations. Finally, the relevance to sunspot formation is discussed.

We investigate the dynamics of solar activity using a nonlinear one-dimensional dynamo model and a phenomenological equation for the evolution of Wolf numbers. This system of equations is solved numerically. We take into account the algebraic and dynamic nonlinearities of the alpha effect. The dynamic nonlinearity is related to the evolution of a small-scale magnetic helicity, and it leads to a complicated behavior of solar activity. The evolution equation for the Wolf number is based on a mechanism of formation of magnetic spots as a result of the negative effective magnetic pressure instability (NEMPI). This phenomenon was predicted 25 years ago and has been investigated intensively in recent years through direct numerical simulations and mean-field simulations. The evolution equation for the Wolf number includes the production and decay of sunspots. Comparison between the results of numerical simulations and observational data of Wolf numbers shows a 70 % correlation over all intervals of observation (about 270 years). We determine the dependence of the maximum value of the Wolf number versus the period of the cycle and the asymmetry of the solar cycles versus the amplitude of the cycle. These dependencies are in good agreement with observations.

This work presents an extensive study of the previously discovered formation of bipolar flux concentrations in a two-layer model. We interpret the formation process in terms of negative effective magnetic pressure instability (NEMPI), which is a possible mechanism to explain the origin of sunspots. In our simulations, we use a Cartesian domain of isothermal stratified gas that is divided into two layers. In the lower layer, turbulence is forced with transverse nonhelical random waves, whereas in the upper layer no flow is induced. A weak uniform magnetic field is imposed in the entire domain at all times. In this study we vary the stratification by changing the gravitational acceleration, magnetic Reynolds number, strength of the imposed magnetic field, and size of the domain to investigate their influence on the formation process. Bipolar magnetic structure formation takes place over a large range of parameters. The magnetic structures become more intense for higher stratification until the density contrast becomes around $100$ across the turbulent layer. For the Reynolds numbers considered, magnetic flux concentrations are generated at magnetic Prandtl number between 0.1 and 1. The magnetic field in bipolar regions increases with higher imposed field strength until the field becomes comparable to the equipartition field strength of the turbulence. A larger horizontal extent enables the flux concentrations to become stronger and more coherent. The size of the bipolar structures turns out to be independent of the domain size. In the case of bipolar region formation, we find an exponential growth of the large-scale magnetic field, which is indicative of a hydromagnetic instability. Additionally, the flux concentrations are correlated with strong large-scale downward and converging flows. These findings imply that NEMPI is responsible for magnetic flux concentrations.

We consider strongly stratified forced turbulence in a plane-parallel layer with helicity and corresponding large-scale dynamo action in the lower part and non-helical turbulence in the upper. The magnetic field is found to develop strongly concentrated bipolar structures near the surface. They form elongated bands with a sharp interface between opposite polarities. Unlike earlier experiments with imposed magnetic field, the inclusion of rotation does not strongly suppress the formation of these structures. We perform a systematic numerical study of this phenomenon by varying magnetic Reynolds number, scale separation ratio, and Coriolis number. We focus on the formation of a current sheet between bipolar regions where reconnection of oppositely oriented field lines occurs. We determine the reconnection rate by measuring either the inflow velocity in the vicinity of the current sheet or by measuring the electric field in the reconnection region. We demonstrate that for large Lundquist numbers, S>10^3, the reconnection rate is nearly independent of S in agreement with results of recent numerical simulations performed by other groups in simpler settings.

(abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation.

We derive equations for the mean entropy and the mean internal energy in the low-Mach-number temperature stratified turbulence (i.e., for turbulent convection or stably stratified turbulence), and show that turbulent flux of entropy is given by ${\bf F}_s=\overline{\rho} \, \overline{{\bf u} s}$, where $\overline{\rho}$ is the mean fluid density, $s$ are fluctuations of entropy and overbars denote averaging over an ensemble of turbulent velocity field, ${\bf u}$. We demonstrate that the turbulent flux of entropy is different from the turbulent convective flux, ${\bf F}_c=\overline{T} \, \overline{\rho} \, \overline{{\bf u} s}$, of the fluid internal energy, where $\overline{T}$ is the mean fluid temperature. This turbulent convective flux is well-known in the astrophysical and geophysical literature, and it cannot be used as a turbulent flux in the equation for the mean entropy. This result is exact for low-Mach-number temperature stratified turbulence and is independent of the model used. We also derive equations for the velocity-entropy correlation, $\overline{{\bf u} s}$, in the limits of small and large Peclet numbers, using the quasi-linear approach and the spectral tau approximation, respectively. This study is important in view of different applications to the astrophysical and geophysical temperature stratified 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.

Recent work by Mitra et al. (2014) has shown that in strongly stratified forced two-layer turbulence with helicity and corresponding large-scale dynamo action in the lower layer, a magnetic field occurs in the upper layer in the form of sharply bounded bipolar magnetic spots. Here we extend this model to spherical wedge geometry covering the northern hemisphere up to 75{\deg} latitude and an azimuthal extent of 180{\deg}. The kinetic helicity and therefore also the large-scale magnetic field are strongest at low latitudes. For moderately strong stratification, several bipolar spots form that fill eventually the full longitudinal extent. At early times, the polarity of spots reflects the orientation of the underlying azimuthal field, as expected from {\Omega}-shaped flux loops. At late times their tilt changes such that there is a radial field of opposite orientation at different latitudes separated by about 10{\deg}. Our model demonstrates the spontaneous formation of spots of sizes much larger than the pressure scale height. Their tendency to produce filling factors close to unity is argued to be reminiscent of highly active stars. We confirm that strong stratification and strong scale separation are essential ingredients behind magnetic spot formation, which appears to be associated with downflows at larger depths.

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 mean-field theory of magnetized stellar convection gives rise to the two possibility of distinct instabilities: the large-scale dynamo instability, operating in the bulk of the convection zone, and a negative effective magnetic pressure instability (NEMPI) operating in the strongly stratified surface layers. The latter might be important in connection with magnetic spot formation, but the growth rate of NEMPI is suppressed with increasing rotation rates, although recent direct numerical simulations (DNS) have shown a subsequent increase in the growth rate. We examine quantitatively whether this increase in the growth rate of NEMPI can be explained by an alpha squared mean-field dynamo, and whether both NEMPI and the dynamo instability can operate at the same time. We use both DNS and mean-field simulations (MFS) to solve the underlying equations numerically either with or without an imposed horizontal field. We use the test-field method to compute relevant dynamo coefficients. DNS show that magnetic flux concentrations are still possible up to rotation rates above which the large-scale dynamo effect produces mean magnetic fields. The resulting DNS growth rates are quantitatively well reproduced with MFS. As expected, for weak or vanishing rotation, the growth rate of NEMPI increases with increasing gravity, but there is a correction term for strong gravity and large turbulent magnetic diffusivity. Magnetic flux concentrations are still possible for rotation rates above which dynamo action takes over. For the solar rotation rate, the corresponding turbulent turnover time is about 5 hours, with dynamo action commencing in the layers beneath.

We perform direct numerical simulations of the equations of magnetohydrodynamics with external random forcing and in the presence of gravity. The domain is divided into two parts: a lower layer where the forcing is helical and an upper layer where the helicity of the forcing is zero with a smooth transition in between. At early times, a large-scale helical dynamo develops in the bottom layer. At later times the dynamo saturates, but the vertical magnetic field continues to develop and rises to form dynamic bipolar structures at the top, which later disappear and reappear. Some of the structures look similar to $\delta$ spots observed in the Sun. This is the first example of magnetic flux concentrations, owing to strong density stratification, from self-consistent dynamo simulations that generate bipolar, super-equipartition strength, magnetic structures whose energy density can exceeds the turbulent kinetic energy by even a factor of ten.

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.

In the presence of strong density stratification, turbulence can lead to a large-scale instability of a horizontal magnetic field if its strength is in a suitable range (within a few percent of the turbulent equipartition value). This instability is related to a suppression of the turbulent pressure so that the turbulence contribution to the mean magnetic pressure becomes negative. This results in the excitation of a negative effective magnetic pressure instability (NEMPI). This instability has so far only been studied for an imposed magnetic field. We want to know how NEMPI works when the mean magnetic field is generated self-consistently by an $\alpha^2$ dynamo, whether it is affected by global spherical geometry, and whether it can influence the properties of the dynamo itself. We adopt the mean-field approach which has previously been shown to provide a realistic description of NEMPI in direct numerical simulations. We assume axisymmetry and solve the mean-field equations with the Pencil-Code for an adiabatic stratification at a total density contrast in the radial direction of approximately 4 orders of magnitude. NEMPI is found to work when the dynamo-generated field is about 4% of the equipartition value, which is achieved through strong $\alpha$ quenching. This instability is excited in the top 5% of the outer radius provided the density contrast across this top layer is at least 10. NEMPI is found to occur at lower latitudes when the mean magnetic field is stronger. For weaker fields, NEMPI can make the dynamo oscillatory with poleward migration. NEMPI is a viable mechanism for producing magnetic flux concentrations in a strongly stratified spherical shell in which a magnetic field is generated by a strongly quenched $\alpha$ effect dynamo.

Strongly stratified hydromagnetic turbulence has recently been identified as a candidate for explaining the spontaneous formation of magnetic flux concentrations by the negative effective magnetic pressure instability (NEMPI). Much of this work has been done for isothermal layers for which the density scale height is constant throughout. We now study the validity of earlier conclusions about the size and growth rate of magnetic structures in the case of polytropic layers, which scale height decreases sharply towards the surface. To allow for a continuous transition from isothermal to polytropic layers, we employ a generalization of the exponential function known as the q-exponential. Now, the top of the polytropic layer shifts with the polytropic index such that the scale height at some reference height is always the same. We use both mean-field and direct numerical simulations of forced stratified turbulence to determine the resulting flux concentrations. Magnetic structures begin to form at a depth where magnetic field strength is about 3-4% the local equipartition field strength with respect to the turbulent kinetic energy. Unlike the isothermal case where stronger fields can give rise to magnetic flux concentrations at larger depths, in the polytropic one the growth rates decreases for structures deeper down. For vertical fields, magnetic structures of super-equipartition strengths are formed because such fields survive downward advection, unlike NEMPI under horizontal magnetic fields. The horizontal cross-section of such structures is approximately circular. Results based on isothermal models can be applied locally to polytropic layers. For vertical fields, magnetic flux concentrations of super-equipartition strengths form, which supports suggestions that sunspot formation might be a shallow phenomenon.

Strongly stratified hydromagnetic turbulence has previously been found to produce magnetic flux concentrations if the domain is large enough compared with the size of turbulent eddies. Mean-field simulations (MFS) using parameterizations of the Reynolds and Maxwell stresses show a negative effective magnetic pressure instability and have been able to reproduce many aspects of direct numerical simulations (DNS) regarding the growth rate of this large-scale instability, shape of the resulting magnetic structures, and their height as a function of magnetic field strength. Unlike the case of an imposed horizontal field, for a vertical one, magnetic flux concentrations of equipartition strength with the turbulence can be reached. This results in magnetic spots that are reminiscent of sunspots. Here we want to find out under what conditions magnetic flux concentrations with vertical field occur and what their internal structure is. We use a combination of MFS, DNS, and implicit large-eddy simulations to characterize the resulting magnetic flux concentrations in forced isothermal turbulence with an imposed vertical magnetic field. We confirm earlier results that in the kinematic stage of the large-scale instability the horizontal wavelength of structures is about 10 times the density scale height. At later times, even larger structures are being produced in a fashion similar to inverse spectral transfer in helically driven turbulence. Using turbulence simulations, we find that magnetic flux concentrations occur for different values of the Mach number between 0.1 and 0.7. DNS and MFS show magnetic flux tubes with mean-field energies comparable to the turbulent kinetic energy. The resulting vertical magnetic flux tubes are being confined by downflows along the tubes and corresponding inflow from the sides, which keep the field concentrated.

We report the spontaneous formation of bipolar magnetic structures in direct numerical simulations of stratified forced turbulence with an outer coronal envelope. The turbulence is forced with transverse random waves only in the lower (turbulent) part of the domain. Our initial magnetic field is either uniform in the entire domain or confined to the turbulent layer. After about 1--2 turbulent diffusion times, a bipolar magnetic region of vertical field develops with two coherent circular structures that live during one turbulent diffusion time, and then decay during 0.5 turbulent diffusion times. The resulting magnetic field strengths inside the bipolar region are comparable to the equipartition value with respect to the turbulent kinetic energy. The bipolar magnetic region forms a loop-like structure in the upper coronal layer. We associate the magnetic structure formation with the negative effective magnetic pressure instability in the two-layer model.

Recent studies have demonstrated that in fully developed turbulence, the effective magnetic pressure of a large-scale field (non-turbulent plus turbulent contributions) can become negative. In the presence of strongly stratified turbulence, this was shown to lead to a large-scale instability that produces spontaneous magnetic flux concentrations. Furthermore, using a horizontal magnetic field, elongated flux concentrations with a strength of a few per cent of the equipartition value were found. Here we show that a uniform vertical magnetic field leads to circular magnetic spots of equipartition field strengths. This could represent a minimalistic model of sunspot formation and highlights the importance of two critical ingredients: turbulence and strong stratification. Radiation, ionization, and supergranulation may be important for realistic simulations, but are not critical at the level of a minimalistic model of magnetic spot formation.

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.

In a strongly stratified turbulent layer, a uniform horizontal magnetic field can become unstable to spontaneously form local flux concentrations due to a negative contribution of turbulence to the large-scale (mean-field) magnetic pressure. This mechanism, called the negative effective magnetic pressure instability (NEMPI), is of interest in connection with dynamo scenarios where most of the magnetic field resides in the bulk of the convection zone, and not at the bottom. Recent work using the mean-field hydromagnetic equations has shown that NEMPI becomes suppressed at rather low rotation rates with Coriolis numbers as low as 0.1.}{Here we extend these earlier investigations by studying the effects of rotation both on the development of NEMPI and on the effective magnetic pressure. We also quantify the kinetic helicity from direct numerical simulations (DNS) and compare with earlier work.}{To calculate the rotational effect on the effective magnetic pressure we consider both DNS and analytical studies using the $\tau$ approach. To study the effects of rotation on the development of NEMPI we use both DNS and mean-field calculations of the 3D hydromagnetic equations in a Cartesian domain.}{We find that the growth rates of NEMPI from earlier mean-field calculations are well reproduced with DNS, provided the Coriolis number is below about 0.06. In that case, kinetic and magnetic helicities are found to be weak. For faster rotation, dynamo action becomes possible. However, there is an intermediate range of rotation rates where dynamo action on its own is not yet possible, but the rotational suppression of NEMPI is being alleviated.}{Production of magnetic flux concentrations through the suppression of turbulent pressure appears to be possible only in the upper-most layers of the Sun, where the convective turnover time is less than 2 hours.}

We investigate to what extent the current helicity distribution observed in solar active regions is compatible with solar dynamo models. We use an advanced 2D mean-field dynamo model with dynamo action largely concentrated near the bottom of the convective zone, and dynamo saturation based on the evolution of the magnetic helicity and algebraic quenching. For comparison, we also studied a more basic 2D mean-field dynamo model with simple algebraic alpha quenching only. Using these numerical models we obtain butterfly diagrams for both the small-scale current helicity and the large-scale magnetic helicity, and compare them with the butterfly diagram for the current helicity in active regions obtained from observations. This comparison shows that the current helicity of active regions, as estimated by $-A \cdot B$ evaluated at the depth from which the active region arises, resembles the observational data much better than the small-scale current helicity calculated directly from the helicity evolution equation. Here $B$ and $A$ are respectively the dynamo generated mean magnetic field and its vector potential.

The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.

Using simulations of slowly rotating stratified turbulence, we show that the alpha effect responsible for the generation of astrophysical magnetic fields is proportional to the logarithmic gradient of kinetic energy density rather than that of momentum, as was previously thought. This result is in agreement with a new analytic theory developed in this paper for large Reynolds numbers. Thus, the contribution of density stratification is less important than that of turbulent velocity. The alpha effect and other turbulent transport coefficients are determined by means of the test-field method. In addition to forced turbulence, we also investigate supernova-driven turbulence and stellar convection. In some cases (intermediate rotation rate for forced turbulence, convection with intermediate temperature stratification, and supernova-driven turbulence) we find that the contribution of density stratification might be even less important than suggested by the analytic theory.

We present preliminary results from high resolution magneto-convection simulations where we find the formation of flux concentrations from an initially uniform magnetic field. We compute the effective magnetic pressure but find that the concentrations appear also in places where it is positive. The structures appear in roughly ten convective turnover times and live close to a turbulent diffusion time. The time scales are compatible with the negative effective magnetic pressure instability (NEMPI), although structure formation is not restricted to regions where the effective magnetic pressure is negative.

Using direct numerical simulations (DNS) and mean-field simulations (MFS), the effects of non-uniformity of the magnetic field on the suppression of the turbulent pressure is investigated. This suppression of turbulent pressure can lead to an instability which, in turn, makes the mean magnetic field even more non-uniform. This large-scale instability is caused by a resulting negative contribution of small-scale turbulence to the effective (mean-field) magnetic pressure. We show that enhanced mean current density increases the suppression of the turbulent pressure. The instability leads to magnetic flux concentrations in which the normalized effective mean-field pressure is reduced to a certain value at all depths within a structure.

We investigate the effects of weakly and strongly stratified turbulent convection on the mean effective Lorentz force, and especially on the mean effective magnetic pressure. Earlier studies with isotropically forced non-stratified and stratified turbulence have shown that the contribution of the turbulence to the mean magnetic pressure is negative for mean horizontal magnetic fields that are smaller than the equipartition strength, so that the effective mean magnetic pressure that takes into account the turbulence effects, can be negative. Compared with earlier cases of forced turbulence with an isothermal equation of state, we find that the turbulence effect is similar to or even stronger in the present case of turbulent convection. This is argued to be due to the anisotropy of turbulence in the vertical direction. Another important difference compared with earlier studies is the presence of an evolution equation for the specific entropy. Mean-field modelling with entropy evolution indicates that the negative effective magnetic pressure can still lead to a large-scale instability which forms local flux concentrations, even though the specific entropy evolution tends to have a stabilizing effect when applied to a stably stratified (e.g., isothermal) layer. It is argued that this large-scale instability could be important for the formation of solar large-scale magnetic structures such as active regions.

In this study we discuss two key issues related to a small-scale dynamo instability at low magnetic Prandtl numbers and large magnetic Reynolds numbers, namely: (i) the scaling for the growth rate of small-scale dynamo instability in the vicinity of the dynamo threshold; (ii) the existence of the Golitsyn spectrum of magnetic fluctuations in small-scale dynamos. There are two different asymptotics for the small-scale dynamo growth rate: in the vicinity of the threshold of the excitation of the small-scale dynamo instability, $\lambda \propto \ln({\rm Rm}/ {\rm Rm}^{\rm cr})$, and when the magnetic Reynolds number is much larger than the threshold of the excitation of the small-scale dynamo instability, $\lambda \propto {\rm Rm}^{1/2}$, where ${\rm Rm}^{\rm cr}$ is the small-scale dynamo instability threshold in the magnetic Reynolds number ${\rm Rm}$. We demonstrated that the existence of the Golitsyn spectrum of magnetic fluctuations requires a finite correlation time of the random velocity field. On the other hand, the influence of the Golitsyn spectrum on the small-scale dynamo instability is minor. This is the reason why it is so difficult to observe this spectrum in direct numerical simulations for the small-scale dynamo with low magnetic Prandtl numbers.

We demonstrate that the current helicity observed in solar active regions traces the magnetic helicity of the large-scale dynamo generated field. We use an advanced 2D mean-field dynamo model with dynamo saturation based on the evolution of the magnetic helicity and algebraic quenching. For comparison, we also studied a more basic 2D mean-field dynamo model with simple algebraic alpha quenching only. Using these numerical models we obtained butterfly diagrams both for the small-scale current helicity and also for the large-scale magnetic helicity, and compared them with the butterfly diagram for the current helicity in active regions obtained from observations. This comparison shows that the current helicity of active regions, as estimated by $-{\bf A \cdot B}$ evaluated at the depth from which the active region arises, resembles the observational data much better than the small-scale current helicity calculated directly from the helicity evolution equation. Here ${\bf B}$ and ${\bf A}$ are respectively the dynamo generated mean magnetic field and its vector potential. A theoretical interpretation of these results is given.

We show that an alpha effect is driven by the cosmic ray Bell instability exciting left-right asymmetric turbulence. Alfven waves of a preferred polarization have maximally helical motion, because the transverse motion of each mode is parallel to its curl. We show how large-scale Alfven modes, when rendered unstable by cosmic ray streaming, can create new net flux over any finite region, in the direction of the original large-scale field. We perform direct numerical simulations (DNS) of an MHD fluid with a forced cosmic ray current and use the test-field method to determine the alpha effect and the turbulent magnetic diffusivity. As follows from DNS, the dynamics of the instability has the following stages: (i) in the early stage, the small-scale Bell instability that results in a production of small-scale turbulence is excited; (ii) in the intermediate stage, there is formation of larger-scale magnetic structures; (iii) finally, quasi-stationary large-scale turbulence is formed at a growth rate that is comparable to that expected from the dynamo instability, but its amplitude over much longer timescales remains unclear. The results of DNS are in good agreement with the theoretical estimates. It is suggested that this dynamo is what gives weakly magnetized relativistic shocks such as those from gamma ray bursts a macroscopic correlation length. It may also be important for large-scale magnetic field amplification associated with cosmic ray production and diffusive shock acceleration in supernova remnants (SNR) and blast waves from gamma ray bursts. Magnetic field amplification by Bell turbulence in SNR is found to be significant, but it is limited owing to the finite time available to the super-Alfvenicly expanding remnant. The effectiveness of the mechanisms is shown to be dependent on the shock velocity.

The negative effective magnetic pressure instability operates on scales encompassing many turbulent eddies and is here discussed in connection with the formation of active regions near the surface layers of the Sun. This instability is related to the negative contribution of turbulence to the mean magnetic pressure that causes the formation of large-scale magnetic structures. For an isothermal layer, direct numerical simulations and mean-field simulations of this phenomenon are shown to agree in many details in that their onset occurs at the same depth. This depth increases with increasing field strength, such that the maximum growth rate of this instability is independent of the field strength, provided the magnetic structures are fully contained within the domain. A linear stability analysis is shown to support this finding. The instability also leads to a redistribution of turbulent intensity and gas pressure that could provide direct observational signatures.

To understand the basic mechanism of the formation of magnetic flux concentrations, we determine by direct numerical simulations the turbulence contributions to the mean magnetic pressure in a strongly stratified isothermal layer with large plasma beta, where a weak uniform horizontal mean magnetic field is applied. The negative contribution of turbulence to the effective mean magnetic pressure is determined for strongly stratified forced turbulence over a range of values of magnetic Reynolds and Prandtl numbers. Small-scale dynamo action is shown to reduce the negative effect of turbulence on the effective mean magnetic pressure. However, the turbulence coefficients describing the negative effective magnetic pressure phenomenon are found to be converged for magnetic Reynolds numbers between 60 and 600, which is the largest value considered here. In all these models the turbulent intensity is arranged to be nearly independent of height, so the kinetic energy density decreases with height due to the decrease in density. In a second series of numerical experiments, the turbulent intensity increases with height such that the turbulent kinetic energy density is nearly independent of height. Turbulent magnetic diffusivity and turbulent pumping velocity are determined with the test-field method for both cases. The vertical profile of the turbulent magnetic diffusivity is found to agree with what is expected based on simple mixing length expressions. Turbulent pumping is shown to be down the gradient of turbulent magnetic diffusivity, but it is twice as large as expected. Corresponding numerical mean-field models are used to show that a large-scale instability can occur in both cases, provided the degree of scale separation is large enough and hence the turbulent magnetic diffusivity small enough.

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large-scale instability. In this study, hydromagnetic mean-field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two-dimensional. The growth rate is found to depend on a parameter beta_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter beta_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI.

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNS) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in formation of large-scale inhomogeneous magnetic structures from initial uniform magnetic field. Earlier investigations of this instability in DNS of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L^2/eta_t, where eta_t is the turbulent diffusivity and L is the scale of the domain. DNS are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number.

The phenomenon of turbulent thermal diffusion in temperature-stratified turbulence causing a non-diffusive turbulent flux of inertial and non-inertial particles in the direction of the turbulent heat flux is found using direct numerical simulations (DNS). In simulations with and without gravity, this phenomenon is found to cause a peak in the particle number density around the minimum of the mean fluid temperature for Stokes numbers less than 1, where the Stokes number is the ratio of particle Stokes time to turbulent Kolmogorov time at the viscous scale. Turbulent thermal diffusion causes the formation of large-scale inhomogeneities in the spatial distribution of inertial particles. The strength of this effect is maximum for Stokes numbers around unity, and decreases again for larger values. The dynamics of inertial particles is studied using Lagrangian modelling in forced temperature-stratified turbulence, whereas non-inertial particles and the fluid are described using DNS in an Eulerian framework.

We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed.

Using different analytical methods (the quasi-linear approach, the path-integral technique and tau-relaxation approximation) we develop a comprehensive mean-field theory for a pumping effect of the mean magnetic field in homogeneous non-rotating helical turbulence with imposed large-scale shear. The effective pumping velocity is proportional to the product of alpha effect and large-scale vorticity associated with the shear, and causes a separation of the toroidal and poloidal components of the mean magnetic field along the direction of the mean vorticity. We also perform direct numerical simulations of sheared turbulence in different ranges of hydrodynamic and magnetic Reynolds numbers and use a kinematic test-field method to determine the effective pumping velocity. The results of the numerical simulations are in agreement with the theoretical predictions.

In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed for simulating the upper solar convection zone and lower atmosphere. To this end we derived the energy budget equations, which include the effects of the subgrid-scale turbulence on the Lorentz-force, and implemented the new subgrid-scale turbulence model (TELF-Model) in a three-dimensional nonlinear MHD LES code. Using imposed initial vertical and horizontal uniform magnetic fields in LES with the TELF-Model, we have shown that the magnetic flux tubes formation is started when the initial mean magnetic field is larger than a threshold value (about 100 G). This is in agreement with the theoretical studies by Rogachevskii and Kleeorin (2007). We have determined the vertical profiles of the velocity and magnetic fluctuations, total MHD energy and anisotropy of turbulent magneto-convection, kinetic and current and cross helicities.

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.

A reduction of total mean turbulent pressure due to the presence of magnetic fields was previously shown to be a measurable effect in direct numerical simulations. However, in the studied parameter regime the formation of large-scale structures, as anticipated from earlier mean-field simulations, was not found. An analysis of the relevant mean-field parameter dependency and the parameter domain of interest is conducted in order to clarify this apparent discrepancy.

Is the scaling, \lambda \sim Rm^{1/2}, for the growth rate of small-scale dynamo instability at low magnetic Prandtl numbers and large magnetic Reynolds numbers, Rm, valid in the vicinity of the threshold? Our analysis and even numerical solution (Malyshkin and Boldyrev, 2010) of the dynamo equations for a Gaussian white-noise velocity field (the Kazantsev-Kraichnan model) imply that the answer is negative. Contrary to the claim by Malyshkin and Boldyrev (2010), there are two different asymptotics for the dynamo growth rate: in the vicinity of the threshold and far from the threshold.

While the rising flux tube paradigm is an elegant theory, its basic assumptions, thin flux tubes at the bottom of the convection zone with field strengths two orders of magnitude above equipartition, remain numerically unverified at best. As such, in recent years the idea of a formation of sunspots near the top of the convection zone has generated some interest. The presence of turbulence can strongly enhance diffusive transport mechanisms, leading to an effective transport coefficient formalism in the mean-field formulation. The question is what happens to these coefficients when the turbulence becomes anisotropic due to a strong large-scale mean magnetic field. It has been noted in the past that this anisotropy can also lead to highly non-diffusive behaviour. In the present work we investigate the formation of large-scale magnetic structures as a result of a negative contribution of turbulence to the large-scale effective magnetic pressure in the presence of stratification. In direct numerical simulations of forced turbulence in a stratified box, we verify the existence of this effect. This phenomenon can cause formation of large-scale magnetic structures even from initially uniform large-scale magnetic field.

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.

In this study we provide the first numerical demonstration of the effects of turbulence on the mean Lorentz force and the resulting formation of large-scale magnetic structures. Using three-dimensional direct numerical simulations (DNS) of forced turbulence we show that an imposed mean magnetic field leads to a decrease of the turbulent hydromagnetic pressure and tension. This phenomenon is quantified by determining the relevant functions that relate the sum of the turbulent Reynolds and Maxwell stresses with the Maxwell stress of the mean magnetic field. Using such a parameterization, we show by means of two-dimensional and three-dimensional mean-field numerical modelling that an isentropic density stratified layer becomes unstable in the presence of a uniform imposed magnetic field. This large-scale instability results in the formation of loop-like magnetic structures which are concentrated at the top of the stratified layer. In three dimensions these structures resemble the appearance of bipolar magnetic regions in the Sun. The results of DNS and mean-field numerical modelling are in good agreement with theoretical predictions. We discuss our model in the context of a distributed solar dynamo where active regions and sunspots might be rather shallow phenomena.

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.

We study magnetic field evolution in flows with fluctuating in time governing parameters in electrically conducting fluid. We use a standard mean-field approach to derive equations for large-scale magnetic field for the fluctuating ABC-flow as well as for the fluctuating Roberts flow. The derived mean-field dynamo equations have growing solutions with growth rate of the large-scale magnetic field which is not controlled by molecular magnetic diffusivity. Our study confirms the Zeldovich idea that the nonstationarity of the fluid flow may remove the obstacle in large-scale dynamo action of classic stationary flows.

Numerical simulations of forced turbulence in elongated shearing boxes are carried out to demonstrate that a nonhelical turbulence in conjunction with a linear shear can give rise to a mean-field dynamo. Exponential growth of magnetic field at scales larger than the outer (forcing) scale of the turbulence is found. Over a range of values of the shearing rate S spanning approximately two orders of magnitude, the growth rate of the magnetic field is proportional to the imposed shear, gamma ~ S, while the characteristic spatial scale of the field is l_b ~ S^(-1/2). The effect is quite general: earlier results for the nonrotating case by Yousef et al. 2008 (PRL 100, 184501) are extended to shearing boxes with Keplerian rotation; it is also shown that the shear dynamo mechanism operates both below and above the threshold for the fluctuation dynamo. The apparently generic nature of the shear dynamo effect makes it an attractive object of study for the purpose of understanding the generation of magnetic fields in astrophysical systems.

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.

Mechanisms of nonhelical large-scale dynamos (shear-current dynamo and effect of homogeneous kinetic helicity fluctuations with zero mean) in a homogeneous turbulence with large-scale shear are discussed. We have found that the shear-current dynamo can act even in random flows with small Reynolds numbers. However, in this case mean-field dynamo requires small magnetic Prandtl numbers (i.e., ${\rm Pm} < {\rm Pm}^{\rm cr}<1$). The threshold in the magnetic Prandtl number, ${\rm Pm}^{\rm cr} = 0.24$, is determined using second order correlation approximation (or first-order smoothing approximation) for a background random flow with a scale-dependent viscous correlation time $\tau_c=(\nu k^2)^{-1}$ (where $\nu$ is the kinematic viscosity of the fluid and $k$ is the wave number). For turbulent flows with large Reynolds numbers shear-current dynamo occurs for arbitrary magnetic Prandtl numbers. This dynamo effect represents a very generic mechanism for generating large-scale magnetic fields in a broad class of astrophysical turbulent systems with large-scale shear. On the other hand, mean-field dynamo due to homogeneous kinetic helicity fluctuations alone in a sheared turbulence is not realistic for a broad class of astrophysical systems because it requires a very specific random forcing of kinetic helicity fluctuations that contains, e.g., low-frequency oscillations.

The feasibility of a mean-field dynamo in nonhelical turbulence with superimposed linear shear is studied numerically in elongated shearing boxes. Exponential growth of magnetic field at scales much larger than the outer scale of the turbulence is found. The charateristic scale of the field is l_B ~ S^{-1/2} and growth rate is gamma ~ S, where S is the shearing rate. This newly discovered shear dynamo effect potentially represents a very generic mechanism for generating large-scale magnetic fields in a broad class of astrophysical systems with spatially coherent mean flows.

We study effects of kinetic helicity fluctuations in a turbulence with large-scale shear using two different approaches: the spectral tau-approximation and the second order correlation approximation (or first-order smoothing approximation). These two approaches demonstrate that homogeneous kinetic helicity fluctuations alone with zero mean value in a sheared homogeneous turbulence cannot cause large-scale dynamo. Mean-field dynamo can be possible when kinetic helicity fluctuations are inhomogeneous which cause a nonzero mean alpha effect in a sheared turbulence. On the other hand, shear-current effect can generate large-scale magnetic field even in a homogeneous nonhelical turbulence with large-scale shear. This effect was investigated previously for large hydrodynamic and magnetic Reynolds numbers. In this study we examine the threshold required for the shear-current dynamo versus Reynolds number. We demonstrate that there is no need for a developed inertial range in order to maintain the shear-current dynamo (e.g., the threshold in the Reynolds number is of the order of 1).

Magnetic fluctuations generated by a tangling of the mean magnetic field by velocity fluctuations are studied in a developed turbulent convection with large magnetic Reynolds numbers. We show that the energy of magnetic fluctuations depends on magnetic Reynolds number only when the mean magnetic field is smaller than $B_{eq} / 4 Rm^{1/4}$, where $B_{eq}$ is the equipartition mean magnetic field determined by the turbulent kinetic energy and Rm is magnetic Reynolds number. Generation of magnetic fluctuations in a turbulent convection with a nonzero mean magnetic field results in a decrease of the total turbulent pressure and may cause formation of the large-scale inhomogeneous magnetic structures even in an originally uniform mean magnetic field. This effect is caused by a negative contribution of the turbulent convection to the effective mean Lorentz force. The inhomogeneous large-scale magnetic fields are formed due to the excitation of the large-scale instability. The energy for this instability is supplied by the small-scale turbulent convection. The discussed effects might be useful for understanding the origin of the solar nonuniform magnetic fields, e.g., sunspots.

A review of recent studies on a new mechanism of generation of large-scale magnetic field in a sheared turbulent plasma is presented. This mechanism is associated with the shear-current effect which is related to the W x J-term in the mean electromotive force. This effect causes the generation of the large-scale magnetic field even in a nonrotating and nonhelical homogeneous sheared turbulent convection whereby the alpha effect vanishes. It is found that turbulent convection promotes the shear-current dynamo instability, i.e., the heat flux causes positive contribution to the shear-current effect. However, there is no dynamo action due to the shear-current effect for small hydrodynamic and magnetic Reynolds numbers even in a turbulent convection, if the spatial scaling for the turbulent correlation time is k^{-2}, where k is the small-scale wave number. We discuss here also the nonlinear mean-field dynamo due to the shear-current effect and take into account the transport of magnetic helicity as a dynamical nonlinearity. The magnetic helicity flux strongly affects the magnetic field dynamics in the nonlinear stage of the dynamo action. When the magnetic helicity flux is not small, the saturated level of the mean magnetic field is of the order of the equipartition field determined by the turbulent kinetic energy. The obtained results are important for elucidation of origin of the large-scale magnetic fields in astrophysical and cosmic sheared turbulent plasma.

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 study a nonlinear quenching of turbulent magnetic diffusion and effective drift velocity of large-scale magnetic field in a developed two-dimensional MHD turbulence at large magnetic Reynolds numbers. We show that transport of the mean-square magnetic potential strongly changes quenching of turbulent magnetic diffusion. In particularly, the catastrophic quenching of turbulent magnetic diffusion does not occur for the large-scale magnetic fields $B \gg B_{\rm eq} / \sqrt{\rm Rm}$ when a divergence of the flux of the mean-square magnetic potential is not zero, where $B_{\rm eq}$ is the equipartition mean magnetic field determined by the turbulent kinetic energy and Rm is the magnetic Reynolds number. In this case the quenching of turbulent magnetic diffusion is independent of magnetic Reynolds number. The situation is similar to three-dimensional MHD turbulence at large magnetic Reynolds numbers whereby the catastrophic quenching of the alpha effect does not occur when a divergence of the flux of the small-scale magnetic helicity is not zero.

The shear-current effect in a nonrotating homogeneous turbulent convection with a large-scale constant shear is studied. The large-scale velocity shear causes anisotropy of turbulent convection, which produces the mean electromotive force $\bec{\cal E}^{(W)} \propto {\bf W} {\bf \times} {\bf J}$ and the mean electric current along the original mean magnetic field, where ${\bf W}$ is the background mean vorticity due to the shear and ${\bf J}$ is the mean electric current. This results in a large-scale dynamo even in a nonrotating and nonhelical homogeneous sheared turbulent convection, whereby the $\alpha$ effect vanishes. It is found that turbulent convection promotes the shear-current dynamo instability, i.e., the heat flux causes positive contribution to the shear-current effect. However, there is no dynamo action due to the shear-current effect for small hydrodynamic and magnetic Reynolds numbers even in a turbulent convection, if the spatial scaling for the turbulent correlation time is $\tau(k) \propto k^{-2}$, where $k$ is the small-scale wave number.

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.

The nonlinear mean-field dynamo due to a shear-current effect in a nonhelical homogeneous turbulence with a mean velocity shear is discussed. The transport of magnetic helicity as a dynamical nonlinearity is taken into account. The shear-current effect is associated with the ${\bf W} {\bf \times} {\bf J}$ term in the mean electromotive force, where ${\bf W}$ is the mean vorticity due to the large-scale shear motions and ${\bf J}$ is the mean electric current. This effect causes the generation of large-scale magnetic field in a turbulence with large hydrodynamic and magnetic Reynolds numbers. The dynamo action due to the shear-current effect depends on the spatial scaling of the correlation time $\tau(k)$ of the background turbulence, where $k$ is the wave number. For Kolmogorov scaling, $\tau(k) \propto k^{-2/3}$, the dynamo instability occurs, while when $\tau(k) \propto k^{-2}$ (small hydrodynamic and magnetic Reynolds numbers) there is no the dynamo action in a sheared nonhelical turbulence. The magnetic helicity flux strongly affects the magnetic field dynamics in the nonlinear stage of the dynamo action. Numerical solutions of the nonlinear mean-field dynamo equations which take into account the shear-current effect, show that if the magnetic helicity flux is not small, the saturated level of the mean magnetic field is of the order of the equipartition field determined by the turbulent kinetic energy. Turbulence with a large-scale velocity shear is a universal feature in astrophysics, and the obtained results can be important for elucidation of origin of the large-scale magnetic fields in astrophysical sheared turbulence.

We determine the nonlinear drift velocities of the mean magnetic field and nonlinear turbulent magnetic diffusion in a turbulent convection. We show that the nonlinear drift velocities are caused by the three kinds of the inhomogeneities, i.e., inhomogeneous turbulence; the nonuniform fluid density and the nonuniform turbulent heat flux. The inhomogeneous turbulence results in the well-known turbulent diamagnetic and paramagnetic velocities. The nonlinear drift velocities of the mean magnetic field cause the small-scale magnetic buoyancy and magnetic pumping effects in the turbulent convection. These phenomena are different from the large-scale magnetic buoyancy and magnetic pumping effects which are due to the effect of the mean magnetic field on the large-scale density stratified fluid flow. The small-scale magnetic buoyancy and magnetic pumping can be stronger than these large-scale effects when the mean magnetic field is smaller than the equipartition field. We discuss the small-scale magnetic buoyancy and magnetic pumping effects in the context of the solar and stellar turbulent convection. We demonstrate also that the nonlinear turbulent magnetic diffusion in the turbulent convection is anisotropic even for a weak mean magnetic field. In particular, it is enhanced in the radial direction. The magnetic fluctuations due to the small-scale dynamo increase the turbulent magnetic diffusion of the toroidal component of the mean magnetic field, while they do not affect the turbulent magnetic diffusion of the poloidal field.

An effect of the differential rotation on the nonlinear electromotive force in MHD turbulence is found. It includes a nonhelical $\alpha$ effect which is caused by a differential rotation, and it is independent of a hydrodynamic helicity. There is no quenching of this effect contrary to the quenching of the usual $\alpha$ effect caused by a hydrodynamic helicity. The nonhelical $\alpha$ effect vanishes when the rotation is constant on the cylinders which are parallel to the rotation axis. The mean differential rotation creates also the shear-current effect which changes its sign with the nonlinear growth of the mean magnetic field. However, there is no quenching of this effect. These phenomena determine the nonlinear evolution of the mean magnetic field. An effect of a uniform rotation on the nonlinear electromotive force is also studied. A nonlinear theory of the ${\bf \Omega} {\bf \times} \bar{\bf J}$ effect is developed, and the quenching of the hydrodynamic part of the usual $\alpha$ effect which is caused by a uniform rotation and inhomogeneity of turbulence, is found. Other contributions of a uniform rotation to the nonlinear electromotive force are also determined. All these effects are studied using the spectral $\tau$ approximation (the third-order closure procedure). An axisymmetric mean-field dynamo in the spherical and cylindrical geometries is considered. The nonlinear saturation mechanism based on the magnetic helicity evolution is discussed. It is shown that this universal mechanism is nearly independent of the form of the flux of magnetic helicity, and it requires only a nonzero flux of magnetic helicity. Astrophysical applications of these effects are discussed.

A new mechanism of generation of large-scale magnetic fields in colliding protogalactic clouds and merging protostellar clouds is discussed. Interaction of the colliding clouds produces large-scale shear motions which are superimposed on small-scale turbulence. Generation of the large-scale magnetic field is due to a ''shear-current" effect (or "vorticity-current" effect), and the mean vorticity is caused by the large-scale shear motions of colliding clouds. This effect causes the generation of the mean magnetic field even in a nonrotating and nonhelical homogeneous turbulence. There is no quenching of the nonlinear shear-current effect contrary to the quenching of the nonlinear alpha effect, the nonlinear turbulent magnetic diffusion, etc. During the nonlinear growth of the mean magnetic field, the shear-current effect only changes its sign at some value of the mean magnetic field which determines the level of the saturated mean magnetic field. Numerical study shows that the saturated level of the mean magnetic field is of the order of the equipartition field determined by the turbulent kinetic energy. The estimated large-scale magnetic field for merging protogalactic clouds is about several microgauss, and for merging protostellar clouds is of the order of several tenth of microgauss.

A fraction of solar active regions are observed to have current helicity of a sign that contradicts the polarity law for magnetic helicity; this law corresponds to the well-known Hale polarity law for sunspots. A significant excess of active regions with the "wrong" sign of helicity is seen to occur just at the beginning of the cycle. We compare these observations with predictions from a dynamo model based on principles of helicity conservation, discussed by Zhang et al. (2006). This model seems capable of explaining only a fraction of the regions with the wrong sign of the helicity. We attribute the remaining excess to additional current helicity production from the twisting of rising magnetic flux tubes, as suggested by Choudhuri et al. (2004). We estimate the relative contributions of this effect and that connected with the model based on magnetic helicity conservation.

We studied the effect of the turbulent heat flux on the Reynolds stresses in a rotating turbulent convection. To this end we solved a coupled system of dynamical equations which includes the equations for the Reynolds stresses, the entropy fluctuations and the turbulent heat flux. We used a spectral $\tau$ approximation in order to close the system of dynamical equations. We found that the ratio of the contributions to the Reynolds stresses caused by the turbulent heat flux and the anisotropic eddy viscosity is of the order of $\sim 10 (L_\rho / l_0)^2$, where $l_{0}$ is the maximum scale of turbulent motions and $L_\rho$ is the fluid density variation scale. This effect is crucial for the formation of the differential rotation and should be taken into account in the theories of the differential rotation of the Sun, stars and planets. In particular, we demonstrated that this effect may cause the differential rotation which is comparable with the typical solar differential rotation.

We continue our attempt to connect observational data on current helicity in solar active regions with solar dynamo models. In addition to our previous results about temporal and latitudinal distributions of current helicity (Kleeorin et al. 2003), we argue that some information concerning the radial profile of the current helicity averaged over time and latitude can be extracted from the available observations. The main feature of this distribution can be presented as follows. Both shallow and deep active regions demonstrate a clear dominance of one sign of current helicity in a given hemisphere during the whole cycle. Broadly speaking, current helicity has opposite polarities in the Northern and Southern hemispheres, although there are some active regions that violate this polarity rule. The relative number of active regions violating the polarity rule is significantly higher for deeper active regions. A separation of active regions into `shallow', `middle' and `deep' is made by comparing their rotation rate and the helioseismic rotation law. We use a version of Parker's dynamo model in two spatial dimensions, that employs a nonlinearity based on magnetic helicity conservation arguments. The predictions of this model about the radial distribution of solar current helicity appear to be in remarkable agreement with the available observational data; in particular the relative volume occupied by the current helicity of "wrong" sign grows significantly with the depth.

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.

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.

The nonlinear theory of a "shear-current" effect in a nonrotating and nonhelical homogeneous turbulence with an imposed mean velocity shear is developed. The ''shear-current" effect is associated with the $\bar{\bf W} {\bf \times} \bar{\bf J}$-term in the mean electromotive force and causes the generation of the mean magnetic field even in a nonrotating and nonhelical homogeneous turbulence (where $\bar{\bf W}$ is the mean vorticity and $\bar{\bf J}$ is the mean electric current). It is found that there is no quenching of the nonlinear "shear-current" effect contrary to the quenching of the nonlinear $\alpha$-effect, the nonlinear turbulent magnetic diffusion, etc. During the nonlinear growth of the mean magnetic field, the ''shear-current" effect only changes its sign at some value $\bar{\bf B}_\ast$ of the mean magnetic field. The magnitude $\bar{\bf B}_\ast$ determines the level of the saturated mean magnetic field which is less than the equipartition field. It is shown that the background magnetic fluctuations due to the small-scale dynamo enhance the "shear-current" effect, and reduce the magnitude $\bar{\bf B}_\ast$. When the level of the background magnetic fluctuations is larger than 1/3 of the kinetic energy of the turbulence, the mean magnetic field can be generated due to the "shear-current" effect for an arbitrary exponent of the energy spectrum of the velocity fluctuations.

We study a simple model for the solar dynamo in the framework of the Parker migratory dynamo, with a nonlinear dynamo saturation mechanism based on magnetic helicity conservation arguments. We find a parameter range in which the model demonstrates a cyclic behaviour with properties similar to that of Parker dynamo with the simplest form of algebraic alpha-quenching. We compare the nonlinear current helicity evolution in this model with data for the current helicity evolution obtained during 10 years of observations at the Huairou Solar Station of China. On one hand, our simulated data demonstrate behaviour comparable with the observed phenomenology, provided that a suitable set of governing dynamo parameters is chosen. On the other hand, the observational data are shown to be rich enough to reject some other sets of governing parameters. We conclude that, in spite of the very preliminary state of the observations and the crude nature of the model, the idea of using observational data to constrain our ideas concerning magnetic field generation in the framework of the solar dynamo appears promising.

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.

An effect of sheared large-scale motions on a mean electromotive force in a nonrotating turbulent flow of a conducting fluid is studied. It is demonstrated that in a homogeneous divergence-free turbulent flow the alpha-effect does not exist, however a mean magnetic field can be generated even in a nonrotating turbulence with an imposed mean velocity shear due to a new ''shear-current" effect. A contribution to the electromotive force related with the symmetric parts of the gradient tensor of the mean magnetic field (the kappa-effect) is found in a nonrotating turbulent flows with a mean shear. The kappa-effect and turbulent magnetic diffusion reduce the growth rate of the mean magnetic field. It is shown that a mean magnetic field can be generated when the exponent of the energy spectrum of the background turbulence (without the mean velocity shear) is less than 2. The ''shear-current" effect was studied using two different methods: the Orszag third-order closure procedure and the stochastic calculus. 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.

We have extended our previous mean-field galactic dynamo model which included algebraic and dynamic alpha nonlinearities (Kleeorin et al., A&A, v. 387, 453, 2002), to include also a quenching of turbulent diffusivity. We readily obtain equilibrium states for the large-scale magnetic field in the local disc dynamo model, and these fields have strengths that are comparable to the equipartition field strength. We find that the algebraic nonlinearity alone (i.e. quenching of both the alpha effect and turbulent magnetic diffusion) cannot saturate the growth of the mean magnetic field; only the combined effect of algebraic and dynamic nonlinearities can limit the growth of the mean magnetic field. However, in contrast to our earlier work without quenching of the turbulent diffusivity, we cannot now find satisfactory solutions in the no-z approximation to the axisymmetric galactic dynamo problem.

An effect of rotation on a developed turbulent stratified convection is studied. Dependencies of the hydrodynamic helicity, the alpha-tensor and the effective drift velocity of the mean magnetic field on the rate of rotation and an anisotropy of turbulent convection are found. It is shown that in an anisotropic turbulent convection the alpha-effect can change its sign depending on the rate of rotation. Different properties of the effective drift velocity of the mean magnetic field in a rotating turbulent convection are found, e.g., there is a difference in the effective drift velocities for the toroidal and poloidal magnetic fields; a toroidal effective drift velocity can play a role of an additional differential rotation. These effects and an effect of a nonzero divergence of the effective drift velocity of the toroidal magnetic field on a magnetic dynamo are studied. Astrophysical applications of the obtained results are discussed.

The mean electromotive force in the framework of mean--field magnetohydrodynamics is found with a particular attention to the effect of a mean rotation of the fluid on the mean electromotive force. Anisotropy of the turbulence due to gradients of its intensity or its helicity are admitted. The mean magnetic field is considered to be weak enough to exclude quenching effects. Magnetic fluctuations with a zero mean magnetic field are taken into account. The Omega x J -effect works in the same way with velocity fluctuations and magnetic fluctuations. A contribution to the electromotive force connected with the symmetric parts of the gradient tensor of the mean magnetic field was found in the case of rotation, resulting from velocity or magnetic fluctuations. The implications of the results for the mean electromotive force for mean--field dynamo models are discussed with special emphasis to dynamos working without alpha-effect.

The nonlinear coefficients defining the mean electromotive force (i.e., the nonlinear turbulent magnetic diffusion, the nonlinear effective velocity, the nonlinear kappa-tensor, etc.) are calculated for an anisotropic turbulence. A particular case of an anisotropic background turbulence (i.e., the turbulence with zero mean magnetic field) with one preferential direction is considered. It is shown that the toroidal and poloidal magnetic fields have different nonlinear turbulent magnetic diffusion coefficients. It is demonstrated that even for a homogeneous turbulence there is a nonlinear effective velocity which exhibits diamagnetic or paramagnetic properties depending on anisotropy of turbulence and level of magnetic fluctuations in the background turbulence. Analysis shows that an anisotropy of turbulence strongly affects the nonlinear mean electromotive force. Two types of nonlinearities (algebraic and dynamic) are also discussed. The algebraic nonlinearity implies a nonlinear dependence of the mean electromotive force on the mean magnetic field. The dynamic nonlinearity is determined by a differential equation for the magnetic part of the alpha-effect. It is shown that for the alpha-Omega axisymmetric dynamo the algebraic nonlinearity alone cannot saturate the dynamo generated mean magnetic field while the combined effect of the algebraic and dynamic nonlinearities limits the mean magnetic field growth. Astrophysical applications of the obtained results are discussed.

The evolution of the magnetic helicity tensor for a nonzero mean magnetic field and for large magnetic Reynolds numbers in an anisotropic turbulence is studied. It is shown that the isotropic and anisotropic parts of the magnetic helicity tensor have different characteristic times of evolution. The time of variation of the isotropic part of the magnetic helicity tensor is much larger than the correlation time of the turbulent velocity field. The anisotropic part of the magnetic helicity tensor changes for the correlation time of the turbulent velocity field. The mean turbulent flux of the magnetic helicity is calculated as well. It is shown that even a small anisotropy of turbulence strongly modifies the flux of the magnetic helicity. It is demonstrated that the tensor of the magnetic part of the alpha-effect for weakly inhomogeneous turbulence is determined only by the isotropic part of the magnetic helicity tensor. Astrophysical applications of the obtained results are discussed.

It is demonstrated that the two-point correlation function of the magnetic helicity in the case of zero mean magnetic field has anomalous scalings for both, compressible and incompressible turbulent helical fluid flow. The magnetic helicity in the limit of very high electrical conductivity is conserved. This implies that the two-point correlation function of the conserved property does not necessarily have normal scaling. The reason for the anomalous scalings of the magnetic helicity correlation function is that the magnetic field in the equation for the two-point correlation function of the magnetic helicity plays a role of a pumping with anomalous scalings. It is shown also that when magnetic fluctuations with zero mean magnetic field are generated the magnetic helicity is very small even if the hydrodynamic helicity is large. Astrophysical applications of the obtained results are discussed.

We demonstrate that the inclusion of the helicity flux in the magnetic helicity balance in the nonlinear stage of galactic dynamo action results in a radical change in the magnetic field dynamics. The equilibrium value of the large-scale magnetic field is then approximately the equipartition level. This is in contrast to the situation without the flux of helicity, when the magnetic helicity is conserved locally, which leads to substantially subequipartition values for the equilibrium large-scale magnetic field.

Magnetic fluctuations with a zero mean field in a random flow with a finite correlation time and a small yet finite magnetic diffusion are studied. Equation for the second-order correlation function of a magnetic field is derived. This equation comprises spatial derivatives of high orders due to a non-local nature of magnetic field transport in a random velocity field with a finite correlation time. For a random Gaussian velocity field with a small correlation time the equation for the second-order correlation function of the magnetic field is a third-order partial differential equation. For this velocity field and a small magnetic diffusion with large magnetic Prandtl numbers the growth rate of the second moment of magnetic field is estimated. The finite correlation time of a turbulent velocity field causes an increase of the growth rate of magnetic fluctuations. It is demonstrated that the results obtained for the cases of a small yet finite magnetic diffusion and a zero magnetic diffusion are different. Astrophysical applications of the obtained results are discussed.

We consider the magnetic helicity balance for the galactic dynamo in the framework of the local dynamo problem, as well as in the no-z model (which includes explicitly the radial distribution of the magnetic fields). When calculating the magnetic helicity balance we take into account the redistribution of the small-scale and large-scale magnetic fields between the magnetic helicities, as well as magnetic helicity transport and diffusion due to small-scale turbulence. We demonstrate that the magnetic helicity flux through the galactic disc boundaries leads to a steady-state magnetic field with magnetic energy comparable to the equipartition energy of the turbulent motions of the interstellar medium. If such flux is ignored, the steady-state magnetic field is found to be much smaller than the equipartition field. The total magnetic helicity flux through the boundaries consists of both an advective flux and a diffusive flux. The exact ratio of these contributions seems not to be crucial for determining the strength of the steady-state magnetic field and its structure. However at least some diffusive contribution is needed to smooth the magnetic helicity profile near to the disc boundaries. The roles of various transport coefficients for magnetic helicity are investigated, and the values which lead to magnetic field configurations comparable with those observed are determined.