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

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

湍流扩散火焰的非稳态火焰面模拟

采用稳态的和非稳态的火焰面模型同时对一个湍流甲烷射流扩散火焰进行了数值模拟,比较了两者对湍流平均火焰结构、活性自由基和污染物(氮氧化物)排放的模拟效果。速度场采用κ-ε模型计算,守恒标量混合物分数的分布通过其概率密度函数(PDF)输运方程的求解得到。稳态的火焰面结构由查询火焰面数据库得到,而非稳态的火焰面结构由火焰面方程和流场方程耦合求解来计算。采用详细的GRI—Mech 3.0机理描述甲烷的氧化和氮氧化物的形成。数值模拟结果和实验数据作了广泛的对比,验证了火焰面模型对湍流扩散燃烧的定量模拟能力。     

湍流扩散火焰局部熄火和再燃现象的PDF模拟

对一个值班湍流CH4／O2／N2射流扩散火焰(Sandia Flame D)进行了数值模拟研究．所采用的数学物理模型包括双尺度的k—ε湍流模型，标量联合的概率密度函数(PDF)输运方程方法，甲烷氧化的ARM简化化学反应机理(包含16种组分，12步总包反应)和欧几里德最小生成树(EMST)小尺度混合模型．将计算结果和实验数据进行了比较，不仅对于平均量，对于标量的散点分布和条件概率密度分布也是如此．计算结果表明文中采用的模型不仅能够预测宏观的火焰结构，而且预测了湍流燃烧中复杂的局部熄火和再燃过程．     

基于FGM的湍流射流扩散火焰超大涡模拟

为了研究超大涡模拟(VLES)预测湍流燃烧问题的能力,本文采用VLES结合基于假定概率密度函数的火焰面生成流型(FGM)建表湍流燃烧模型对值班甲烷/空气湍流射流扩散火焰(Sandia Flame D)开展了高精度数值研究,并与实验结果进行了详细比较。结果表明本文发展的VLES-FGM方法可以较准确地预测出湍流射流扩散火焰中的非稳态燃烧过程,且VLES湍流模拟方法对于湍流燃烧问题具有较大的应用潜力。     

抬举湍流H2/N2射流火焰的PDF模拟

采用数值目的研究了一个高温燃烧产物环境中的抬举湍流H₂/N₂射流火焰,对火焰的自然和抬举特性进行了研究.采用标量联合概率密度函数(PDF)目的处理详细的化学动力学过程,而湍流流场采用一个多时间尺度(MTS)k-ε湍流模型计算.计算中结合了一套描述氢气氧化的详细化学反应动力学机理.计算结果和实验数据进行了对比,表明所采用的模型可以精确的模拟火焰抬举高度和自然的过程.     

受限湍流射流扩散火焰的PDF模拟

采用k-ε双方程模型和概率密度函数(PDF)相结合的办法,研究受限条件下的湍流射流扩散火焰,着重考虑受限条件下固体壁面、压力梯度等因素对速度场和标量场求解的影响,并在此基础上对两个不同尺寸的受限湍流燃烧场进行计算,分别研究了受限湍流射流扩散火焰的流场结构、火焰结构和火焰形状.最后给出结果,定性分析,并得出结论.     

PDF方法模拟钝体驻定的湍流扩散火焰

采用标量联合的概率密度函数方法,对钝体驻定的湍流射流扩散Sydney火焰HM1进行数值模拟,结合当地自适应建表方法加速化学反应计算,用修正的LRR-IP雷诺应力模型求解速度场.首次对3种不同规模的甲烷化学反应动力学机理进行研究,并与实验数据进行比较,结果表明,模型和反应机理很好地预测了速度场和标量场的变化及局部熄火现象,而考虑反应机理中的C2化学对火焰HM1的影响不大.     

甲醇喷雾湍流反应流的JPDF数值模拟

喷雾湍流燃烧过程中,液滴、湍流和化学反应之间强烈耦合,物理化学机理非常复杂。本文将速度-标量-频率联合概率密度函数JPDF输运方程方法应用于两相喷雾湍流反应流问题,利用火焰面模型解耦流动和化学反应动力学的耦合关系,建立起相应的数值计算模型。采用Monte-Carlo数值计算方法,针对澳大利亚悉尼大学Masri等人以甲醇为燃料所进行的湍流喷雾燃烧值班火焰这一试验进行了数值模拟,通过与Fluent下的计算结果及试验结果的对比分析,验证了本文所建模型的准确性。     

轴对称钝体后湍流扩散燃烧的PDF模拟

在无结构网格中，对轴对称钝体驻定的湍流扩散火焰进行数值模拟．采用有限容积(FV)／Monte Carlo(MC)混合算法求解湍流燃烧问题的混合算法．Monte Carlo法求解脉动速度一标量．频率的联合概率密度函数方程，有限容积法求解平均质量、动量和能量方程．求解的两组方程是相容的，合理的耦合方式可以减少统计偏差，计算精度和效率显著优于单独的颗粒方法．文中对化学反应采用层流火焰面模型，并将数值计算结果与实验结果作了比较和分析．     

氢气扩散火焰中辐射源项湍流脉动特征的PDF模拟

采用κ-ε湍流模型、标量联合的概率密度函数(PDF)输运方程和层流火焰面模型相结合,模拟氢气自由扩散火焰中辐射源项湍流脉动特征.给出了主燃区内辐射源项湍流脉动的频率图.辐射源项的样本点分布集中,大约95%以上的样本落在其系综的±3倍方差以内,频谱图为单峰.     

A comparison of different approaches to integrate flamelet tables with presumed-shape PDF in flamelet models for turbulent flames

Flamelet models for turbulent combustion modelling make use of presumed-shape probability density functions (PDFs) for integrating laminar flamelet solutions to obtain an integrated flamelet table that can readily be used for turbulent flame calculations. The existence of non-unique approaches for such an integration has rarely been investigated before. For the first time, this work studies systematically the non-uniqueness of the flamelet table integration approaches. A flamelet model called the flamelet/progress variable model is used in the study, although the issue exists generally in many other flamelet models. Two classes of table integration approaches are investigated, one preserving the laminar flamelet structures during integration and the other not. Three different table integration approaches are examined and compared in detail to provide a thorough understanding of the different approaches. A partially stirred reactor is used as a test case for examining the different approaches. A method based on the transported PDF method is also employed to provide a reference for the assessment of the different flamelet table integration approaches. It is found in general that the flamelet preserving integration approach yields a more reasonable joint PDF of the mixture fraction and the progress variable, and the prediction results are closer to the referenced transported PDF results.     

Prediction of NO in turbulent diffusion flames using Eulerian particle flamelet model

The combustion characteristics for the turbulent diffusion flames using the unsteady flamelet concept have been numerically investigated. The Favre-averaged Navier–Stokes equations are solved by a finite volume method of SIMPLE type that incorporates the laminar flamelet concept with a modified k ? ε turbulence model. The NO formation is estimated by solving the Eulerian particle transport equations in a postprocessing mode. Two test problems are considered: CH₄/H₂/N₂ jet flame and CH₄/H₂ stabilised bluff body flame. The temperature and species profiles are well captured by the flamelet model. Two different chemical mechanisms (GRI 2.11 and 3.0) give nearly identical results for temperature and species except NO. The GRI 3.0 gives significantly higher NO levels compared to the GRI 2.11. This is mainly attributed to the difference in NO formation by the prompt mechanism. The NO formation is sensitive to the number of flamelet particles. The NO levels for two test flames do not change when the flamelet particle number exceeds six.     

Flamelet modeling of coal particle ignition

This paper numerically investigates the ignition of a single coal particle during the devolatilization phase in a laminar entrained-flow reactor, for which experimental data are available from Molina and Shaddix [3]. Different numerical approaches are combined to evaluate the non-premixed flamelet approach for coal particle ignition. First, the particle trajectory and the particle heating are simulated with a Lagrangian–Eulerian approach using a detailed pyrolysis model. In a second step, these results are used as transient boundary conditions for a simulation fully resolving the flow, the mixing field and the chemical reactions around the particle. Finally, in combination with the boundary conditions the time-dependent scalar dissipation rate profiles from the resolved particle calculation are used in a flamelet calculation for the particle up- and downstream directions. Very good agreement is obtained in terms of ignition delay as well as temperature and chemical species distributions in the mixture fraction space when the resolved particle calculation and the unsteady flamelet calculation are compared in the downstream direction. Good agreement is obtained when the numerical results for the ignition time and the time-averaged OH distribution are compared with the available experimental data. The results show the capability of the laminar flamelet approach to correctly predict coal particle ignition during devolatilization using accurate scalar dissipation rate profiles.     

小火焰模型在贫燃预混火焰中的研究

由层流小火焰库引入详细化学反应机理,通过简化的PDF方法计算组分浓度、平均温度和密度等变量,以钝体火焰稳定燃烧室和某燃气轮机上的燃烧室为例,模拟甲烷/空气贫燃条件下预混燃烧的平均火焰位置和火焰厚度,计算结果与实验结果吻合良好,这表明此方法能够较好计算出平均湍流火焰的主要特征。     

A quadrature-based LES/transported probability density function approach for modeling supersonic combustion

The joint-scalar probability density function (PDF) approach provides a comprehensive framework for large eddy simulation (LES) based combustion modeling. However, currently available stochastic approaches for solving the high-dimensional PDF transport equation can be error prone and numerically unstable in highly compressible shock-containing flows. In this work, a novel Eulerian approach called the direct quadrature method of moments (DQMOM) is developed for evolving the PDF-based supersonic combustion model. The DQMOM technique uses a set of scalar transport equations with specific source terms to recover the PDF. The new technique is coupled to a compressible LES solver through the energy equation. The DQMOM approach is then used to simulate two practical flow configurations: a supersonic reacting jet and a cavity-stabilized supersonic combustor. Comparisons with experimental data demonstrate the predictive accuracy of the method.     

PDF mixing time scales for premixed combustion in the laminar flame limit

A new mixing time scale for PDF calculations of premixed combustion in the laminar flame limit is introduced. It is based on the characteristics of the “random walk” diffusion process and accurately captures scalar micro-mixing such that physical features of the premixed flames are preserved. The speed of a freely-propagating laminar flame can be captured accurately. Extreme stochastic events may, however, lead to flame acceleration, but flame stability can be recovered with the introduction of a conditioning variable that relaxes towards an Eulerian reference field. With conditioning the predicted flame speed is quite insensitive towards the exact value of the only modelling parameter. This modelling parameter then allows to control the minor scalar fluctuations and deviations from a flamelet structure which is thought to be one of the key features of conditioning methods.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.