Helicity generation in turbulent MHD flows

We study turbulent flow of a conducting liquid in a uniform external magnetic field. It is shown that intense helicity generation is possible in the presence of a mean shear flow. It is noted that even though the mean helicity of the initial flow can be zero, the presence of internal topological structure of the flow, for example the presence of helicity of different signs at different scales, is nevertheless necessary for helicity generation. Zh. éksp. Teor. Fiz. 114, 946–955 (September 1998)     

On the role of helicity in complex fluid flows

Results of direct numerical simulations of the Taylor-Green vortex are analysed by conditional sampling. In regions of small energy dissipation, there are tendencies for (1) velocity, u, and vorticity, ω, to be aligned and (2) vorticity and curl of vorticity, ? × ω, to be nearly orthogonal. The fields of dissipation, enstrophy, turbulence production, and vortex stretching exhibit a striking similarity in their spatial structure.     

Effect of hyperdiffusivity on turbulent dynamos with helicity

In numerical studies of turbulence, hyperviscosity is often used as a tool to extend the inertial subrange and to reduce the dissipative subrange. By analogy, hyperdiffusivity (or hyperresistivity) is sometimes used in magnetohydrodynamics. The underlying assumption is that only the small scales are affected by this manipulation. In the present paper, possible side effects on the evolution of the large-scale magnetic field are investigated. It is found that for turbulent flows with helicity, hyperdiffusivity causes the dynamo-generated magnetic field to saturate at a higher level than normal diffusivity. This result is successfully interpreted in terms of magnetic helicity conservation, which also predicts that full saturation is reached only after a time comparable to the large-scale magnetic (hyper)diffusion time.     

The role of directionality on the structure and dynamics of strongly anisotropic turbulent flows

Strong anisotropy in turbulent flows may be induced by body forces, Coriolis, buoyancy, Lorentz, and/or by large-scale gradients. These effects combined to the redistribution pressure terms are first identified by an angle dependence of the wave vector in Fourier space, the directionality. The resulting anisotropic structure is not taken into account in classical phenomenological theory, using essentially ‘isotropised’ dimensional analysis. Besides, it is generally hidden in practical engineering models by means of tuned constants, which may vary if the flow changes of nature. In this paper, different examples of anisotropic turbulence are revisited and compared to each other in order to shed light on fundamental aspects of this specific turbulence. To begin with, flows without energy production like rotating turbulence are considered. In this case, isotropy is broken by mean of third-order correlations in the equations. These correlations quantify the interscale energy transfer, and must be investigated at three-point, or triad by triad in Fourier space. This allows to account for the role of typical anisotropic frequency 2Ωcos θ _k , with θ _k the angle of the wave vector to the axis of rotation, and to simultaneously restore the role of phase coherence. We pursue the discussion with a second flow case, with production, quasi-static magnetohydrodynamics. This illustrates turbulence forced towards two-dimensional structure by an explicit Ohmic dissipation term. Linear dynamics displays an angle (called Moreau, or Shebalin) capable of reflecting the basic anisotropy in models as simple as . In the final phase of transition towards 2D structure, however, dynamics are essentially driven by third-order velocity correlations, and both successive linear and nonlinear phases yield counter-intuitive anisotropic results. The last case considered here is the turbulent mixing induced by a Rayleigh-Taylor instability. It is shown that anisotropy plays a central role in the dynamics of the mixing zone by means of an angular dimensionality parameter similar to the Moreau angle but for the density field, and appearing in a global model of buoyancy-drag equation.     

Nonlinear current helicity fluxes in turbulent dynamos and alpha quenching

Large scale dynamos produce small scale current helicity as a waste product that quenches the large scale dynamo process (alpha effect). This quenching can be catastrophic (i.e., intensify with magnetic Reynolds number) unless one has fluxes of small scale magnetic (or current) helicity out of the system. We derive the form of helicity fluxes in turbulent dynamos, taking also into account the nonlinear effects of Lorentz forces due to fluctuating fields. We confirm the form of an earlier derived magnetic helicity flux term, and also show that it is not renormalized by the small scale magnetic field, just like turbulent diffusion. Additional nonlinear fluxes are identified, which are driven by the anisotropic and antisymmetric parts of the magnetic correlations. These could provide further ways for turbulent dynamos to transport out small scale magnetic helicity, so as to avoid catastrophic quenching.     

Helical turbulent spectra in the presence of energy and helicity fluxes

On the basis of asymptotic approach, the behavior of helicity spectra in turbulent flows has been studied. The obtained expressions relate these spectra to energy spectra on the condition that energy flux ε and helicity flux η are nonzero both in the inertial and dissipative intervals.     

Large-scale amplification of a magnetic field by turbulent helicity fluctuations

In this paper the procedure of large-scale averaging of the magnetic-field diffusion equation with the α-term curlα(r,t)B(r,t) is used to show that a nonuniform distribution of the turbulent helicity fluctuations (more precisely, the fluctuations of the coefficient α) with a zero average value gives rise to large-scale amplification of the initial magnetic field. A detailed study is carried out of the dependence of the resulting large-scale α effect on the characteristics of the correlator 〈〈α(r, t)α(r″,t″)〉〉 in a rotating medium with a nonuniform distribution of the angular velocity ω=ω(ρ,z) (ρ is the distance for the rotation axis z). The effect of helicity fluctuations and the diffusion coefficient on the turbulent diffusion process is also investigated. Zh. éksp. Teor. Fiz. 116, 85–104 (July 1999)     

Evaporating droplets in turbulent reacting flows

Three-dimensional direct numerical simulations are carried out to determine the effects of turbulence on the preferential segregation of an evaporating spray and then to study the evolution of the resulting mixture fraction topology and propagating flame. First, the mixing between an initially randomly dispersed phase and the turbulent gaseous carrier phase is studied with non-evaporating particles. According to their inertia and the turbulence properties, the formation of clusters of particles is analyzed (formation delay, cluster characteristic size and density). Once the particles are in dynamical equilibrium with the surrounding turbulent flow, evaporation is considered through the analysis of the mixture fraction evolution. Finally, to mimic ignition, a kernel of burnt gases is generated at the center of the domain and the turbulent flame evolution is described.     

Effects of forcing in three-dimensional turbulent flows

We present the results of a numerical investigation of three-dimensional homogeneous and isotropic turbulence, stirred by a random forcing with a power-law spectrum, E(f)(k) approximately k(3-y). Numerical simulations are performed at different resolutions up to 512(3). We show that at varying the spectrum slope y, small-scale turbulent fluctuations change from a forcing independent to a forcing dominated statistics. We argue that the critical value separating the two behaviors, in three dimensions, is y(c)=4. When the statistics is forcing dominated, for yy(c), we find the same anomalous scaling measured in flows forced only at large scales. We connect these results with the issue of universality in turbulent flows.     

Statistics of Lagrangian velocities in turbulent flows

We present a generalized Fokker-Planck equation for the joint position-velocity probability distribution of a single fluid particle in a turbulent flow. Based on a simple estimate, the diffusion term is related to the two-point two-time Eulerian acceleration-acceleration correlation. Dimensional analysis yields a velocity increment probability distribution with normal scaling v approximately t(1/2). However, the statistics need not be Gaussian.     

Micromixing models for turbulent flows

In transported probability density function and filtered density function methods, micromixing models are required to close the molecular mixing term. The accuracy and computational efficiency of improved versions of the parameterized scalar profile (PSP) model are assessed and compared with commonly used mixing models such as Curl, modified Curl, interaction by exchange with the mean and Euclidean minimum spanning tree. Different generalizations of the PSP mixing model for spatially inhomogeneous flow configurations are presented. The selected test cases focus on molecular mixing and avoid interference with other models. Simulation results for a three-stream problem, involving two inert scalars, and a multi-scalar test case with mean-scalar-gradients are presented.     

Solutions for steady plane orthogonal MHD flows

Steady, plane flow of an inviscid, electrically conducting incompressible fluid with infinite electrical conductivity is considered and a single partial differential equation, in the case of orthogonal flow, is obtained which involves two functions. Appropriate specialization of these functions generates new exact solutions of the original equations.     

A fourth-order divergence-free method for MHD flows

This paper extends our previous third-order method [S. Li, High order central scheme on overlapping cells for magneto-hydrodynamic flows with and without constrained transport method, J. Comput. Phys. 227 (2008) 7368–7393] to the fourth-order. Central finite-volume schemes on overlapping grid are used for both the volume-averaged variables and the face-averaged magnetic field. The magnetic field at the cell boundaries falls within the dual grid and is naturally continuous so that our method eliminates the instability triggered by the discontinuity in the normal component of the magnetic field. Our fourth-order scheme has much smaller numerical dissipation than the third-order scheme. The divergence-free condition of the magnetic field is preserved by our fourth-order divergence-free reconstruction and the constrained transport method. Numerical examples show that the divergence-free condition is essential to the accuracy of the method when a limiter is used in the reconstruction. The high-order, low-dissipation, and divergence-free properties of this method make it an ideal tool for direct magneto-hydrodynamic turbulence simulations.     

Multiscale model of gradient evolution in turbulent flows

A multiscale model for the evolution of the velocity gradient tensor in turbulence is proposed. The model couples "restricted Euler" (RE) dynamics describing gradient self-stretching with a cascade model allowing energy exchange between scales. We show that inclusion of the cascade process is sufficient to regularize the finite-time singularity of the RE dynamics. Also, the model retains geometrical features of real turbulence such as preferential alignments of vorticity and joint statistics of gradient tensor invariants. Furthermore, gradient fluctuations are non-Gaussian, skewed in the longitudinal case, and derivative flatness coefficients are in good agreement with experimental data.     

Properties of steady states in turbulent axisymmetric flows

We experimentally study the properties of mean and most probable velocity fields in a turbulent von Kármán flow. These fields are found to be described by two families of functions, as predicted by a recent statistical mechanics study of 3D axisymmetric flows. We show that these functions depend on the viscosity and on the forcing. Furthermore, when the Reynolds number is increased, we exhibit a tendency for Beltramization of the flow, i.e., a velocity-vorticity alignment. This result provides a first experimental evidence of nonlinearity depletion in nonhomogeneous nonisotropic turbulent flow.     

Modelling of turbulent reacting flows in furnace devices

The spatial combustion of turbulent jets in furnace devices is modelled numerically basing on equations for multi-component turbulent reacting gaseous mixtures. The dependencies are obtained for the influence of the secondary air velocity and composition of gaseous components on torch configuration at a diffusion combustion process. The effect of regime parameters on the increase in torch sizes, which arises at the interaction of secondary air with gaseous components, has been elucidated.