Prandtl number effects in decaying homogeneous isotropic turbulence with a mean scalar gradient

Decaying homogeneous isotropic turbulence with an imposed mean scalar gradient is investigated numerically, thanks to a specific eddy-damped quasi-normal Markovian closure developed recently for passive scalar mixing in homogeneous anisotropic turbulence (BGC). The present modelling is compared successfully with recent direct numerical simulations and other models, for both very large and small Prandtl numbers. First, scalings for the cospectrum and scalar variance spectrum in the inertial range are recovered analytically and numerically. Then, at large Reynolds numbers, the decay and growth laws for the scalar variance and mixed velocity–scalar correlations, respectively, derived in BGC, are shown numerically to remain valid when the Prandtl number strongly departs from unity. Afterwards, the normalised correlation ρ_wθ is found to decrease in magnitude at a fixed Reynolds number when Pr either increases or decreases, in agreement with earlier predictions. Finally, the small scales return to isotropy of the scalar second-order moments is found to depend not only on the Reynolds number, but also on the Prandtl number.     

Reassessment of the classical closures for scalar turbulence

In deducing the consequences of the Direct Interaction Approximation, Kraichnan was sometimes led to consider the properties of special classes of nonlinear interactions in degenerate triads in which one wavevector is very small. Such interactions can be described by simplified models closely related to elementary closures for homogeneous isotropic turbulence such as the Heisenberg and Leith models. These connections can be exploited to derive considerably improved versions of the Heisenberg and Leith models that are only slightly more complicated analytically. This paper applies this approach to derive some new simplified closure models for passive scalar advection and investigates the consistency of these models with fundamental properties of scalar turbulence. Whereas some properties, such as the existence of the Kolmogorov–Obukhov range and the existence of thermal equilibrium ensembles, follow the velocity case closely, phenomena special to the scalar case arise when the diffusive and viscous effects become important at different scales of motion. These include the Batchelor and Batchelor–Howells–Townsend ranges pertaining, respectively, to high and low molecular Schmidt number. We also consider the spectrum in the diffusive range that follows the Batchelor range. We conclude that improved elementary models can be made consistent with many nontrivial properties of scalar turbulence, but that such models have unavoidable limitations.     

Closure Approximations for Passive Scalar Turbulence: A Comparative Study on an Exactly Solvable Model with Complex Features

Some standard closure approximations used in turbulence theory are analyzed by examining systematically the predictions these approximations produce for a passive scalar advection model consisting of a shear flow with a fluctuating cross sweep. This model has a general geometric structure of a jet flow with transverse disturbances, which occur in a number of contexts, and it encompasses a wide variety of possible spatio-temporal statistical structures for the velocity field, including strong long-range correlations. Even though the Eulerian and Lagrangian velocity statistics are not equal and the passive scalar statistics exhibit broader-than-Gaussian intermittency, this model is nevertheless simple enough so that many passive scalar statistics can be computed exactly and compared systematically with the predictions of the closure approximations. Our comparative study illustrates the strength and weaknesses of the closure approximations and points out the physical phenomena that these approximations are able or not able to describe properly. In particular it is shown that the direct interaction approximation (DIA), one of the most sophisticated closure approximations available, fails to reproduce adequately the statistical features of the scalar and may even lead to absdurd predictions, even though the equations it produces are rather complicated and difficult to analyze. Two alternative closure approximations, the Modified DIA (MDIA) and the Renormalized Lagrangian Approximation (RLA), with different levels of sophistication, both are simpler to use than the DIA and perform better. In particular, it is shown that both closure approximations always reproduce exactly the second order statistics for the scalar and that the MDIA is even able to capture intermittency effects.     

Turbulent scalar fluxes in H2-air premixed flames at low and high Karlovitz numbers

The statistical behaviour and closure of sub-grid scalar fluxes in the context of turbulent premixed combustion have been assessed based on an a priori analysis of a detailed chemistry Direct Numerical Simulation (DNS) database consisting of three hydrogen-air flames spanning the corrugated flamelets (CF), thin reaction zones (TRZ) and broken reaction zones (BRZ) regimes of premixed turbulent combustion. The sub-grid scalar fluxes have been extracted by explicit filtering of DNS data. It has been found that the conventional gradient hypothesis model is not appropriate for the closure of sub-grid scalar flux for any scalar in the context of a multispecies system. However, the predictions of the conventional gradient hypothesis exhibit a greater level of qualitative agreement with DNS data for the flame representing the BRZ regime. The aforementioned behaviour has been analysed in terms of the properties of the invariants of the anisotropy tensor in the Lumley triangle. The flames in the CF and TRZ regimes are characterised by a pronounced two-dimensional anisotropy due to strong heat release whereas a three-dimensional and more isotropic behaviour is observed for the flame located in the BRZ regime. Two sub-grid scalar flux models which are capable of predicting counter-gradient transport have been considered for a priori DNS assessment of multispecies systems and have been analysed in terms of both qualitative and quantitative agreements. By combining the latter two sub-grid scalar flux closures, a new modelling strategy is suggested where one model is responsible for properly predicting the conditional mean accurately and the other model is responsible for the correlations between model and unclosed term. Detailed physical explanations for the observed behaviour and an assessment of existing modelling assumptions have been provided. Finally, the classical Bray–Moss–Libby theory for the scalar flux closure has been extended to address multispecies transport in the context of large eddy simulations.     

标量湍流的能量传输特性

应用直接数值模拟数据,从标量湍流传输的三波关系出发,进行湍流及标量湍流传输谱的多尺度分析,研究不同尺度间的能量传输性质,证实标量能量的传输与湍动能传输具有不同性质,大尺度速度脉动对标量传输有较大贡献,尤其是与标量小尺度脉动的相互作用,使标量模拟需要有比速度场更高的网格分辨率;并发现标量湍流的能量传输具有明显的非局部性;另外,定义了能量传输系数,发现在相同的Re数和Pe数条件下,标量湍流的对流惯性较速度脉动的惯性子区宽.     

Fundamental weak interaction studies using polarised nuclei and ion traps

Two experiments to search for new physics beyond the standard model for electroweak interactions by measuring correlations between different spin and momentum vectors in nuclear β-decay are discussed. In the first experiment the correlation between the emission asymmetry and the longitudinal polarisation of positrons emitted by polarised nuclei is determined. This type of measurement is sensitive to the presence of right-handed currents but also to possible scalar and tensor-type currents in the weak interaction. The aim of the second experiment is to determine the βν-correlation in β-decay by measuring the energy spectrum of the recoil ions, using a Penning trap and a retardation spectrometer. In this case the focus is on the search for scalar currents in the weak interaction. The results of the experiments presented here are complementary to results from experiments in muon decay and at high-energy colliders. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modelling of turbulent scalar flux in turbulent premixed flames based on DNS databases

A transport equation for scalar flux in turbulent premixed flames was modelled on the basis of DNS databases. Fully developed turbulent premixed flames were obtained for three different density ratios of flames with a single-step irreversible reaction, while the turbulent intensity was comparable to the laminar burning velocity. These DNS databases showed that the countergradient diffusion was dominant in the flame region. Analyses of the Favre-averaged transport equation for turbulent scalar flux proved that the pressure related terms and the velocity–reaction rate correlation term played important roles on the countergradient diffusion, while the mean velocity gradient term, the mean progress variable gradient term and dissipation terms suppressed it. Based on these analyses, modelling of the combustion-related terms was discussed. The mean pressure gradient term and the fluctuating pressure term were modelled by scaling, and these models were in good agreement with DNS databases. The dissipation terms and the velocity–reaction rate correlation term were also modelled, and these models mimicked DNS well.     

Gradient trajectory analysis in a Jet flow for turbulent combustion modelling

Based on planar high-speed Rayleigh scattering measurements of the mixture fraction Z of propane discharging from a turbulent round jet into co-flowing carbon dioxide at nozzle-based Reynolds numbers Re ₀ = 3000–8600, we use scalar gradient trajectories to investigate the local structure of the turbulent scalar field with a focus on the scalar turbulent/non-turbulent interface. The latter is located between the fully turbulent part of the jet and the outer flow. Using scalar gradient trajectories, we partition the turbulent scalar field into these three regions according to an approach developed by Mellado et al. (J.P. Mellado, L. Wang, and N. Peters, Gradient trajectory analysis of a scalar field with external intermittency, J. Fluid Mech. 626 (2009), pp. 333–365.). Based on these different regions, we investigate in a next step zonal statistics of the scalar probability density function (pdf) P(Z) as well as the scalar difference along the trajectory ΔZ and its mean scalar value Z_m , where the latter two quantities are used to parameterize the scalar profile along gradient trajectories. We show that the scalar pdf P(Z) can be reconstructed from zonal gradient trajectory statistics of the joint pdf P(Z_m , ΔZ). Furthermore, on the one hand we relate our results for the scalar turbulent/non-turbulent interface to the findings made in other experimental and numerical studies of the turbulent/non-turbulent interface, and on the other hand discuss them in the context of the flamelet approach and the modelling of pdfs in turbulent non-premixed combustion. Finally, we compare the zonal statistics for P(Z) with the composite model of Effelsberg and Peters (E. Effelsberg and N. Peters, A composite model for the conserved scalar pdf, Combust. Flame 50 (1983), pp. 351–360) and observe a very good qualitative and quantitative agreement.     

Conditional velocity statistics in the double scalar mixing layer – A mapping closure approach

In this work we use 3D direct numerical simulations (DNS) to investigate the average velocity conditioned on a conserved scalar in a double scalar mixing layer (DSML). The DSML is a canonical multistream flow designed as a model problem for the extensively studied piloted diffusion flames. The conditional mean velocity appears as an unclosed term in advanced Eulerian models of turbulent non-premixed combustion, like the conditional moment closure and transported probability density function (PDF) methods. Here it accounts for inhomogeneous effects that have been found significant in flames with relatively low Damköhler numbers. Today there are only a few simple models available for the conditional mean velocity and these are discussed with reference to the DNS results. We find that both the linear model of Kutznetzov and the Li and Bilger model are unsuitable for multi stream flows, whereas the gradient diffusion model of Pope shows very close agreement with DNS over the whole range of the DSML. The gradient diffusion model relies on a model for the conserved scalar PDF and here we have used a presumed mapping function PDF, that is known to give an excellent representation of the DNS. A new model for the conditional mean velocity is suggested by arguing that the Gaussian reference field represents the velocity field, a statement that is evidenced by a near perfect agreement with DNS. The model still suffers from an inconsistency with the unconditional flux of conserved scalar variance, though, and a strategy for developing fully consistent models is suggested.     

Conserved scalar mixing in a confined-opposed-jets flow

The focus of this paper is on the mixing of a conserved passive scalar for Sc = 1 (Sc is the Schmidt number) in axisymmetric turbulence for which the initial injections of turbulent kinetic energy and scalar variance are similar. Two confined-opposed-jets (COJ) are experimentally studied through simultaneous PIV (particle image velocimetry) and PLIF (planar laser induced fluorescence) measurements, for different flow regimes. One-point transport equation for the scalar variance is assessed through experimental data, along the common axis of the two opposed jets, and different physical phenomena are revealed (production, diffusion, dissipation). The production of scalar variance is equilibrated by the diffusion term (～75%) and the mean dissipation of the scalar variance (～25%). To further assess the scalar behaviour at each scale in this anisotropic, but axisymmetric, flow, a scale-by-scale scalar variance budget equation is derived for axisymmetric turbulence. This equation reduces to Yaglom's 4/3 law, under additional restrictions. The equation is assessed through experimental data, in the impingement region between the two COJ. In particular, the anisotropic energy transfer along different directions is quantified. It is shown that for scales smaller than the size of the central region, Δ, the cascade of the scalar variance is completely inhibited, independently of the particular direction. For scales larger than Δ, the apparent aspect of the energy transfer is that of an inverse cascade, with positive values of the scalar variance transfer. Nonetheless, inhomogeneity of the flow and mixing at those scales is directly responsible for these positive values.     

Mixing of a passive scalar across a thin shearless layer: concentration of intermittency on the sides of the turbulent interface

The advection of a passive scalar through an initial flat interface separating two different isotropic decaying turbulent fields is investigated in two and three dimensions. Simulations have been performed for a range of Taylor’s microscale Reynolds numbers from 45 to 250 and for a Schmidt number equal to 1. Different to the case where the transport involves the momentum and kinetic energy only and one intermittency layer is formed in the low-turbulent energy side of the system, in the passive scalar concentration field two intermittent layers are observed to develop at the sides of the interface. The layers move normally to the interface in opposite directions. The dimensionality produces different time scaling of the passive scalar diffusion, which is much faster in the two-dimensional case. In two dimensions, the propagation of the intermittent layers exhibits a significant asymmetry with respect to the initial position of the interface and is deeper for the layer which moves towards the high kinetic energy side of the system. In three dimensions, the two intermittent layers propagate nearly symmetrically with respect the centre of the mixing region. During the temporal decay, inside the mixing, which is both inhomogeneous and anisotropic but devoid of a mean velocity shear, the passive scalar spectra are computed. In three dimensions, the exponent in the scaling range gets in time a value close to that of the kinetic energy spectrum of isotropic turbulence (?5/3). In two dimensions, instead the exponent settles down to a value that is about one-half of the corresponding isotropic case. By means of an analysis based on simple wavy perturbations of the interface we show that the formation of the double layer of intermittency is a dynamic general feature not specific to the turbulent transport. These results of our numerical study are discussed in the context of experimental results and numerical simulations.     

Influence of velocity spatiotemporal correlations on the anomalous scaling exponents of passive scalars

In this paper, we consider spatial-temporal correlation functions of the turbulent velocities. With numerical simulations on the Gledzer--Ohkitani--Yamada (GOY) shell model, we show that the correlation function decays exponentially. The advecting velocity field is regarded as a colored noise field, which is spatially and temporally correlative. For comparison, we are also given the scaling exponents of passive scalars obtained by the Gaussian random velocity field, the multi-dimensional normal velocity field and the She--Leveque velocity field, introduced by She, et al. We observe that extended self-similarity scaling exponents H(p)/ H(2) of passive scalar obtained by the colored noise field are more anomalous than those obtained by the other three velocity fields.     

EXACT SOLUTIONS OF CLASSICAL SCALAR FIELD EQUATIONS

We give a class of exact solutions of quartic scalar field theories. These solutions prove to be interesting as are characterized by the production of mass contributions arising from the nonlinear terms while maintaining a wave-like behavior. So, a quartic massless equation has a nonlinear wave solution with a dispersion relation of a massive wave and a quartic scalar theory gets its mass term renormalized in the dispersion relation through a term depending on the coupling and an integration constant. When spontaneous breaking of symmetry is considered, such wave-like solutions show how a mass term with the wrong sign and the nonlinearity give rise to a proper dispersion relation. These latter solutions do not change the sign maintaining the property of the selected value of the equilibrium state. Then, we use these solutions to obtain a quantum field theory for the case of a quartic massless field. We get the propagator from a first-order correction showing that is consistent in the limit of a very large coupling. The spectrum of a massless quartic scalar field theory is then provided. From this we can conclude that, for an infinite countable number of exact classical solutions, there exist an infinite number of equivalent quantum field theories that are trivial in the limit of the coupling going to infinity.     

Measurement of small spin-spin splittings in the presence of chemical exchange: case of deuteriated water

A new implementation of a described data processing for measuring small scalar coupling constant J is proposed. It is applied particularly to situations where the coupling constant is around the linewidth at half-height, and in the presence of chemical exchange. This modified data processing uses only a simple "one-pulse" experiment instead of a series of spin echoes experiments required by the previous processing. The FID recorded in the one pulse experiment is used by a reconstitution program to generate a set of signals, which are analyzed in the time domain to obtain a spectrum where the scalar coupling constant is apparently multiplied by n+1, where n is positive. The new processing is tested with simulated spectra. The coupling constant between proton and deuterium is measured in the proton spectrum of a solution of 80% of D(2)O and 20% of H(2)O. It was found to be J=1.54+/-0.01 Hz.     

Influence of Velocity Shell Correlations on Anomalous Scaling Exponents of Passive Scalars in Shell Models

We propose a new approach to the old-standing problem of the anomaly of the scaling exponents of passive scalars of turbulence. Different to the original problem, the distribution function of the prescribed random velocity field is multi-dimensional normal and delta-correlated in time. Here, our random velocity field is spatially correlative. For comparison, we also give the result obtained by the Gaussian random velocity field without spatial correlation. The anomalous scaling exponents H(p) of passive scalar advected by two kinds of random velocity above are determined for structure function up to p=15 by numerical simulations of the random shell model with Runge-Kutta methods to solve the stochastic differential equations. We observed that the H(p)'s obtained by the multi-dimensional normal distribution random velocity are more anomalous than those obtained by the independent Gaussian random velocity.     

One self-consistent closure for chains of equations for Green's functions. An anisotropic heisenberg model

This article present a self-consistent approach to computation of the correlation function in the method of Green's functions. The basis for the approach is representation of the desired Green's function in the form of a chain of fractions that is subsequently closed. The closure is based on the use of concrete relations imposed on the higher-order correlation functions. General expressions for the correlation functions and conditions for self-consistency of the computations are presented. The method has been tested by computing the magnetization and critical temperature in a Heisenberg model with arbitrary anisotropy parameter. We obtain general expressions for these quantities. The critical temperature obtained is less than the corresponding value given by the molecular-field approximation. The latter approximation also leads to an overestimate of the magnetization values. It is shown that no critical transition is possible for any value of the anisotropy parameter. The corresponding inequalities are obtained. The method is compared with a method for self-consistent computation of correlation functions that was proposed earlier by the author.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii Fizika, No. 1, pp. 46–52, January, 1996.     

Influence of advective heat flux on extinction scalar dissipation rate and velocity in negative edge flames

A counterflow flame geometry, which has previously been experimentally shown to produce stable negative edge flames, was studied using numerical simulations. In this geometry, the flame edge is formed off the counterflow centreline owing to a local increase in scalar dissipation rate. Hot products from the stable nonpremixed flame on the centreline flow through the edge at velocities of ～ 1–5 m/s. The size of the counterflow burner and the gas flowrates are varied in the simulations to alter the flame strength and velocity at the flame edge. The advection of products through the edge is shown to extend the flame extinction to higher scalar dissipation rates than required for centreline extinction. For high velocities, the scalar dissipation rate required for flame extinction can be related to the centreline extinction value by considering only the effect of energy addition to the flame edge via advection. However, for lower edge flame velocities, the effects of increased thermal and species diffusion through the edge must also be included. Since the advection at the edge is a product of both the local velocity and temperature gradient, a single correlation between the scalar dissipation rate and the negative edge flame velocity does not exist.     

Slow Modes in Passive Advection

The anomalous scaling in the Kraichnan model of advection of the passive scalar by a random velocity field with nonsmooth spatial behavior is traced to the presence of slow resonance-type collective modes of the stochastic evolution of fluid trajectories. We show that the slow modes are organized into infinite multiplets of descendants of the primary conserved modes. Their presence is linked to the nondeterministic behavior of the Lagrangian trajectories at high Reynolds numbers caused by the sensitive dependence on initial conditions within the viscous range where the velocity fields are more regular. Revisiting the Kraichnan model with smooth velocities, we describe the explicit solution for the stationary state of the scalar. The properties of the probability distribution function of the smeared scalar in this state are related to a quantum mechanical problem involving the Calogero–Sutherland Hamiltonian with a potential.Partially supported by NSF Grant DMS-9205296     

The double Lagrangian its applications with U（1, J） symmetry and to the cosmology

By introducing the double complex scalar field, this paper constructs the double Lagrangian with the potential V（Ф（J）,Ф（J）） = V（│Ф（J）│）, which not only can describe the evolution of quintom Universe, but also can naturally give the spintessence and hessence Universe. Furthermore, the U（1, J） symmetry of the double complex Lagrangian is verified, and the total conserved charge within the physical volume is derived by means of the Norther theorem. Moreover, it can point out that the ＇imaginary motion＇ of the angular parameter in Ref.[14], in fact, is only a real phase displacement in our model.