钝体燃烧器湍流预混燃烧的数值模拟

采用湍流火焰封闭燃烧模型(TFC)模拟了钝体燃烧器的湍流预混燃烧,比较了基于火焰褶皱率和湍流燃烧速度2种源项解法对钝体预混燃烧的预测,对3个不同湍流燃烧速度表达式模拟的性能进行了比较,采用粒子成像测速技术(PIV)测量了燃烧器中心射流出口的速度分布,并将其作为边界条件代入计算.结果表明:不同湍流燃烧速度公式的计算结果在火焰刷厚度、位置及火焰前锋位置方面存在较大差别;Gulder公式的计算结果最接近试验数据,火焰刷厚度与试验结果吻合较好,但火焰刷位置与试验结果差别较大;Dinkelacker的火焰褶皱率模型主要模拟燃烧器在高压条件下的燃烧,在运行压力接近标准大气压的情况下,计算结果与试验值存在较大误差.     

甲烷/空气预混射流火焰的大涡模拟

在不同射流速度条件下,对甲烷/空气预混射流火焰进行了大涡模拟.甲烷/空气预混射流气体按化学当量比混合,计算采用两步简化反应机理和WALE亚格子湍流输运模型,3个算例下流场特征和火焰结构计算结果与前人实验结果一致,中心线轴向速度和温度场结果与实验数据相符.通过对不同Karlovitz数条件下甲烷/空气预混射流火焰结构进行分析,并计算Takeno指数,研究了湍流涡对预混火焰的影响.研究发现:在Ka100(Ka=37)条件下,预混射流火焰会出现预热区的增厚,放热区保持完整,湍流火焰保持为预混燃烧;在Ka100(Ka=112)条件下,湍流火焰进入分布反应区模式.Takeno指数显示,由于卷吸和小尺度涡的作用,湍流火焰出现局部的部分预混燃烧.甲烷/空气预混射流湍流火焰的大涡模拟证实了湍流火焰分布反应区模式的特点:未燃气体与燃后气体之间不再有明显的界面,火焰面模型不再适用;反应区增厚,放热区展宽,放热率降低;由于卷吸和小尺度涡对火焰的作用,湍流火焰局部出现部分预混燃烧;湍流火焰温度降低,放热区附近温度场趋向均匀.     

基于CO-DAC反应机理简化的非预混射流火焰的大涡模拟

为了研究自相关自适应化学(CO-DAC)在湍流燃烧中的特性,计算了湍流非预混射流火焰Sandia FlameD的大涡模拟,使用了GRI-Mech 3.0详细化学反应机理以及CO-DAC简化方法.计算结果与实验值的对比表明,使用详细化学反应机理和CO-DAC机理简化,可以有效捕捉湍流火焰瞬时结构;可以根据燃烧场的特性,自适应地减少化学反应计算量,有效提高计算速度;可以准确模拟火焰的速度、温度分布特性;可以准确模拟火焰中重要中间组分和自由基的组分分布.     

稳态湍流非预混燃烧的小火焰模拟

以甲烷/空气的湍流射流非预混燃烧为对象,建立二维稳态湍流非预混火焰的小火焰模型.利用湍流流动模型和小火焰模型耦合求解,计算出速度、混合分数、温度以及反应标量的摩尔分数在燃烧室内的分布,模拟结果表明小火焰模型能够用来描述燃烧室内燃烧机理.     

煤粉复合旋流火焰燃烧特性研究（2）：数值模拟

本文数值模拟了煤粉旋流火焰燃烧过程，燃烧数值计算包括理论物理模型建立，数值方法两个大部分，计算模型处理了气相湍流与燃烧、气固两相流动、煤颗粒燃烧过程和辐射传热等物理化学过程，以ｋ－ε模型模拟湍流流动；ＰＤＦ法模拟气相扩散火焰燃烧；颗粒运动计算颗粒运动少颗粒湍流浓度方程模拟颗粒湍流扩散；通量法计算火焰辐射传热，煤粉颗粒复杂燃烧模型计算了颗粒尺寸、形状变化和颗粒孔隙内部燃烧、表面平度对整个颗粒的燃烧过程影响。计算获得了气相速度分布场、气相ｋ和ε分布场、气相温度场、气相组份场和颗粒浓度场及运动过程，揭示了煤粉复合旋流燃烧特性。     

配风比对微型燃气轮机燃烧室燃烧性能影响的数值模拟

以采用燃料和空气预混燃烧方式的微型燃气轮机燃烧室为研究对象,根据设计参数对燃烧室进行建模和模拟计算,模拟不同工况下预混燃料在燃烧室内经过湍流流动并发生燃烧化学反应的过程,进而得到燃烧室内的热态流场、温度分布以及出口烟气中污染物的排放量.结果 表明:在总过量空气系数为3.01的情况下,随着旋流器进口当量比的增大以及助燃风质量流量比例的升高,预混火焰的锋面温度有所升高,出口NOx质量浓度与出口温度分布因子呈正相关.     

600MW超临界W火焰锅炉燃烧特性及水冷壁分布参数混合模拟

以东方锅炉600MW超临界W火焰锅炉为研究对象,从非预混燃烧、气相湍流、颗粒相轨道模型、辐射传热、煤粉挥发分燃烧等模型入手,运用计算流体力学软件FLUENT进行了计算流体力学(CFD)模拟,同时建立了炉膛水冷壁一维分布参数模型,将其得到的水冷壁温度分布作为CFD模拟的边界条件,通过两种模型的混合模拟,得到了更为准确的基础工况和变负荷工况下炉膛内温度场、流场、组分浓度场的分布特性,并分析了炉内煤粉燃烧规律的变化。     

超声速湍流燃烧G/Z方程模型验证

为了研究超声速条件下的部分预混燃烧,引入一种基于Level set重构方法和稳态火焰面数据库的 G／Z方程模型,并利用德国宇航中心的DLR支板算例对 G／Z方程模型进行了验证。结果显示,超声速湍流燃烧G／Z方程模型可以捕捉到部分预混燃烧现象,数值模拟结果与实验结果吻合较好,验证了 G／Z方程模型运用到超声速部分预混条件下湍流燃烧流场计算的可行性。同时,超声速湍流燃烧 G／Z方程模型依赖于运用到的火焰传播速度模型与火焰面模型,模型精确度的提高有待进一步探究。     

一维湍流预混火焰的数值模拟

建立湍流燃烧的“双流体”数学模型，用于一维湍流预混稳态火焰的描述，假定燃烧火焰由冷的反应物（预混气体）和热的生成物（燃烧产物）组成，它们既有各自的属性，又相互作用，进行热量，质量和动量的交换，采用Patankar和Spalding的Phoenics计算程序来求解该数学模型，成功地模拟了一维湍流参混火焰的压力场，密度场，速度场。     

甲烷预混多喷嘴阵列燃烧器热声振荡模态实验研究

为研究燃气轮机模型燃烧室的非预混燃烧流场,采用大涡模拟方法分别结合火焰面生成流形模型（FGM）和部分预混稳态火焰面模型（PSFM）对甲烷/空气同轴射流非预混燃烧室开展了数值模拟研究,并与试验结果进行对比。结果表明：FGM所预测的速度分布、混合分数分布、燃烧产物及CO分布与试验结果更符合;两种模型均能捕捉到燃烧室中的火焰抬举现象;燃烧过程中的火焰结构较为复杂,同时存在预混燃烧区域和扩散燃烧区域,扩散燃烧主要分布在化学恰当比等值线附近,预混燃烧区域主要分布在贫油区。     

燃烧模型在NexGen燃烧器出口流场模拟中的应用

通过求解三维定常雷诺平均的N-S方程,对NexGen燃烧器的出口流场进行数值模拟。首先,利用冷态流场的实验数据验证燃烧器几何模型和数值计算方法的有效性。然后,在计算中分别选取火焰面模型、混合分数PDF模型和涡耗散模型3种燃烧模型,比较燃烧模型对燃烧器出口流场模拟结果的影响。研究结果表明:燃烧模型对Nex Gen燃烧器出口的速度场、火焰形状和热流密度分布基本没有影响,但是对火焰长度、火焰最高温度、最高热流密度及校准面上温度分布和温度值有较大影响。相比火焰面模型和涡耗散模型,混合分数PDF模型的计算结果与实验结果吻合较好,可以为防火试验方案设计提供参考。     

采用考虑详细化学反应机理的火焰面模型模拟湍流扩散火焰

采用详细的甲烷氧化化学反应动力学机理(GRI-Mech3．0)对不同拉伸率条件下的拉伸层流扩散火焰面结构进行了数值计算，建立了一个包含一系列拉伸层流火焰面结构的火焰面数据库．将这些层流火焰面结构和美国Sandia国家实验室测得的湍流扩散火焰(FlameD)的平均火焰结构进行了对比，发现层流火焰面所覆盖的范围基本包含了所考虑的湍流火焰中不同位置的平均火焰结构，这表明火焰面模型是合理的．然后，采用火焰面模型对该湍流扩散火焰进行了数值模拟并和实验数据进行了比较，考察了火焰面模型的精确程度和模拟深度.     

LES of a premixed jet flame DNS using a strained flamelet model

Turbulent premixed flames in the thin and broken reaction zones regimes are difficult to model with Large Eddy Simulation (LES) because turbulence strongly perturbs subfilter scale flame structures. This study addresses the difficulty by proposing a strained flamelet model for LES of high Karlovitz number flames. The proposed model extends a previously developed premixed flamelet approach to account for turbulence’s perturbation of subfilter premixed flame structures. The model describes combustion processes by solving strained premixed flamelets, tabulating the results in terms of a progress variable and a hydrogen radical, and invoking a presumed PDF framework to account for subfilter physics. The model is validated using two dimensional laminar flame studies, and is then tested by performing an LES of a premixed slot-jet direct numerical simulation (DNS). In the premixed regime diagram this slot-jet is found at the edge of the broken reaction zones regime. Comparisons of the DNS, the strained flamelet model LES, and an unstrained flamelet model LES confirm that turbulence perturbs flame structure to leading order effect, and that the use of an unstrained flamelet LES model under-predicts flame height. It is shown that the strained flamelet model captures the physics characterizing interactions of mixing and chemistry in highly turbulent regimes.     

Zone conditional assessment of flame-generated turbulence with DNS database of a turbulent premixed flame

Zone conditional two-fluid equations are derived and validated against a DNS database for a turbulent premixed flame. The conditional statistics of major flow variables are investigated to understand the mechanism of flame-generated turbulence. The flow field in the burned region shows substantially increased, highly anisotropic turbulence to conserve mass through a flamelet surface. The transverse component may be larger than the axial component for a distributed pdf of the flamelet orientation angle in the middle of the flame brush. The opposite occurs due to redistribution of turbulent kinetic energy and flamelet orientation mostly normal with respect to the mean flow at the end of the flame brush. The major source or sink terms of turbulent kinetic energy are the interfacial transfer by the mean reaction rate and the work terms induced by fluctuating pressure and velocity on the flame surface. Ad hoc modeling of some interfacial terms may be required for further application of the two-fluid model for modeling turbulence in turbulent premixed combustion simulations.     

Evaluation of the laminar diffusion flamelet model in the calculation of an axisymmetric coflow laminar ethylene-air diffusion flame

A numerical study of an axisymmetric coflow laminar ethylene-air diffusion flame at atmospheric pressure was conducted using detailed chemistry and complex thermal and transport properties and two different methodologies: (1) the direct simulation method of solving the two-dimensional axisymmetric elliptic governing equations, and (2) the steady-state stretched diffusion flamelet model. Soot formation and radiative heat transfer were not taken into account in these calculations, both for simplicity and to avoid the complications associated with the issues of how to incorporate these chemical and physical processes into the flamelet model. The same reaction mechanism and thermal and transport properties were used in the 2D direct simulation and the generation of the flamelet library. The flamelet library was generated from the solutions of counterflow ethylene-air diffusion flames at a series of stretch rates. Results of the 2D direct simulation and the flamelet model are compared in physical space. Although the overall results of the flamelet model are qualitatively similar to those of the direct simulation, significant differences exist between the results of the two methods even for temperature and major species. The direct simulation method predicts that the peak concentrations of CO₂ and H₂O occur in different regions in the flame, while the flamelet model results show that the peak concentrations of CO₂ and H₂O occur in the same region. The flamelet model predicts an overly rapid approach to the equilibrium structure in the downstream region, leading to significantly higher flame temperatures. The main reason for the failure of the flamelet model in the downstream region is due to the neglect of the effects of multidimensional convection and diffusion and the fundamental difference in the chemical structure between a coflow diffusion flame and a counterflow diffusion flame. The findings of this paper are highly relevant to understanding the flamelet model results in the calculations of multidimensional turbulent diffusion flames.     

Large-eddy simulation of turbulent autoigniting hydrogen lifted jet flame with a multi-regime flamelet approach

In this work, the combustion model is focused on to describe a multitude of reaction regimes that are deemed to affect the flame stabilization. For this purpose, an efficient flame indicator is formulated to differentiate the differing flame structures and make use of flamelet chemistry that accounts for autoignition and multi-regime reactions. The large eddy simulation with this methodology is carried out to compute a turbulent lifted hydrogen-nitrogen flame in vitiated coflow. The canonical flame models of a laminar premixed flame and an unsteady counterflowing flame have been used to simulate the flamelet structure at different regimes. Present model improves the prediction of mean and rms profiles for temperature and species mass fraction in the comparison with experiments and a reference simulation, adopting the single-regime flamelet. The computed results also demarcate the formation of a triple flame structure at the flame base, where combustion develops into the premixed reaction that extends to the fuel-lean and rich branches. The counterflow mixing mode with autoignition is identified as the major mechanism for stabilization and is responsible for the propagating premixed zone above the liftoff height. The developed multi-regime flamelet approach properly accounts for the interactive different modes of burning.     

Turbulent consumption speed via local dilatation rate measurements in a premixed bunsen jet

The mean local reaction rate related to the average expansion across the front and computed from the mean velocity divergence is evaluated in this work. Measurements are carried out in a air/methane premixed jet flame by combined PIV/LIF acquisitions. The procedure serves the purpose of obtaining values of a turbulent flame speed, namely the local turbulent consumption speed S_LC, as a function of the position along the bunsen flame. With the further position that the flamelet assumption provides a proportionality between turbulent burning speed normalized with the laminar unstretched one and the turbulent to average flame surface ratio, the proportionality constant, i.e., the stretching factor becomes available. The results achieved so far show the existence of a wide region along which the bunsen flame front has a constant stretching factor which apparently depends only on the ratio between turbulent fluctuations and laminar flame speed and on the jet Reynolds number.     

A new Lagrangian flamelet model for local flame extinction and reignition

Flamelet-based modeling for the evolution of chemically reactive scalars originally designed assuming large Damköhler numbers is extended for situations where local flame extinction and reignition occur. A Lagrangian extension of the flamelet model is formulated, considering a fluid particle that belongs to the same local flame structure as time evolves, where the temporal evolution of the flame structure is determined by, among other things, the scalar dissipation rate associated with the particle. The flamelet model, which is capable of predicting flame extinction but not reignition, is modified by mimicking known reignition scenarios. Model performance is examined using Lagrangian data obtained from direct numerical simulations of non-premixed reactions in isotropic turbulence; the results of the model agree both qualitatively and quantitatively with results obtained from the direct numerical simulations.     

甲烷/空气湍流扩散燃烧的小火焰模拟

以甲烷/空气的湍流射流扩散燃烧为基础,对通用的反应标量方程在火焰面上进行坐标变换,建立二维稳态湍流扩散火焰的小火焰模型。利用湍流流动模型、甲烷/空气半详细化学反应机理和小火焰模型耦合求解,分别计算出过量空气系数为1.2和1.4的速度在燃烧室内的分布状况以及混合分数、温度和组分的径向分布,模拟结果表明小火焰模型能够用来描述燃烧室内燃烧机理。     

Influence of evaporation on spray flamelet structures

The structure of laminar spray flames considerably differs from their gaseous counterpart. However, most often flamelet models employed in the simulation of turbulent spray combustion are based on laminar gas flame structures neglecting the influence of spray evaporation in the laminar spray flamelets. In this work, a combined theoretical and numerical study of the impact of spray evaporation on the structure of laminar spray flames is presented. Spray flamelet equations are derived, which explicitly take into account evaporation effects – the classical gas flamelet equations are recovered for non-evaporating conditions. Two new terms accounting for evaporation and for combined mixing and evaporation, respectively, are identified, and their relative importance is evaluated by means of numerical simulations of an axisymmetric laminar mono-disperse ethanol/air counterflow spray flame. The results show that the distribution of the spray evaporation rate plays a key role in the characterization of the spray flame structure. The new source terms overweigh the dissipation term of the gas phase in most situations even for non-evaporating species. Therefore, spray evaporation should always be considered. The relevance of the present formulation for turbulent spray modeling is evaluated and discussed, and a novel spray flamelet formulation is suggested.