高温空气燃烧炉内湍流混合特性的数值研究

应用自行研发的三维流动、燃烧、传热和污染物NO_x湍流生成的数值模拟程序,对高温空气燃烧实验模型炉进行了湍流扩散燃烧混合特性的数值模拟.数值预报了燃烧室内气体燃料和空气的混合物分数及其湍流脉动的三维分布.数值研究结果表明:在一定的几何条件和气体动力学条件下,高温空气燃烧的湍流混合在更广泛的区域内以较小梯度的进行;混合物分数的脉动主要分布在燃烧区,这表明高温空气燃烧的火焰厚度更大,具有燃烧释热更趋均匀的特性.数值模拟结果与相关的实验结果有相同的规律.     

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

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

湍流分层燃烧的直接数值模拟和小火焰模型研究

湍流分层燃烧广泛应用于工业燃烧装置,但是目前还比较缺乏适用于湍流分层燃烧的高精度数值模型。本文利用直接数值模拟数据库,对高Karlovitz数分层射流火焰的小火焰模型表现进行了先验性评估。考虑了两种小火焰模型,一种是基于自由传播层流预混火焰的小火焰模型M1,另一种是基于分层对冲小火焰的小火焰模型M2。研究发现M1和M2在c-Z空间的结果与直接数值模拟在定性上是一致的。在物理空间,M2对过程变量反应速率脉动值的预测结果要优于M1.     

湍流扩散火焰的数值分析

本文基于GM80通用计算机程序运用K-ε-g湍流燃烧模型对丙烷-空气湍流扩散火焰进行数值分析,结果合理。 简要介绍了用湍流模型方法估计燃烧现象的基本思想和途径。建立了用于计算浓度脉动的子程序、解浓度脉动均方值方程的程序和便于输入非均匀性初值条件的程序,把这些并入GM80,用实际计算进行了检验。     

湍流脉动对自由火焰温度场反演的影响

发射-透射方法是火焰温度测量的一种重要方法。相对于湍流脉动的时间尺度而言,火焰红外光谱测量中探测器的时间分辨率一般较低。本文用迭代法分析了轴对称自由火焰温度的发射-透射测量中湍流脉动对温度场反演的影响。用简化的概率密度函数模拟温度和吸收系数的脉动。分析结果表明,在湍流脉动微弱的情况下,湍流脉动对时均温度场的反演影响很小,但在强湍流脉动的情况下湍流脉动对时均温度场的反演有很大影响。     

基于欧拉喷雾与火焰面/反应进度变量湍流燃烧模型的PRF燃料喷雾燃烧模拟研究

在开源计算流体软件OpenFOAM环境下,将基于欧拉方法的Σ-Y喷雾模型与非稳态火焰面/反应进度变量湍流燃烧模型相耦合,发展用于高温高压环境中液体燃料喷雾湍流燃烧高精度计算模型,分别对非燃烧和燃烧工况下的五种典型参比燃料(Primary Reference Fuel,PRF)的燃油喷射雾化与湍流燃烧过程开展数值研究。结果表明：该新型耦合模型能够准确的预测PRF燃料的喷雾和着火燃烧特性;所开发的重构数值喷雾纹影图像和燃烧OH^*图像处理方法能够很好地捕捉到试验的滞燃期、火焰浮起长度及喷雾火焰结构;研究揭示了不同比例的PRF燃料对喷雾及着火燃烧过程的影响特性,为替代燃料在发动机上的高效应用提供了理论指导。     

燃烧管内悬浮铝粉燃烧爆炸过程的研究

 铝粉的燃烧与加速机理一直是相关部门研究的热点。为了深入研究其在受限空间内的燃烧与传播特性，基于双流体模型、湍流k-ε模型以及铝粉燃烧等模型，选用SIMPLE格式，对大型卧式燃烧管内铝粉颗粒与空气的两相悬浮流湍流燃烧加速转爆炸现象进行了数值研究，得到了铝粉火焰在管内的燃烧传播过程中管内各相关参数的详细变化情况，并与相关的实验结果吻合。     

煤粉炉内辐射图象处理和燃烧过程数值模拟相结合初探

1引言大型电站锅炉炉膛内的燃烧过程是发生在相对较大空间内的、不断脉动的、具有强烈三维特征的复杂物理和化学过程。因为实际炉膛尺寸太大,以至于还没有建立适用的可视化技术手段[1],只能对其缩小了的模型在实验室进行研究[2.3]。作者在炉膛煤粉燃烧二维温度分布检测研究[4]的基础上,提出了以多幅辐射图像处理为基础的三维温度分布检测方法[5]。本文将建立较为严密的辐射图象信息同炉内燃烧过程的关系式,并借助燃烧过程数值模拟技术来估计炉内燃烧介质辐射特性参数非均匀分布,改进以辐射图象处理为基础的炉内三维燃烧温度分布检测方…     

采用一维湍流模型对瞬态煤粉气化火焰的数值模拟

煤的气化火焰中,湍流脉动与煤颗粒气化过程间存在着强烈的相互作用。为了在全尺度范围内直接模拟这种相互作用,本文采用一维湍流模型(ODT)与煤的气化过程相耦合,在Kolrnogorov尺度下对煤气化火焰区域进行数值模拟,得到了二维平面煤气化火焰的瞬态结构。模拟结果表明,大尺度涡团能够显著改变火焰的结构,并诱发局部小尺度涡团的产生。颗粒粒径决定湍流-气化过程作用的尺度范围,粒径较小的煤颗粒容易受到气体温度和速度脉动的影响,从而改变其运动轨迹和气化反应进程。     

湍流燃烧二阶矩模型的大涡模拟验证

对美国Sandia国家实验室测量的甲烷-空气湍流射流火焰进行了三维大涡模拟(LES),其统计的时均量以及脉动量符合实验测量值,且和RANS代数二阶矩湍流燃烧模型模拟结果大致符合。LES统计得到的化学反应率系数脉动和甲烷浓度脉动的互关联与对应的时均量梯度的乘积分布趋势相同,大小幅值的位置也相对应,与二阶矩湍流燃烧模型的封闭假设相符,验证了二阶矩代数燃烧模型的合理性。     

Investigation of non-premixed flame combustion characters in GO_2/GH_2 shear coaxial injectors using non-intrusive optical diagnostics

Single-element combustor experiments are conducted for three shear coaxial geometry configuration injectors by using gaseous oxygen and gaseous hydrogen(GO2/GH2) as propellants. During the combustion process, several spatially and time- resolved non-intrusive optical techniques, such as OH planar laser induced fluorescence(PLIF), high speed imaging, and infrared imaging, are simultaneously employed to observe the OH radical concentration distribution, flame fluctuations, and temperature fields. The results demonstrate that the turbulent flow phenomenon of non-premixed flame exhibits a remarkable periodicity, and the mixing ratio becomes a crucial factor to influence the combustion flame length. The high speed and infrared images have a consistent temperature field trend. As for the OH-PLIF images, an intuitionistic local flame structure is revealed by single-shot instantaneous images. Furthermore, the means and standard deviations of OH radical intensity are acquired to provide statistical information regarding the flame, which may be helpful for validation of numerical simulations in future. Parameters of structure configurations, such as impinging angle and oxygen post thickness, play an important role in the reaction zone distribution. Based on a successful flame contour extraction method assembled with non-linear anisotropic diffusive filtering and variational level-set, it is possible to implement a fractal analysis to describe the fractal characteristics of the non-premixed flame contour. As a result, the flame front cannot be regarded as a fractal object. However, this turbulent process presents a self-similarity characteristic.     

The impact of dilatation,scrambling, and pressure transport in turbulent premixed flames

Premixed turbulent flames feature strong interactions between chemical reactions and turbulence that affect scalar and turbulence statistics. The focus of the present work is on clarifying the impact of pressure dilatation/flamelet scrambling effects with a comprehensive second-moment closure used for evaluation purposes. Model extensions that take into account flamelet orientation and molecular diffusion are derived. Isothermal pressure transport is included with an additional variable density contribution derived for the flamelet regime of combustion. Full closure is assessed by comparisons with Direct Numerical Simulations (DNSs) of statistically ‘steady’ fully developed premixed turbulent planar flames at different expansion ratios. Subsequently, the prediction of lean premixed turbulent methane–air flames featuring fractal grid generated turbulence in an opposed jet geometry is considered. The overall agreement shows that ‘dilatation’ effects contribute to counter-gradient transport and can also increase the turbulent kinetic energy significantly. Levels of anisotropy are broadly consistent with the DNS data and key aspects of opposed jet flames are well predicted. However, it is also shown that complications arise due to interactions between the imposed pressure gradient and combustion and that redistribution is affected along with the scalar flux at the leading edge. The latter is strongly affected by the reaction rate closure and, potentially, by pressure transport. Overall, the derived models offer significant improvements and can readily be applied to the modelling of premixed turbulent flames at practical rates of heat release.     

Problems of predicting turbulent burning rates

Different approaches to the modelling of turbulent combustion first are reviewed briefly. A unified, stretched flamelet approach then is presented. With Reynolds stress modelling and a generalized probability density function (PDF) of strain rate, it enables a source term, in the form of a probability of burning function, P_b, to be expressed as a function of Markstein numbers and the Karlovitz stretch factor. When P_b is combined with some turbulent flame fractal considerations, an expression is obtained for the turbulent burning velocity. When it is combined with the profile of the unstretched laminar flame volumetric heat release rate plotted against the reaction progress variable and the PDF of the latter, an expression is obtained for the mean volumetric turbulent heat release rate. Through these relationships experimental values of turbulent burning velocity might be used to evaluate P_b and hence the CFD source term, the mean volumetric heat release rate.

Different theoretical expressions for the turbulent burning velocity, including the present one, are compared with experimental measurements. The differences between these are discussed and this is followed by a review of CFD applications of these flamelet concepts to premixed and non-premixed combustion. The various assumptions made in the course of the analyses are scrutinized in the light of recent direct numerical simulations of turbulent flames and the applications to the flames of laser diagnostics. Remaining problem areas include a sufficiently general combination of strain rate and flame curvature PDFs to give a single PDF of flame stretch rate, the nature of flame quenching under positive and negative stretch rates, flame responses to changing stretch rates and the effects of flame instabilities.     

Fractal structures in turbulence

We present a qualitative overview of our work on the issue of fractal structures in turbulence. We explain why fully developed turbulence is not space filling and describe how its fractal dimension can be estimated theoretically. The implications of the fractal nature of turbulence on transport processes like turbulent diffusion and on fluctuations in passive scalars are discussed. The latter affect wave propagation in turbulent media and these effects are examined. In addition we consider clouds in the atmosphere which are claimed to have fractal perimeters (or surfaces) and outline the physical reasons for this phenomenon. The fractal dimension of clouds is tied to the theory of turbulent diffusion and is computed theoretically. Indications of the road ahead are given.     

Modelling of the subgrid scale wrinkling factor for large eddy simulation of turbulent premixed combustion

We propose a model for assessing the unresolved wrinkling factor in the large eddy simulation of turbulent premixed combustion. It relies essentially on a power-law dependence of the wrinkling factor on the filter size and an original expression for the ‘active’ corrugating strain rate. The latter is written as the turbulent strain multiplied by an efficiency function that accounts for viscous effects and the kinematic constraint of Peters. This yields functional expressions for the fractal dimension and the inner cut-off length scale, the latter being (i) filter-size independent and (ii) consistent with the Damköhler asymptotic behaviours at both large and small Karlovitz numbers. A new expression for the wrinkling factor that incorporates finite Reynolds number effects is further proposed. Finally, the model is successfully assessed on an experimental filtered database.     

煤多相燃烧分形增长模型的初步研究

本文运用分形理论对煤的多相燃烧过程进行了探索性研究。根据分形动力学的概念认为；煤的多相燃烧中，反应面积的分形增长有DLC模式和KLC模式，一般情况是这两种模式的迭加。模型还引入了孔洞的合并及煤种的影响因素，初步提出了较有通用意义的燃烧模型。该模型把对碳焦的表面分形与燃烧速率的计算结合起来初步揭示了碳焦结构对燃烧过程的影响。     

The fractal measurement of experimental images of supersonic turbulent mixing layer

Flow visualization of supersonic mixing layer has been studied based on the high spatiotemporal resolution Nano-based Planar Laser Scattering(NPLS) method in SML-1 wind tunnel. The corresponding images distinctly reproduced the flow structure of laminar,transitional and turbulent region,with which the fractal measurement can be implemented. Two methods of measuring fractal dimension were introduced and compared. The fractal dimension of the transitional region and the fully developing turbulence region of supersonic mixing layer were measured based on the box-counting method. In the transitional region,the fractal dimension will increase with turbulent intensity. In the fully developing turbulent region,the fractal dimension will not vary apparently for different flow structures,which em-bodies the self-similarity of supersonic turbulence.     

Small scales,many species and the manifold challenges of turbulent combustion

A major goal of combustion research is to develop accurate, tractable, predictive models for the phenomena occurring in combustion devices, which predominantly involve turbulent flows. With the focus on gas-phase, non-premixed flames, recent progress is reviewed, and the significant remaining challenges facing models of turbulent combustion are examined. The principal challenges are posed by the small scales, the many chemical species involved in hydrocarbon combustion, and the coupled processes of reaction and molecular diffusion in a turbulent flow field. These challenges, and how different modeling approaches face them, are examined from the viewpoint of low-dimensional manifolds in the high-dimensional space of chemical species. Most current approaches to modeling turbulent combustion can be categorized as flamelet-like or PDF-like. The former assume or imply that the compositions occurring in turbulent combustion lie on very-low-dimensional manifolds, and that the coupling between turbulent mixing and reaction can be parameterized by at most one or two variables. PDF-like models do not restrict compositions in this way, and they have proved successful in describing more challenging combustion regimes in which there is significant local extinction, or in which the turbulence significantly disrupts flamelet structures. Advances in diagnostics, the design of experiments, computational resources, and direct numerical simulations are all contributing to the continuing development of more accurate and general models of turbulent combustion.     

A scale-entropy diffusion equation to describe the multi-scale features of turbulent flames near a wall

Multi-scale features of turbulent flames near a wall display two kinds of scale-dependent fractal features. In scale-space, an unique fractal dimension cannot be defined and the fractal dimension of the front is scale-dependent. Moreover, when the front approaches the wall, this dependency changes: fractal dimension also depends on the wall-distance. Our aim here is to propose a general geometrical framework that provides the possibility to integrate these two cases, in order to describe the multi-scale structure of turbulent flames interacting with a wall. Based on the scale-entropy quantity, which is simply linked to the roughness of the front, we thus introduce a general scale-entropy diffusion equation. We define the notion of “scale-evolutivity” which characterises the deviation of a multi-scale system from the pure fractal behaviour. The specific case of a constant “scale-evolutivity” over the scale-range is studied. In this case, called “parabolic scaling”, the fractal dimension is a linear function of the logarithm of scale. The case of a constant scale-evolutivity in the wall-distance space implies that the fractal dimension depends linearly on the logarithm of the wall-distance. We then verified experimentally, that parabolic scaling represents a good approximation of the real multi-scale features of turbulent flames near a wall.