Transition from Countergradient to Gradient Scalar Transport in Developing Premixed Turbulent Flames

A simple model of turbulent scalar flux developed recently by the present authors is applied to determine the direction of the flux in a statistically planar one-dimensional premixed flame that does not affect turbulence and has self-similar mean structure. Results obtained in the case of statistically stationary turbulence indicate that transition from countergradient to gradient turbulent scalar transport may occur during flame development, as the peak mean rate of product creation moves to the trailing edge of the flame brush. In the case of decaying turbulence, the opposite transition (from gradient to countergradient transport) was simulated in line with available DNS data. In both cases, transition instant depends strongly on turbulence and mixture characteristics. In particular, countergradient transport is suppressed by an increase in the rms turbulent velocity and by a decrease in the laminar flame speed or density ratio, in line with available experimental and DNS data. The obtained results lend qualitative support to the model of turbulent scalar flux addressed in the present work.     

Structure-based transient models for scalar dissipation rate in homogeneous turbulence

We present, for the first time, transient models for scalar-dissipation rate that carry information about the morphology of the turbulence structures. The proposed models are derived from a structure-sensitized transient model for the large-scale field of enstrophy. Special attention is given to the model equation of Yoshizawa (1988), from which alternative closures for scalar-dissipation rate are derived. Assuming negligible transient variations of the turbulent anisotropies, lead to an expression similar to the one used by traditional transient models for scalar dissipation rate. To reduce model uncertainties, the turbulent scalar-flux vector and the turbulence structure tensors are obtained from the same closure models for all cases considered. The estimation performance of the models is evaluated on eight cases of deformation in both fixed and rotating frames, showing encouraging results. The proposed models achieve fair agreement with homogeneous Direct Numerical Simulations (DNS) predictions in the presence of either transverse or streamwise mean scalar gradient.     

Statistics and Modelling of Turbulent Kinetic Energy Transport in Different Regimes of Premixed Combustion

The statistical behaviour of turbulent kinetic energy transport in turbulent premixed flames is analysed using data from three-dimensional Direct Numerical Simulation (DNS) of freely propagating turbulent premixed flames under decaying turbulence. For flames within the corrugated flamelets regime, it is observed that turbulent kinetic energy is generated within the flame brush. By contrast, for flames within the thin reaction zones regime it has been found that the turbulent kinetic energy decays monotonically through the flame brush. Similar trends are observed also for the dissipation rate of turbulent kinetic energy. Within the corrugated flamelets regime, it is demonstrated that the effects of the mean pressure gradient and pressure dilatation within the flame are sufficient to overcome the effects of viscous dissipation and are responsible for the observed augmentation of turbulent kinetic energy in the flame brush. In the thin reaction zones regime, the effects of the mean pressure gradient and pressure dilatation terms are relatively much weaker than those of viscous dissipation, resulting in a monotonic decay of turbulent kinetic energy across the flame brush. The modelling of the various unclosed terms of the turbulent kinetic energy transport equation has been analysed in detail. The predictions of existing models are compared with corresponding quantities extracted from DNS data. Based on this a-priori DNS assessment, either appropriate models are identified or new models are proposed where necessary. It is shown that the turbulent flux of turbulent kinetic energy exhibits counter-gradient (gradient) transport wherever the turbulent scalar flux is counter-gradient (gradient) in nature. A new model has been proposed for the turbulent flux of turbulent kinetic energy, and is found to capture the qualitative and quantitative behaviour obtained from DNS data for both the corrugated flamelets and thin reaction zones regimes without the need to adjust any of the model constants.     

Accounting for Buoyancy Effects in the Explicit Algebraic Stress Model: Homogeneous Turbulent Shear Flows

Most explicit algebraic stress models are formulated for turbulent shear flows without accounting for external body force effects, such as the buoyant force. These models yield fairly good predictions of the turbulence field generated by mean shear. As for thermal turbulence generated by the buoyant force, the models fail to give satisfactory results. The reason is that the models do not explicitly account for buoyancy effects, which interact with the mean shear to enhance or suppress turbulent mixing. Since applicable, coupled differential equations have been developed describing these thermal turbulent fields, it is possible to develop corresponding explicit algebraic stress models using tensor representation theory. While the procedure to be followed has been employed previously, unique challenges arise in extending the procedure for developing the algebraic representations to turbulent buoyant flows. In this paper the development of an explicit algebraic stress model (EASM) is confined to the homogeneous buoyant shear flow case to illustrate the methodology needed to develop the proper polynomial representations. The derivation is based on the implicit formulation of the Reynolds stress anisotropy at buoyant equilibrium. A five-term representation is found to be necessary to account properly for the effect of the thermal field. Thus derived, external buoyancy effects are represented in the scalar coefficients of the basis tensors, and structural buoyancy effects are accounted for in additional terms in the stress anisotropy tensor. These terms will not vanish even in the absence of mean shear. The performance of the new EASM, together with a two-equation (2-Eq) model, the non-buoyant EASM of Gatski and Speziale (1993) and a full second-order model, is assessed against direct numerical simulations of homogeneous, buoyant shear flows at two different Richardson numbers representing weak and strong buoyancy effects. The calculations show that this five-term representation yields better results than the 2-Eq model and the EASM of Gatski and Speziale where buoyancy effects are not explicitly accounted for. Received 5 March 2001 and accepted 15 January 2002     

The decay of weak,homogeneous turbulence in the presence of a vertical body force and concentration gradient with a first-order reaction

The effects of buoyancy, produced by a uniform vertical concentration gradient and body force, on a homogeneous turbulent field accompanied by a first-order chemical reaction, are analysed by considering a simplified model. A system of two-point correlation equations, which contains mean concentration gradient and body force terms, is constructed from the Navier-Stokes, convective diffusion and continuity equations. By well-known methods, these equations are converted into equations for the spectrum functions in the wave-number space and solutions for different spectral tensors are obtained by neglecting the contributions of the triple correlation terms. For carrying out the numerical calculations, it is assumed that the turbulence is initially isotropic and the concentration fluctuations initially zero. It turns out that the turbulence decays with time, although the buoyancy forces do alter the rate of decay. The buoyancy forces can either extract energy from the turbulent field or feed energy into it, depending upon the direction of the body force and the concentration gradient. Spectra are displayed graphically for several values of the reaction rate parameter for stabilizing, as well as destabilizing, buoyancy forces.     

Mixing in a stably stratified medium by horizontal shear near vertical walls

Stratified environmental flows near boundaries can have a horizontal mean shear component, orthogonal to the vertical mean density gradient. Vertical transport, against the stabilizing force of gravity, is possible in such situations if three-dimensional turbulence is sustained by the mean shear. A model problem, water with a constant mean density gradient flowing in a channel between parallel vertical walls, is examined here using the technique of large eddy simulation (LES). It is found that, although the mean shear is horizontal, the fluctuating velocity field has significant vertical shear and horizontal vorticity, thereby causing small-scale vertical mixing of the density field. The vertical stirring is especially effective near the boundaries where the mean shear is large and, consequently, the gradient Richardson number is small. The mean stratification is systematically increased between cases in our study and, as expected, the buoyancy flux correspondingly decreases. Even so, horizontal mean shear is found to be more effective than the well-studied case of mean vertical shear in inducing vertical buoyancy transport as indicated by generally larger values of vertical eddy diffusivity and mixing efficiency.     

The Scalar Gradient Alignment Statistics of Flame Kernels and its Modelling Implications for Turbulent Premixed Combustion

The effects of mean flame curvature on reaction progress variable gradient, $\nabla c$ , alignment with local turbulent strain rate are studied based on three-dimensional Direct Numerical Simulation (DNS) data of turbulent premixed flame kernels with different initial radii under decaying turbulence. A statistically planar flame is also considered in order to compare the results obtained from the kernels with a flame of zero mean curvature. It is found that the dilatation rate effects diminish with decreasing kernel radius due to defocusing of heat in the positively curved regions. This gives rise to a decrease in the extent of reaction progress variable gradient alignment with most extensive principal strain rate with decreasing kernel radius. The modelling implications of the statistics of the alignment of $\nabla c$ with local strain rate have been studied in terms of scalar dissipation rate transport. A new modelling methodology for the contribution of the scalar-turbulence interaction term in the transport equation for the mean scalar dissipation is suggested addressing the reduced effects of dilatation rate for flame kernels and the diminished value of turbulent straining at the small length scales at which turbulence interacts with small flame kernels. The performance of the new models is found to be satisfactory while comparing to DNS results. The existing models for the dilatation contribution and the combined chemical reaction and molecular dissipation contributions to the transport of mean scalar dissipation, which were originally proposed for statistically planar flames, are found to satisfactorily predict the corresponding quantities for turbulent flame kernels.     

Numerical Study and Lagrangian Modelling of Turbulent Heat Transport

A stochastic Lagrangian model for both fluid velocities and temperature fluctuations is evaluated from Direct Numerical Simulation of heat transport in homogeneous isotropic turbulence submitted to a linear mean temperature gradient. The first stage lies on the study of the Lagrangian fluid turbulence statistics (Lagrangian correlations functions) computed from predictions of DNS. They are crucial for the analysis and the modelling of the fluid turbulent properties along discrete particle trajectories. In the second stage, a velocity-scalar Lagrangian stochastic model is proposed and evaluated from the DNS data. The coefficients of the drift and diffusion terms of the model are determined by only Lagrangian timescales, temperature variance and turbulent flux. The shapes of correlation functions present a good agreement between DNS results and stochastic modelling approach.     

Modeling plane turbulent Couette flow

Among the salient features of shear-driven plane Couette flow is the constancy of the total shear stress (viscous and turbulent) across the flow. This constancy gives rise to a quasi-homogenous core region, which makes the bulk of the flow substantially different from pressure-driven Poiseuille flow. The present second-moment closure study addresses the conflicting hypotheses relating to turbulent Couette flow. The inclusion of a new wall-proximity function in the wall-reflection part of the pressure-strain model seems mandatory, and the greement with recent experimental and direct numerical simulation (DNS) results is encouraging. Analysis of model computations in the range 750 ≤ Re ≤ 35,000 and comparisons with low-Re DNS data suggest that plane Couette flow exhibits a local-equilibrium core region, in which anisotropic, homogeneous turbulence prevails. However, the associated variation of the mean velocity in the core, as obtained by the model, conflicts with the intuitively appealing assumption of homogeneous mean shear. The constancy of the velocity gradient exhibited by the DNS therefore signals a deficiency in the modeled transport equation for the energy dissipation rate.     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

Large eddy simulation of turbulent statistical and transport properties in stably stratified flows

Three dimensional large eddy simulation （LES） is performed in the investigation of stably stratified turbulence with a sharp thermal interface. Main results are focused on the turbulent characteristic scale, statistical properties, transport properties, and temporal and spatial evolution of the scalar field. Results show that the buoyancy scale increases first, and then goes to a certain constant value. The stronger the mean shear, the larger the buoyancy scale. The overturning scale increases with the flow, and the mean shear improves the overturning scale. The flatness factor of temperature departs from the Gaussian distribution in a fairly large region, and its statistical properties are clearly different from those of the velocity fluctuations in strong stratified cases. Turbulent mixing starts from small scale motions, and then extends to large scale motions.     

Turbulence Budgets in Buoyancy-affected Vertical Backward-facing Step Flow at Low Prandtl Number

The paper investigates buoyancy impact on the vertical flow over a backward-facing step at low Prandtl number by Direct Numerical Simulation. In particular, the very low Prandtl number of liquid sodium, 0.0088, is considered in the regime of mixed convection, i.e. for Richardson numbers below unity. The effects of buoyancy on mean flow, heat transfer and turbulence are assessed. Buoyancy is found to attenuate recirculation and, consequently, increase heat transfer. Turbulence is decreased in the attached boundary layer for moderate buoyancy impact but surpasses the levels found in forced convection at the largest Richardson number investigated. Beyond the mean flow and second moments, the budgets of turbulent kinetic energy, Reynolds shear stress, temperature variance, and turbulent heat flux components are studied and related to the alterations in the mean field quantities. Due to scale separation, production and dissipation nearly balance for temperature variance while this is not the case for turbulent kinetic energy. Similar findings for the turbulent heat fluxes show that the correlation between temperature and pressure gradient is the most important contribution to the budget aside from production and dissipation. In addition to the physical insight into this flow, the data presented may be used for the validation and improvement of turbulence models for liquid metal flows.     

Direct numerical simulation of stable and unstable turbulent thermal boundary layers

This paper presents direct numerical simulations (DNS) of stable and unstable turbulent thermal boundary layers. Since a buoyancy-affected boundary layer is often encountered in an urban environmental space where stable and unstable stratifications exist, exploring a buoyancy-affected boundary layer is very important to know the transport phenomena of the flow in an urban space. Although actual observation may qualitatively provide the characteristics of these flows, the relevant quantitative turbulent quantities are very difficult to measure. Thus, in order to quantitatively investigate a buoyancy-affected boundary layer in detail, we have here carried out for the first time time- and space-developing DNS of slightly stable and unstable turbulent thermal boundary layers. The DNS results show the quantitative turbulent statistics and structures of stable and unstable thermal boundary layers, in which the characteristic transport phenomena of thermally stratified boundary layers are demonstrated by indicating the budgets of turbulent shear stress and turbulent heat flux. Even though the input of buoyant force is not large, the influence of buoyancy is clearly revealed in both stable and unstable turbulent boundary layers. In particular, it is found that both stable and unstable thermal stratifications caused by the weak buoyant force remarkably alter the structure of near-wall turbulence.     

可压缩均匀剪切湍流的二阶矩模拟

基于对可压缩湍流中脉动压力场和脉动速度场特征的理论分析以及ＤＮＳ结果,建立了可均匀剪切湍流中压力－变形率关联的压缩性修正模式,应用这个模式,加上Ｓａｒｋａｒ等建立的脉动体胀率项（ｄｉｌａｔａｔｉｏｎａｌ ｔｅｒｍｓ）的模式,预测可压缩均匀剪切湍流随时间的发展,所得雷诺应力各是性张量的平衡值与Ｂｌａｉｓｄｅｌｌ等的ＤＮＳ数据非常一致。这个模式准确地预测出均匀剪切湍流中压缩性导致的雷诺应力结构的“流向     

Direct numerical simulation of turbulent heat transfer in a wall-normal rotating channel flow

Investigations into the characteristics of turbulent heat transfer and coherent flow structures in a plane-channel subjected to wall-normal system rotation are conducted using direct numerical simulation (DNS). In order to investigate the influence of system rotation on the temperature field, a wide range of rotation numbers are tested, with the flow pattern transitioning from being fully turbulent to being quasilaminar, and eventually, fully laminar. In response to the Coriolis force, secondary flows appear as large vortical structures, which interact intensely with the wall shear layers and have a significant impact on the distribution of turbulence kinetic energy (TKE), turbulence scalar energy (TSE), temperature statistics, and turbulent heat fluxes. The characteristic length scales of turbulence structures responsible for the transport of TSE are the largest at the quasilaminar state, which demands a very large computational domain in order to capture the two-dimensional spectra of temperature fluctuations. The effects of the Coriolis force on the turbulent transport processes of the temperature variance and turbulent heat fluxes are thoroughly examined in terms of their respective budget balances.     

The influence of normal stress anisotropy in predicting scalar dispersion with the v–f model

A numerical study of scalar dispersion is presented to investigate the effectiveness of pairing the v²–f turbulence model with algebraic models for the scalar flux. This approach is contrasted with utilizing a full Second Moment Closure (SMC) as the flow field input to the scalar model. Predictions of scalar transport in a turbulent channel and over a wavy wall are compared to available DNS databases. The latter case includes a scalar release from a point source and therefore detailed comparisons of the three-component turbulent scalar flux are reported. It is found that the transported variable v², representing the near wall turbulent velocity fluctuation scale, can be used to increase the level of normal stress anisotropy provided to algebraic scalar models and thereby improve mean scalar prediction over that of the Standard Gradient Diffusion Hypothesis (SGDH). Improvement is most significant in the near wall region. Three specifications of the normal stresses, derived from v², are considered to provide the link from the v²–f model to the algebraic flux models used to close the scalar transport equation. Barycentric maps are used to examine the state of turbulence anisotropy in each case. As the anisotropy in the normal stress specification becomes more accurate, improvements are realized in the prediction of the spanwise flux as well as the mean concentration.     

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.     

Scalar transport in plane mixing layers

This paper describes the application of the Eulerian, single-point, single-time joint-scalar probability density function (PDF) equation for predicting the scalar transport in mixing layer with a high-speed and a low-speed stream. A finite-volume procedure is applied to obtain the velocity field with the k-ε closure being used to describe turbulent transport. The scalar field is represented through the modelled evolution equation for the scalar PDF and is solved using a Monte Carlo simulation. The PDF equation employs gradient transport modelling to represent the turbulent diffusion, and the molecular mixing term is modelled by the LMSE closure. There is no source term for chemical reaction as only an inert mixing layer is considered here. The experimental shear layer data published by Batt is used to validate the computational results despite the fact that comparisons between experiments and computational results are difficult because of the high sensitivity of the shear layer to initial conditions and free stream turbulence phenomena. However, the bimodal shape of the RMS scalar fluctuation as was measured by Batt can be reproduced with this model, whereas standard gradient diffusion calculations do not predict the dip in this profile. In this work for the first time an explanation is given for this phenomenon and the importance of a micromixing model is stressed. Also it is shown that the prediction of the PDF shape by the LMSE model is very satisfactory. Received on 27 October 1998     

压力梯度对湍流逆梯度输运的影响

研究Birkhoff系统的一般Lie对称性导致的非Noether守恒量. 得到非Noether守恒 量的存在定理，举例说明结果的应用.     

植被层湍流的大涡模拟

研究植被层湍流的大涡模拟,发展了一个TSF（transientstructurefunction）亚格于模式,尽可能真实地处理植被湍流这种既有强剪切,又有热对流的流动．我们建立了植被湍流数据库,并进行了较为详细的分析研究．湍流统计量如平均风速剖面、雷诺应力、湍流脉动能等等,与有关观测结果作了对比,符合较好．大涡模拟计算同样发现已由现场观测到的、在强对流情况时出现的温度场斜坡型有组织结构．