湍流射流与扩散火焰大涡拟序结构的波动特性研究

１前言湍流射流扩散燃烧方式提供了射流火焰与涡团相互作用的最基本形式,对研究在湍流射流剪切边界层内的反应物的卷吸混合、热量及动量的输运及湍流射流结构方面具有其特殊的意义。Ｋａｔｔａ［１］利用浮力与单步反应机理模型对Ｎ２－Ｈ２射流扩散火焰进行了直接数值模...     

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

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

变密度随机涡模型在湍流射流扩散火焰中的应用

采用变密度随机涡模型,对H₂/O₂/N₂湍流射流扩散火焰进行数值模拟,湍流过程通过涡的采样、涡的抑制和涡的翻转实现.其中,针对变密度反应流问题,提出一种大涡抑制的新机制,并详细讨论各种参数对模型预测效果的影响.计算结果表明,修改后的模型可以合理预测H₂/O₂/N₂射流火焰结构,能够反映湍流的涡特性;模型中与涡采样和涡抑制有关的参数对预测结果有一定影响.     

介质阻挡体放电对甲烷扩散火焰特性的影响

本文在高频交流激励模式下,采用同轴圆柱构型激励器,开展了介质阻挡体放电对空气/甲烷同轴剪切扩散火焰燃烧特性影响实验研究。激励器敷设在外喷嘴环缝以电离空气,采用纹影系统和B型热电偶分别获取流场形态和火焰温度,激励频率为8 kHz,通过改变气体流量和放电电压,分析了不同工况下射流流场、火焰结构和火焰温度在等离子体作用下的变化规律。结果表明:等离子体气动效应能有效增强射流湍流强度,强化空气/甲烷掺混,增大射流角,并随激励电压提高作用效果逐渐增强,实验中未形成明显扩张流动的初始射流在放电电压30 kV时其射流角最大为23.5°。贫燃条件下等离子体激励会改善火焰形态,增强燃烧稳定性,并在流量较低时缩短火焰长度。此外,富燃火焰下游温度会随着激励强度增大不断升高,而贫燃火焰下游温度变化受上游燃烧强度影响存在升高和降低两种情况。     

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

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

两相扩散火焰结构探讨

利用鼓风式煤油喷灯,开展了不同空气射流流量下的两相扩散火焰结构的试验研究.采用工业相机及红外热像仪对稳定燃烧状态下的火焰参数进行了测量,获得了射流扩散火焰形状、高度和火焰温度随射流流量的变化规律.结果表明,射流扩散火焰呈现下小中间大上面小的形状,在顶端会出现分层或分叉现象;且随着射流流量的增加,火焰面呈现由光滑向褶皱再光滑的变化规律;火焰高度随着射流流量的增加,起初不断减小,直至基本保持不变,流量增大至180L/h时,出现了熄火;火焰最高温度随着空气射流流量的增加先升高后降低.     

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

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

拉伸流扩散火焰面结构及熄火的研究

对拉伸层流扩散火焰面进行了数值模拟,考察了在以往湍流燃烧的火焰面模型中,假定Lewis数等于1的可靠性,研究了不同分子扩散和火焰辐射对火焰面结构、氮氧化物排放和熄火临界的影响.计算结果表明,Lewis数等于1的假定在火焰面结构的计算中存在很大的近似性,火焰辐射可以引起低拉伸条件下的熄火临界.     

网格湍流对剪切层中预混火焰结构的影响

本文研究网格湍流对射流剪切层以及建立在其中的预混火焰的影响。利用热线风速仪测量射流的速度场,发现网格湍流使剪切层内湍流强度明显降低,抑制了低频速度脉动,同时增加了湍动能在小尺度脉动上的分配,使湍流更趋于各向同性,这表明网格湍流抑制了剪切层内的大涡和拟序结构。用细丝热电偶测量了火焰温度,结果显示网格湍流使火焰前峰的低频大幅摆动减少,小尺度皱褶增加,火焰区平均温度更高,说明网格湍流有利于剪切层中预混火焰的强化和稳定。     

湍流条件下煤粉颗粒群着火特性实验研究

本文搭建了平面火焰携带流反应器系统,研究了不同煤粉条件下煤粉颗粒群湍流射流着火特性。实验表明,湍流射流速度下低挥发份煤着火表现出扩散状暗红色火焰形态;随挥发分含量的增大,火焰亮度增加,并且在一定高度处发展为颗粒群群燃着火。随给粉量的增大,着火延迟呈现先减小后增大的趋势,最佳给粉浓度为0.75~1.00 kg/m~3。随粒径的减小,加热速率更快,挥发分大量析出使得火焰亮度增加;并且颗粒间的相互作用增强,更早出现颗粒群的群燃现象。     

Turbulent flame speed and morphology of pure ammonia flames and blends with methane or hydrogen

Ammonia appears a promising hydrogen-energy carrier as well as a carbon-free fuel. However, there remain limited studies for ammonia combustion especially under turbulent conditions. To that end, using the spherically expanding flame configuration, the turbulent flame speeds of stoichiometric ammonia/air, ammonia/methane and ammonia/hydrogen were examined. The composition of blends studied are currently being investigated for gas turbine application and are evaluated at various turbulent intensities, covering different kinds of turbulent combustion regimes. Mie-scattering tomography was employed facilitating flame structure analysis. Results show that the flame propagation speed of ammonia/air increases exponentially with increasing hydrogen amount. It is less pronounced with increasing methane addition, analogous to the behavior displayed in the laminar regime. The turbulent to laminar flame speed ratio increases with turbulence intensity. However, smallest gains were observed at highest hydrogen content, presumably due to differences in the combustion regime, with the mixture located within the corrugated flamelet zone, with all other mixtures positioned within the thin reaction zone. A good correlation of the turbulent velocity based on the Karlovitz and Damköhler numbers is observable with the present dataset, as well as previous experimental measurements available in literature, suggesting that ammonia-based fuels may potentially be described following the usual turbulent combustion models. Flame morphology and stretch sensitivity analysis were conducted, revealing that flame curvature remains relatively similar for pure ammonia and ammonia-based mixtures. The wrinkling ratio is found to increase with both increasing ammonia fraction and turbulent intensity, in good agreement with measured increases in turbulent flame speed. On the other hand, in most cases, the flame stretch effect does not change significantly with increasing turbulence, whilst following a similar trend to that of the laminar Markstein length.     

Capturing multi-regime combustion in turbulent flames with a virtual chemistry approach

Multiple flame regimes are encountered in industrial combustion chambers, where premixed, stratified and non-premixed flame regions may coexist. To obtain a predictive tool for pollutant formation predictions, chemical flame modeling must take into account the influence of such complex flame structure. The objective of this article is to apply and compare two reduced chemistry models on both laminar and turbulent multi-regime flame configurations in order to analyze their capabilities in predicting flame structure and CO formation. The challenged approaches are (i) a premixed flamelet-based tabulated chemistry method, whose thermochemical variables are parameterized by a mixture fraction and a progress variable, and (ii) a virtual chemical scheme which has been optimized to retrieve the properties of canonical premixed and non-premixed 1-D laminar flames. The methods are first applied to compute a series of laminar partially-premixed methane-air counterflow flames. Results are compared to detailed chemistry simulations. Both approaches reproduced the thermal flame structure but only the virtual chemistry captures the CO formation in all ranges of equivalence ratio from stoichiometry premixed flame to pure non-premixed flame. Finally, the two chemical models combined with the Thickened Flame model for LES are challenged on a piloted turbulent jet flame with inhomogeneous inlet, the Sydney inhomogeneous burner. Mean and RMS of temperature and CO mass fraction radial profiles are compared to available experimental data. Scatter data in mixture fraction space and Wasserstein metric of numerical and experimental data are also studied. The analyses confirm again that the virtual chemistry approach is able to account for the impact of multi-regime turbulent combustion on the CO formation.     

Premixed flames for arbitrary combinations of strain and curvature

Many modeling strategies for combustion rely on laminar flamelet concepts to determine structure and properties of multi-dimensional and turbulent flames. Using flamelet tabulation strategies, the user anticipates certain aspects of the combustion process prior to the simulation and selects a flamelet model which mimics local flame conditions in the more complex configuration. Flame stretch, which can be decomposed into contributions from strain and curvature, is one of the conditions influencing a flame’s properties, structure, and stability. The objective of this work is to study premixed flame structures in the strain-curvature space using a recently published composition space model (CSM) and three physical space models for canonical flame configurations (stagnation flame, spherical expanding flame and inwardly propagating flame). Flames with effective Lewis numbers both smaller and larger than unity are considered. For canonical laminar flames, the stretch components are inherently determined through boundary conditions and their specific flame configuration. Therefore, canonical flames can only represent a certain sub-set of stretch effects experienced by multi-dimensional and turbulent flames. On the contrary, the CSM allows arbitrary combinations of strain and curvature to be prescribed for premixed flames exceeding the conditions attainable with the canonical flame setups. Thereby, also influences of negative strain effects and large curvatures can be studied. A parameter variation with the CSM shows that flame structures still significantly change outside the region of the canonical flame configurations. Furthermore, limits in the strain-curvature space are discussed. The present paper highlights advantages of composition space modeling which is achieved by detaching the representation of the flame structure from a specific canonical flame configuration in physical space.     

CO2 addition and pressure effects on laminar and turbulent lean premixed CH4 air flames

An experimental study on CH₄–CO₂–air flames at various pressures is conducted by using both laminar and turbulent Bunsen flame configurations. The aim of this research is to contribute to the characterization of fuel lean methane/carbon dioxide/air premixed laminar and turbulent flames at different pressures, by studying laminar and turbulent flame propagation velocities, the flame surface density and the instantaneous flame front wrinkling parameters. PREMIX computations and experimental results indicate a decrease of the laminar flame propagation velocities with increasing CO₂ dilution rate. Instantaneous flame images are obtained by Mie scattering tomography. The image analysis shows that although the height of the turbulent flame increases with the CO₂ addition rate, the flame structure is quite similar. This implies that the flame wrinkling parameters and flame surface density are indifferent to the CO₂ addition. However, the pressure increase has a drastic effect on both parameters. This is also confirmed by a fractal analysis of instantaneous images. It is also observed that the combustion intensity S_T/S_L increases both with pressure and the CO₂ rate. Finally, the mean fuel consumption rate decreases with the CO₂ addition rate but increases with the pressure.     

An EDC-based turbulent premixed combustion model

In the present study, a new turbulent premixed combustion model is proposed by integrating the Coherent Flame Model with the modified eddy dissipation concept, and relating the fine structure mass fraction to the flame surface density. First, experimental results of turbulent flame speed available from literature are compared with the predicted results at different turbulence intensities to validate the flame surface density model. It is observed that the model is able to predict the turbulent burning speeds accurately. Then, a comprehensive validation is carried out utilizing data on a turbulent lifted methane flame issuing into a vitiated co-flow. Detailed comparison of temperature and species concentrations between experiment and simulation is performed at different heights of the flame. Overall, the model is found to predict both the spatial variation and peak values of the scalars at various heights satisfactorily.     

Conditional moment closure and transient flamelet modelling for detailed structure and NO x formation characteristics of turbulent nonpremixed jet and recirculating flames

This study has been mainly motivated to assess computationally and theoretically the conditional moment closure (CMC) model and the transient flamelet model for the simulation of turbulent nonpremixed flames. These two turbulent combustion models are implemented into the unstructured grid finite volume method that efficiently handles physically and geometrically complex turbulent reacting flows. Moreover, the parallel algorithm has been implemented to improve computational efficiency as well as to reduce the memory load of the CMC procedure. Example cases include two turbulent CO/H₂/N₂ jet flames having different flow timescales and the turbulent nonpremixed H₂/CO flame stabilized on an axisymmetric bluff-body burner. The Lagrangian flamelet model and the simplified CMC formulation are applied to the strongly parabolic jet flame calculation. On the other hand, the Eulerian particle flamelet model and full conservative CMC formulation are employed for the bluff-body flame with flow recirculation. Based on the numerical results, a detailed discussion is given for the comparative performances of the two combustion models in terms of the flame structure and NO _x formation characteristics.     

Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.     

湍流预混火焰传播速度的分形模型研究

本文利用激光层析技术和数字图像处理技术,在对Ｒｅｄ＝４３３５～１１１００范围内的Ｂｕｎｓｅｎ式湍流预混火焰热图像序列进行分形分析的基础上,提出了一种基于分形理论的湍流预混火焰传播速度模型,该模型将小尺度涡团在火焰锋面的强化湍流扩散效应归结为对锋面结构的改变上。结果表明：利用该模型预测的火焰传播速度与试验结果基本吻合。