Statistics of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence

We present a detailed direct numerical simulation of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work significantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields, and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distinguishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies.

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

Performance analysis of FSO links under strong atmospheric turbulence conditions using various modulation schemes

Strong atmospheric turbulence is a major hindrance in wireless optical communication systems. In this paper, the performance of a wireless optical communication system is analyzed using different modulation formats such as, binary phase shift keying-subcarrier intensity modulation (BPSK-SIM), differential phase shift keying (DPSK), differential phase shift keying-subcarrier intensity modulation (DPSK-SIM), Mary pulse position modulation (M-PPM) and polarization shift keying (PoISK). The atmospheric channel is modeled for strong atmospheric turbulences with combined effect of turbulence and pointing errors. Novel closed-form analytical expressions for average bit error rate (BER), channel capacity and outage probability for the various modulation techniques, viz. BPSK-SIM, DPSK, DPSK-SIM, PoISK and M-PPM are derived. The simulated results for BER, channel capacity and outage probability of various modulation techniques are plotted and analyzed. (C) 2014 Elsevier GmbH. All rights reserved.

Secondary instability as a possible mechanism for clear-air turbulence: a case study

Many theories and mechanisms have been proposed to explain the phenomenon of clear-air turbulence (CAT), and some of them have been successful in predicting light, moderate and, in some cases, severe turbulence. It is only recently that skill in the forecasting of the severe form of CAT, which could lead to injuries to passengers and damage to aircraft, has improved. Recent observations and simulations suggest that some severe to extreme turbulence could be caused by horizontal vortex tubes resulting from secondary instabilities of regions of high shear in the atmosphere. We have conducted direct numerical simulations to understand the scale relationship between primary structures (larger-scale structures related to one of the causes mentioned above) and secondary structures (smaller-sized, shear structures of the size of aircraft). From shear layer simulations, we find that the ratio of sizes of primary and secondary vortices is of the right order to generate aircraft-scale vortex tubes from typical atmospheric shear layers. We have also conducted simulations with a mesoscale atmospheric model, to understand possible causes of turbulence experienced by a flight off the west coast of India. Our simulations show the occurrence of primary flow structures related to synoptic conditions around the time of the incident. The evidence presented for this mechanism also has implications for possible methods of detection and avoidance of severe CAT.

Life of flame particles embedded in premixed flames interacting with near isotropic turbulence

Flame particles are surface points that always remain embedded on, by comoving with a given iso-scalar surface within a flame. Tracking flame particles allow us to study the fate of propagating surface locations uniquely identified throughout their evolution with time. In this work, using Direct Numerical Simulations we study the finite lifetime of such flame particles residing on iso-temperature surfaces of statistically planar H-2-air flames interacting with near-isotropic turbulence. We find that individual flame particles as well as their ensemble, experience progressively increasing tangential straining rate (K-t) and increasing negative curvature (kappa) near the end of their lifetime to finally get annihilated. By studying two different turbulent flow conditions, flame particle tracking shows that such tendency of local flame surfaces to be strained and cusped towards pinch-off from the main surface is a rather generic feature, independent of initial conditions, locations and ambient turbulence intensity levels. The evolution of the alignments between the flame surface normals and the principal components of the local straining rates are also tracked. We find that the surface normals initially aligned with the most extensive principal strain rate components, rotate near the end of flame particles' lifetime to enable preferential alignment between the surface tangent and the most extensive principal strain rate component. This could explain the persistently increasing tangential strain rate, sharp negative curvature formation and eventual detachment. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

On flame-turbulence interaction in constant-pressure expanding flames

In this paper we present one of the first high-speed particle image velocimetry measurements to quantify flame-turbulence interaction in centrally-ignited constant-pressure premixed flames expanding in nearisotropic turbulence. Measurements of mean flow velocity and rms of fluctuating flow velocity are provided over a range of conditions both in the presence and absence of the flame. The distributions of stretch rate contributions from different terms such as tangential straining, normal straining and curvature are also provided. It is found that the normal straining displays non-Gaussian pdf tails whereas the tangential straining shows near Gaussian behavior. We have further tracked the motion of the edge points that reside and co-move with the edge of the flame kernel during its evolution in time, and found that within the measurement conditions, on average the persistence time scales of stretch due to pure curvature exceed that due to tangential straining by at least a factor of two. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Low frequency modulation of jets in quasigeostrophic turbulence

Quasigeostrophic turbulence on a beta-plane with a finite deformation radius is studied numerically, with particular emphasis on frequency and combined wavenumber-frequency domain analyses. Under suitable conditions, simulations with small-scale random forcing and large-scale drag exhibit a spontaneous formation of multiple zonal jets. The first hint of wave-like features is seen in the distribution of kinetic energy as a function of frequency; specifically, for progressively larger deformation scales, there are systematic departures in the form of isolated peaks (at progressively higher frequencies) from a power-law scaling. Concomitantly, there is an inverse flux of kinetic energy in frequency space which extends to lower frequencies for smaller deformation scales. The identification of these peaks as Rossby waves is made possible by examining the energy spectrum in frequency-zonal wavenumber and frequency-meridional wavenumber diagrams. In fact, the modified Rhines scale turns out to be a useful measure of the dominant meridional wavenumber of the modulating Rossby waves; once this is fixed, apart from a spectral peak at the origin (the steady jet), almost all the energy is contained in westward propagating disturbances that follow the theoretical Rossby dispersion relation. Quite consistently, noting that the zonal scale of the modulating waves is restricted to the first few wavenumbers, the energy spectrum is almost entirely contained within the corresponding Rossby dispersion curves on a frequency-meridional wavenumber diagram. Cases when jets do not form are also considered; once again, there is a hint of Rossby wave activity, though the spectral peaks are quite muted. Further, the kinetic energy scaling in frequency domain follows a -5/3 power-law and is distributed much more broadly in frequency-wavenumber diagrams. (C) 2015 AIP Publishing LLC.

Two-dimensional homogeneous isotropic fluid turbulence with polymer additives

We carry out an extensive and high-resolution direct numerical simulation of homogeneous, isotropic turbulence in two-dimensional fluid films with air-drag-induced friction and with polymer additives. Our study reveals that the polymers (a) reduce the total fluid energy, enstrophy, and palinstrophy; (b) modify the fluid energy spectrum in both inverse-and forward-cascade regimes; (c) reduce small-scale intermittency; (d) suppress regions of high vorticity and strain rate; and (e) stretch in strain-dominated regions. We compare our results with earlier experimental studies and propose new experiments.

Reynolds averaged Navier-Stokes simulations of compressible mixing layers of similar and dissimilar gases: Performance of k- turbulence model

The issue of growth rate reduction of high speed mixing layer with convective Mach number is examined for similar and dissimilar gases using Reynolds averaged Navier-Stokes (RANS) methodology with k- turbulence model. It is observed that the growth rate predicted using RANS simulations closely matches with that predicted using model free simulations. Velocity profiles do not depend on the modelled value of Pr-t and Sc-t; while the temperature and species mass fraction distributions depend heavily on them. Although basic k- turbulence model could not capture the reduced growth rate for the mixing layer formed between similar gases, it predicts very well the reduced growth rate for the mixing layer for the dissimilar gases. It appears that density ratio changes caused by temperature changes for the dissimilar gases have profound effect on the growth rate reduction.

Effect of a mesh on boundary layer transitions induced by free-stream turbulence and an isolated roughness element

Streamwise streaks, their lift-up and streak instability are integral to the bypass transition process. An experimental study has been carried out to find the effect of a mesh placed normal to the flow and at different wall-normal locations in the late stages of two transitional flows induced by free-stream turbulence (FST) and an isolated roughness element. The mesh causes an approximately 30% reduction in the free-stream velocity, and mild acceleration, irrespective of its wall-normal location. Interestingly, when located near the wall, the mesh suppresses several transitional events leading to transition delay over a large downstream distance. The transition delay is found to be mainly caused by suppression of the lift-up of the high-shear layer and its distortion, along with modification of the spanwise streaky structure to an orderly one. However, with the mesh well away from the wall, the lifted-up shear layer remains largely unaffected, and the downstream boundary layer velocity profile develops an overshoot which is found to follow a plane mixing layer type profile up to the free stream. Reynolds stresses, and the size and strength of vortices increase in this mixing layer region. This high-intensity disturbance can possibly enhance transition of the accelerated flow far downstream, although a reduction in streamwise turbulence intensity occurs over a short distance downstream of the mesh. However, the shape of the large-scale streamwise structure in the wall-normal plane is found to be more or less the same as that without the mesh.

Origin of nonlinearity and plausible turbulence by hydromagnetic transient growth in accretion disks: Faster growth rate than magnetorotational instability

We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that a mechanism, known as magnetorotational instability (MRI), is responsible for transporting matter in the presence of a weak magnetic field. However, there are some shortcomings, which question the effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g., transient growth (TG), can play an important role in bringing nonlinearity into the system, even at weak magnetic fields. In addition, it should be determined whether MRI or TG is primarily responsible for revealing nonlinearity in order to make the flow turbulent. Our results prove explicitly that the flows with a high Reynolds number (Re), which is the case for realistic astrophysical accretion disks, exhibit nonlinearity via TG of perturbation modes faster than that by modes producing MRI. For a fixed wave vector, MRI dominates over transient effects only at low Re, lower than the value expected to be in astrophysical accretion disks, and low magnetic fields. This calls into serious question the (overall) persuasiveness of MRI in astrophysical accretion disks.

TURBULENCE-TRANSPORT-CHEMISTRY INTERACTION IN STATISTICALLY PLANAR PREMIXED FLAMES AND IGNITION KERNELS IN NEAR ISOTROPIC TURBULENCE

Turbulence-transport-chemistry interaction plays a crucial role on the flame surface geometry, local and global reactionrates, and therefore, on the propagation and extinction characteristics of intensely turbulent, premixed flames encountered in LPP gas-turbine combustors. The aim of the present work is to understand these interaction effects on the flame surface annihilation and extinction of lean premixed flames, interacting with near isotropic turbulence. As an example case, lean premixed H-2-air mixture is considered so as to enable inclusion of detailed chemistry effects in Direct Numerical Simulations (DNS). The work is carried out in two phases namely, statistically planar flames and ignition kernel, both interacting with near isotropic turbulence, using the recently proposed Flame Particle Tracking (FPT) technique. Flame particles are surface points residing and commoving with an iso-scalar surface within a premixed flame. Tracking flame particles allows us to study the evolution of propagating surface locations uniquely identified with time. In this work, using DNS and FPT we study the flame speed, reaction rate and transport histories of such flame particles residing on iso-scalar surfaces. An analytical expression for the local displacement flame speed (SO is derived, and the contribution of transport and chemistry on the displacement flame speed is identified. An examination of the results of the planar case leads to a conclusion that the cause of variation in S-d may be attributed to the effects of turbulent transport and heat release rate. In the second phase of this work, the sustenance of an ignition kernel is examined in light of the S-curve. A newly proposed Damkohler number accounting for local turbulent transport and reaction rates is found to explain either the sustenance or otherwise propagation of flame kernels in near isotropic turbulence.

Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.

A flame particle tracking analysis of turbulence-chemistry interaction in hydrogen-air premixed flames

Interactions of turbulence, molecular transport, and energy transport, coupled with chemistry play a crucial role in the evolution of flame surface geometry, propagation, annihilation, and local extinction/re-ignition characteristics of intensely turbulent premixed flames. This study seeks to understand how these interactions affect flame surface annihilation of lean hydrogen-air premixed turbulent flames. Direct numerical simulations (DNSs) are conducted at different parametric conditions with a detailed reaction mechanism and transport properties for hydrogen-air flames. Flame particle tracking (FPT) technique is used to follow specific flame surface segments. An analytical expression for the local displacement flame speed (S-d) of a temperature isosurface is considered, and the contributions of transport, chemistry, and kinematics on the displacement flame speed at different turbulence-flame interaction conditions are identified. In general, the displacement flame speed for the flame particles is found to increase with time for all conditions considered. This is because, eventually all flame surfaces and their resident flame particles approach annihilation by reactant island formation at the end of stretching and folding processes induced by turbulence. Statistics of principal curvature evolving in time, obtained using FPT, suggest that these islands are ellipsoidal on average enclosing fresh reactants. Further examinations show that the increase in S-d is caused by the increased negative curvature of the flame surface and eventual homogenization of temperature gradients as these reactant islands shrink due to flame propagation and turbulent mixing. Finally, the evolution of the normalized, averaged, displacement flame speed vs. stretch Karlovitz number are found to collapse on a narrow band, suggesting that a unified description of flame speed dependence on stretch rate may be possible in the Lagrangian description. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.