Source author record

I. Rogachevskii

I. Rogachevskii appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2021arXiv

Energy and flux budget closure theory for passive scalar in stably stratified turbulence

The energy and flux budget (EFB) closure theory for a passive scalar (non-buoyant and non-inertial particles or gaseous admixtures) is developed for stably stratified turbulence. The physical background of the EFB turbulence closures is based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, as well as the turbulent flux of particles. The EFB turbulence closure is designed for stratified geophysical flows from neutral to very stable stratification and it implies that turbulence is maintained by the velocity shear at any stratification. In a steady-state, expressions for the turbulent flux of passive scalar and the anisotropic non-symmetric turbulent diffusion tensor are derived, and universal flux Richardson number dependencies of the components of this tensor are obtained. The diagonal component in the vertical direction of the turbulent diffusion tensor is suppressed by strong stratification, while the diagonal components in the horizontal directions are not suppressed, and they are dominant in comparison with the other components of turbulent diffusion tensor. This implies that any initially created strongly inhomogeneous particle cloud is evolved into a thin pancake in horizontal plane with very slow increase of its thickness in the vertical direction. The turbulent Schmidt number increases linearly with the gradient Richardson number. Considering the applications of these results to the atmospheric boundary-layer turbulence, the theoretical relationships are derived which allow to determine the turbulent diffusion tensor as a function of the vertical coordinate measured in the units of the local Obukhov length scale. The obtained relations are potentially useful in modelling applications of particle dispersion in the atmospheric boundary-layer turbulence and free atmosphere turbulence.

preprint2021arXiv

Turbulent transport of radiation in the solar convective zone

A turbulent transport of radiation in the solar convective zone is investigated. The mean-field equation for the irradiation intensity is derived. It is shown that due to the turbulent effects, the effective penetration length of radiation can be increased in several times in comparison with the mean penetration length of radiation (defined as an inverse mean absorption coefficient). Using the model of the solar convective zone based on the mixing length theory, where the mean penetration length of radiation is usually much smaller than the turbulent correlation length, it is demonstrated that the ratio of the effective penetration length to the mean penetration length of radiation increases in 2.5 times in the vicinity of the solar surface. The main reason are the compressibility effects that become important in the vicinity of the solar surface where temperature and density fluctuations increase towards the solar surface, enhancing fluctuations of the radiation absorption coefficient and increasing the effective penetration length of radiation.

preprint2020arXiv

Compressibility effects in a turbulent transport of temperature field

Compressibility effects in a turbulent transport of temperature field are investigated applying the quasi-linear approach for small Péclet numbers and the spectral $τ$ approach for large Péclet numbers. Compressibility of a fluid flow reduces the turbulent diffusivity of the mean temperature field similarly to that for particle number density and magnetic field. However, expressions for the turbulent diffusion coefficient for the mean temperature field in a compressible turbulence are different from those for the mean particle number density and the mean magnetic field. Combined effect of compressibility and inhomogeneity of turbulence causes an increase of the mean temperature in the regions with more intense velocity fluctuations due to a turbulent pumping. Formally, this effect is similar to a phenomenon of compressible turbophoresis found previously [J. Plasma Phys. {\bf 84}, 735840502 (2018)] for non-inertial particles or gaseous admixtures. Gradient of the mean fluid pressure results in an additional turbulent pumping of the mean temperature field. The latter effect is similar to turbulent barodiffusion of particles and gaseous admixtures. Compressibility of a fluid flow also causes a turbulent cooling of the surrounding fluid due to an additional sink term in the equation for the mean temperature field. There is no analog of this effect for particles.

preprint2020arXiv

Detection of standing internal gravity waves in experiments with convection over a wavy heated wall

Convection over a wavy heated bottom wall in the air flow has been studied in experiments with the Rayleigh number $\sim 10^8$. It is shown that the mean temperature gradient in the flow core inside a large-scale circulation is directed upward, that corresponds to the stably stratified flow. In the experiments with a wavy heated bottom wall, we detect large-scale standing internal gravity waves excited in the regions with the stably stratified flow. The wavelength and the period of these waves are much larger than the turbulent spatial and time scales, respectively. In particular, the frequencies of the observed large-scale waves vary from 0.006 Hz to 0.07 Hz, while the turbulent time in the integral scale is about 0.5 s. The measured spectra of these waves contains several localized maxima, that implies an existence of waveguide resonators for the large-scale standing internal gravity waves. For comparisons, experiments with convection over a smooth plane bottom wall at the same mean temperature difference between bottom and upper walls have been also conducted. In these experiments various locations with a stably stratified flow are also found and the large-scale standing internal gravity waves are observed in these regions.

preprint2020arXiv

Estimates of current helicity and tilt of solar active regions and Joy's law

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.

preprint2020arXiv

The mean tilt of sunspot bipolar regions: theory, simulations and comparison with observations

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.

preprint2018arXiv

Mean-field theory of differential rotation in density stratified turbulent convection

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.

preprint2018arXiv

Nonlinear mean-field dynamo and prediction of solar activity

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.

preprint2017arXiv

Multipoint Radiation Induced Ignition of Dust Explosions: Turbulent Clustering of Particles and Increased Transparency

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.

preprint2016arXiv

The dynamics of Wolf numbers based on nonlinear dynamo with magnetic helicity: comparisons with observations

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.

preprint2015arXiv

A hierarchy of energy- and flux-budget (EFB) turbulence closure models for stably stratified geophysical flows

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.

preprint2015arXiv

Acceleration of raindrops formation due to tangling-clustering instability in turbulent stratified atmosphere

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 $μ$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.

preprint2015arXiv

Magnetic flux concentrations from turbulent stratified convection

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

preprint2014arXiv

Intense bipolar structures from stratified helical dynamos

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 $δ$ 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.

preprint2014arXiv

Magnetic flux concentrations in a polytropic atmosphere

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.

preprint2014arXiv

Mean-field and direct numerical simulations of magnetic flux concentrations from vertical field

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.

preprint2014arXiv

Sheared stably stratified turbulence and large-scale waves in a lid driven cavity

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.

preprint2013arXiv

Competition of rotation and stratification in flux concentrations

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 $τ$ 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.}

preprint2012arXiv

Current helicity of active regions as a tracer of large-scale solar magnetic helicity

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.

preprint2012arXiv

Experimental study of temperature fluctuations in forced stably stratified turbulent flows

We study experimentally temperature fluctuations in stably stratified forced turbulence in air flow. In the experiments with an imposed vertical temperature gradient, the turbulence is produced by two oscillating grids located nearby the side walls of the chamber. Particle Image Velocimetry is used to determine the turbulent and mean velocity fields, and a specially designed temperature probe with sensitive thermocouples is employed to measure the temperature field. We found that the ratio $\big[(\ell_x \nabla_x \bar{T})^2 + (\ell_y \nabla_y \bar{T})^2 + (\ell_z \nabla_z \bar{T})^2\big] / <θ^2 >$ is nearly constant, is independent of the frequency of the grid oscillations and has the same magnitude for both, stably and unstably stratified turbulent flows, where $\ell_i$ are the integral scales of turbulence along $x, y, z$ directions, $\bar{T}$ and $θ$ are the mean and turbulent fluctuations components of the fluid temperature. We demonstrated that for large frequencies of the grid oscillations, the temperature field can be considered as a passive scalar, while for smaller frequencies the temperature field behaves as an active field. The theoretical predictions based on the budget equations for turbulent kinetic energy, turbulent potential energy ($\propto <θ^2 >)$ and turbulent heat flux, are in a good agreement with the experimental results. Detailed comparison with experimental results obtained previously in unstably stratified forced turbulence is performed.

preprint2012arXiv

Flux concentrations in turbulent convection

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.

preprint2012arXiv

Growth rate of small-scale dynamo at low magnetic Prandtl numbers

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, $λ\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, $λ\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.

preprint2012arXiv

Negative effective magnetic pressure in turbulent convection

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.

preprint2012arXiv

New scaling for the alpha effect in slowly rotating turbulence

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.

preprint2012arXiv

Properties of the negative effective magnetic pressure instability

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.

preprint2012arXiv

Rotational effects on the negative magnetic pressure instability

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.

preprint2011arXiv

Detection of turbulent thermal diffusion of particles in numerical simulations

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.

preprint2011arXiv

LES of turbulent convection in solar-type stars and formation of large-scale magnetic structures

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.

preprint2011arXiv

Pumping velocity in homogeneous helical turbulence with shear

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.

preprint2011arXiv

Transition phenomena in unstably stratified turbulent flows

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énard 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] / < θ^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 θ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.

preprint2010arXiv

Tangling clustering of inertial particles in stably stratified turbulence

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.

preprint2010arXiv

The negative magnetic pressure effect in stratified turbulence

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.

preprint2009arXiv

Effect of large-scale coherent structures on turbulent convection

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.

preprint2009arXiv

Energy- and flux-budget turbulence closure model for stably stratified flows. Part II: the role of internal gravity waves

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.

preprint2009arXiv

Mixing at the external boundary of a submerged turbulent jet

We study experimentally and theoretically mixing at the external boundary of a submerged turbulent jet. In the experimental study we use Particle Image Velocimetry and an Image Processing Technique based on the analysis of the intensity of the Mie scattering to determine the spatial distribution of tracer particles. An air jet is seeded with the incense smoke particles which are characterized by large Schmidt number and small Stokes number. We determine the spatial distributions of the jet fluid characterized by a high concentration of the particles and of the ambient fluid characterized by a low concentration of the tracer particles. In the data analysis we use two approaches, whereby one approach is based on the measured phase function for the study of the mixed state of two fluids. The other approach is based on the analysis of the two-point second-order correlation function of the particle number density fluctuations generated by tangling of the gradient of the mean particle number density by the turbulent velocity field. This gradient is formed at the external boundary of a submerged turbulent jet. We demonstrate that PDF of the phase function of a jet fluid penetrating into an external flow and the two-point second-order correlation function of the particle number density do not have universal scaling and cannot be described by a power-law function. The theoretical predictions made in this study are in a qualitative agreement with the obtained experimental results.

preprint2009arXiv

Turbulent Diffusion and Turbulent Thermal Diffusion of Aerosols in Stratified Atmospheric Flows

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.

preprint2007arXiv

Clustering of Aerosols in Atmospheric Turbulent Flow

A mechanism of formation of small-scale inhomogeneities in spatial distributions of aerosols and droplets associated with clustering instability in the atmospheric turbulent flow is discussed. The particle clustering is a consequence of a spontaneous breakdown of their homogeneous space distribution due to the clustering instability, and is caused by a combined effect of the particle inertia and a finite correlation time of the turbulent velocity field. In this paper a theoretical approach proposed in Phys. Rev. E 66, 036302 (2002) is further developed and applied to investigate the mechanisms of formation of small-scale aerosol inhomogeneities in the atmospheric turbulent flow. The theory of the particle clustering instability is extended to the case when the particle Stokes time is larger than the Kolmogorov time scale, but is much smaller than the correlation time at the integral scale of turbulence. We determined the criterion of the clustering instability for the Stokes number larger than 1. We discussed applications of the analyzed effects to the dynamics of aerosols and droplets in the atmospheric turbulent flow.

preprint2007arXiv

Energy- and flux-budget (EFB) turbulence closure model for the stably stratified flows. Part I: Steady-state, homogeneous regimes

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

preprint2004arXiv

Turbulent thermal diffusion of aerosols in geophysics and laboratory experiments

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.

preprint2002arXiv

Nonlinear Magnetic Diffusion and Magnetic Helicity Transport in Galactic Dynamos

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.