惯性颗粒所见被动标量统计特性的直接数值模拟研究

本文采用直接数值模拟方法,对惯性颗粒所见均匀各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。模拟结果表明:与颗粒所见流体速度的自相关特性不同,颗粒所见标量的自相关特性随颗粒惯性的增加而单调减少;颗粒所见标量脉动能随颗粒惯性的增加先减少再增大,在St≈1.0的临界颗粒附近达到最小值,而颗粒所见标量耗散率随颗粒惯性的变化行为与颗粒所见标量脉动能的变化行为相反;数值模拟的结果进一步揭示,在St≈1.0时颗粒所见流体标量脉动能和耗散率的极值是因为St≈1.0的临界颗粒聚集于低涡量、高应变区域和标量场在高应变区域形成强耗散的片状结构所致。     

基于拉氏概率密度的两相湍流二阶矩模型

用颗粒运动的拉氏分析和PDF方法,改进了颗粒相的二阶矩模型。由拉氏两相运动的随机微分方程出发,采用随机过程分析和信号分析法得到湍流两相流动的PDF输运方程,双流体模型方程和两相脉动速度相关的基本模式的封闭式,和用其它方法导出的方程与封闭式的结果一致,对封闭式作了重要的改进,在分析颗粒轨道上的流体湍流作用时间时,全面地引入拉氏分析的轨道穿越效应、惯性效应、连续效应和湍流的各向异性。     

两相壁面射流PDF方程有限分析方法及数值模拟

提出求解位置-速度相空间中高维两相流PDF(probability density function)方程的有限分析方法,将位置-速度相空间颗粒PDF方程约化到速度空间,并解析求解,颗粒的位置PDF用轨道方法求解.对壁面射流两相流动进行数值模拟,并与颗粒雷诺应力轨道方法进行比较计算,结果优于颗粒雷诺应力轨道方法.     

自由剪切流动中颗粒扩散的自相似特性

为了研究自由剪切流动中颗粒扩散的统计特性,对空间发展模式的三维气固两相射流和时间发展模式的三维气固两相混合层进行了直接数值模拟.其中对气相不可压缩Navier-Stokes方程的求解分别采用有限容积方法和拟谱方法,对离散惯性颗粒的跟踪采用单向耦合的拉格朗日方法.统计结果显示,当流场进入自相似状态后,以流场速度梯度二阶不变量Q为自变量的颗粒数的概率分布函数也呈现出与颗粒尺寸和时间发展无关的自相似特性.     

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

两相湍流色噪声扩维法PDF模型及数值模拟

由颗粒运动的朗之万方程出发,对流体脉动速度采用扩维方法,得到两个不同层次的PDF输运方程.通过对颗粒运动方程求解和高斯分布假设,解决PDF方程的封闭问题,获得颗粒二阶矩模型,然后将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,对比分析数值结果与实验结果.     

非等温气固两相各向同性湍流的直接数值模拟

本文对非等温气固两相各向同性湍流中颗粒及颗粒所见流体温度的统计行为进行了直接数值模拟研究,系统地讨论了颗粒惯性的影响。对于Tp/Tk<1的颗粒,其温度脉动随惯性增大而减弱;对于Tp/Tk>1的颗粒,其趋势相反。颗粒及颗粒所见流体温度的自相关都随颗粒惯性的增加而减弱,而且随颗粒惯性的增加,颗粒所见流体温度的自相关比颗粒自身温度的自关联下降得快。     

颗粒惯性对标量混合的影响

本文采用直接数值模拟的方法,研究了颗粒惯性对均匀各向同性湍流场中被动标量混合特性的影响.结果表明:与无颗粒算例相比,Tp/Tk<1的小颗粒,增大了标量脉动在高波数上的能量,增强了标量耗散率,标量耗散率最大高出了25%;Tp/Tk>1的大颗粒,在初期显著削弱了标量脉动强度,同时也加快了标量耗散率的衰减,标量脉动强度降低了约15%,标量耗散减少了约5%;而Tp/Tk=1的临界颗粒,对标量场混合的影响较弱,标量脉动强度和标量耗散偏差在5%以内;此外,随着颗粒惯性的增大,标量耗散率的概率密度函数向左平移,但分布特征的改变并不显著.     

基于热质理论的热波传递现象研究

基于热质的概念应用牛顿力学原理建屯了热质输运的控制方程,推导得到了普适的导热定律.在惯性力可以忽略时它即退化为傅立叶导热定律,反映出非傅立叶导热现象的本质是热质的惯性引起的.当热流和温度对空间的惯性以及温度对时间的惯性可以忽略时,所得到的导热定律可以退化为CV模型.对热波传递的数值分析表明:当热扰动和热流都比较小时,热流在空间加速的惯性可以忽略,基丁热质理论的热波方程和CV模型符合得很好;但是,在描述较大的热扰动时,由于热流的空间加速惯性不能够被忽略,CV模型的结果在两个温度波峰叠加时会出现负温度的非物理现象,而基于热质理论的普适导热定律的结果则克服了这一缺陷.     

各向异性颗粒PDF输运方程及其颗粒湍流扩散

根据Minier的模型,获得颗粒所见流体脉动速度的朗之万模型,并推导颗粒所见流体湍流脉动速度自关联函数,该函数具有各向异性的特征,进一步由流体脉动速度朗之万模型和颗粒运动方程获得相应的颗粒PDF输运方程.为描述在网格生成湍流中的颗粒扩散运动,将颗粒PDF输运方程作相应的简化,获得一个具有颗粒湍流各向异性的解析解,将方程的解与Wellsand Stock的颗粒湍流扩散实验结果进行比较,获得了很好的结果.揭示了具有各向异性特征的颗粒PDF输运方程的优越性.     

Is turbulent mixing a self-convolution process?

Experimental results for the evolution of the probability distribution function (PDF) of a scalar mixed by a turbulent flow in a channel are presented. The sequence of PDF from an initial skewed distribution to a sharp Gaussian is found to be nonuniversal. The route toward homogeneization depends on the ratio between the cross sections of the dye injector and the channel. In connection with this observation, advantages, shortcomings, and applicability of models for the PDF evolution based on a self-convolution mechanism are discussed.     

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane–air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.     

Eulerian transported probability density function sub-filter model for large-eddy simulations of turbulent combustion

Reactive flow simulations using large-eddy simulations (LES) require modelling of sub-filter fluctuations. Although conserved scalars like mixture fraction can be represented using a beta-function, the reactive scalar probability density function (PDF) does not follow an universal shape. A one-point one-time joint composition PDF transport equation can be used to describe the evolution of the scalar PDF. The high-dimensional nature of this PDF transport equation requires the use of a statistical ensemble of notional particles and is directly coupled to the LES flow solver. However, the large grid sizes used in LES simulations will make such Lagrangian simulations computationally intractable. Here we propose the use of a Eulerian version of the transported-PDF scheme for simulating turbulent reactive flows. The direct quadrature method of moments (DQMOM) uses scalar-type equations with appropriate source terms to evolve the sub-filter PDF in terms of a finite number of delta-functions. Each delta-peak is characterized by a location and weight that are obtained from individual transport equations. To illustrate the feasibility of the scheme, we compare the model against a particle-based Lagrangian scheme and a presumed PDF model for the evolution of the mixture fraction PDF. All these models are applied to an experimental bluff-body flame and the simulated scalar and flow fields are compared with experimental data. The DQMOM model results show good agreement with the experimental data as well as the other sub-filter models used.     

Statistical properties of passive scalar in a random flow with a strong shear component

We study mixing of passive scalar by a chaotic velocity field with a relatively strong regular shear component. We show that the tail of partition distribution function (PDF) of coarse-grained passive scalar field differs qualitatively from the corresponding asymptotics in the case of isotropic flow statistics.     

Computations of turbulent lean premixed combustion using conditional moment closure

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k–? model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC.     

Finite-rate chemistry effects in turbulent opposed flows: comparison of Raman/Rayleigh measurements and Monte Carlo PDF simulations

This study reports results from experimental and numerical investigations of a partially premixed turbulent opposed methane/air jet flame. Experimentally determined properties of the scalar and the flow field are compared to the results from a Monte Carlo simulation. One-dimensional spatially resolved Raman/Rayleigh scattering serves to quantify the mean species concentrations and temperature, whereas laser Doppler velocimetry is used to measure axial and radial velocity components. The simulation is simplified by using a one-dimensional formulation. It includes a Reynolds-stress turbulence model and a Monte Carlo simulation of the joint scalar probability density function (PDF). A non-uniform Monte Carlo particle distribution is used to minimize stochastic errors. The flame is operated close to extinction with strong interactions between turbulence and chemistry. Comparisons between experimental and numerical results reveal a good agreement of mixture fraction profiles along the centreline. However, species scatter plots and mixture fraction PDFs show discrepancies between experiment and simulation. Numerical simulations over-predict the extinction limits and therefore under-predict the intermittent nature of turbulence and mixing of the scalars.     

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.