旋流燃烧NO生成的USM湍流反应模型

用统一二阶矩(USM)湍流反应模型对不同旋流数下甲烷-空气旋流燃烧NO生成进行了数值模拟,并和EBU-Arrhenius(E-A)燃烧模型对燃烧的模拟结果和简化PDF湍流反应模型对NO生成的模拟结果以及相应的实验结果进行对比。结果表明,USM模型显著地优于E-A模型和简化PDF模型。E-A模型不能合理地模拟有限反应动力学,而简化PDF模型用两个单变量PDF的乘积代替联合PDF,大大地低估了NO湍流反应率。USM模型预报结果和实验结果符合最好。     

湍流燃烧的新二阶矩模型

本文针对湍流燃烧中模拟包括污染物生成的详细有限反应动力学的要求，用二阶矩封闭与简化概率密度的概念相结合的方法，提出了湍流燃烧的新二阶矩模型。此模型可兼顾适用于工程问题中的合理性和经济性。     

柴油机氮氧化物的特征时间模型

本文综合考虑湍流混合、分子扩散的作用,提出了一个新的NO生成计算模型.该模型采用化学动力学时间尺度和湍流时间尺度来考虑化学动力学和湍流混合对NO生成的影响.文中模拟了一台浅盆形燃烧室柴油机,对计算结果和实测结果进行了对比分析.研究表明,新模型计算的NO生成速率与实测结果有较好的一致性.     

旋流对煤粉燃烧NO排放影响的数值模拟

本文用基于HCN释放的简化Solomon模型的NO生成湍流反应的统一二阶矩代数模型(AUSM)和煤粉燃烧的双流体模型,对不同旋流数下煤粉燃烧器内两相流动,煤粉燃烧和NO生成进行了数值模拟。模拟结果和文献中实验结果符合很好。模拟结果指出,随着旋流数的增加,NO的排放先减少后增加,燃尽率先增加后减小,和气体燃烧中得到的规律类似。     

非预混燃烧大涡模拟的代数二阶矩模型

一种改进的代数二阶矩燃烧模型被应用在Flame-D的大涡模拟中。代数二阶矩模型能够将化学反应速率的脉动项进行模拟,更准确地计算化学反应速率。在实际燃烧过程中,化学反应速率不仅与当地的温度以及组分浓度等相关,更与反应物的混合速率有密切的关系。研究考虑了亚网格尺度内混合速率对化学反应速率的影响,提出了改进的代数二阶矩模型。通过与相同模拟条件下的涡耗散概念模型的模拟结果以及实验结果进行对比,改进的代数二阶矩模型的准确性得到了验证,同时亚网格化学反应速率的模拟对总化学反应速率的影响也得到了研究。     

湍流燃烧二阶矩模型的研究进展

湍流燃烧模型在燃烧过程数值模拟中十分重要。商业软件中仍然应用的简单模型,如EBU和预设PDF模型,常常不能很好地模拟有限反应动力学。目前通行的湍流燃烧模型,如层流小火焰模型和条件矩模型,只对一定的火焰类型和火焰结构的效果较好。PDF方程模型更通用,但计算量太大,用于大涡模拟更是如此。另一类是统计矩模型,即基于湍流模型的思路,用雷诺展开和取平均,封闭未知项的二阶矩模型,但是遇到了高度非线性的温度指数函数的困难。不同研究者采取了不同的近似处理,都低估了时平均反应率。作者彻底放弃各种近似方法,构建了终版的二阶矩模型,用于不同的单相和两相燃烧的雷诺平均模拟和大涡模拟,得到了实验验证和直接数值模拟的验证。     

突扩燃烧室湍流预混反应流的数值模拟

对突扩燃烧室这一典型工程燃烧装置内的湍流预混反应流进行了数值模拟。时平均控制方程组的封闭采用k-ε湍流输运模型和EBU-Arhenius湍流反应模型。模拟结果给出了突扩燃烧室内湍流预混反应流的气体时均流场、组分浓度场与温度场的分布。通过数值模拟结果与实验的比较对EBU-Arhenius模型进行了讨论与评价。     

湍流两相流动有燃烧颗粒相概率密度函数输运方程理论

由有燃烧的湍流气粒两相流动的瞬态方程和统计力学概率密度函数概念出发,推导了有燃烧颗粒相的质量－动量－能量联合概率密度函数（ＰＤＦ）输运方程,并对方程中条件期望项用梯度模拟概念进行了模拟封闭。封闭后的ＰＤＦ方程可作为建立颗粒拟流体模型方程和封闭二阶矩模型的基础,也可以通过Ｍｏｎｔｅ－Ｃａｒｌｏ 法求解用以直接计算颗粒雷诺应力和湍流动能,以便和二阶 矩模型的结果相对照,改善二阶矩模型。     

几种二阶矩模型模拟湍流旋流燃烧的比较

本文对几种二阶矩湍流反应模型,包括统一二阶矩模型、只考虑温度脉动的二阶矩模型、只考虑浓度脉动关联的二阶矩模型、及同时考虑温度脉动和浓度脉动关联的二阶矩模型,进行了比较。将上述模型加入到FLUENT6.0软件平台上,模拟了不同旋流数下甲烷-空气的旋流燃烧。模拟结果和实验结果进行了比较,探讨了各关联量大小及其对时平均反应率的影响。结果表明,统一二阶矩湍流反应模型具有最好的模拟效果。其原因是,各种关联矩中,化学反应率系数与浓度间的脉动关联最重要。     

高强度湍流预混火焰NO生成特性的研究

本文采用直接数值模拟方法对不同湍流强度下预混CH_4/H_2/air平面火焰进行了计算,并对湍流火焰中NO的生成规律进行了研究。采用了28组分268步反应的详细化学反应机理。研究表明高Karlovitz数下,湍流微团能穿透火焰面并存在于内反应区中。在强拉伸作用下,局部燃烧热释率达到层流火焰峰值的1.8倍。湍流作用下湍流燃烧速率显著提高,而NO的总生成速率增加不明显。高湍流强度下由于快速型NO的生成率的显著减小,湍流火焰中NO的总生成速率降低。     

旋流数对湍流燃烧中NO生成影响的研究

通过实验和用湍流燃烧二阶矩概率密度模型对不同旋流数下甲烷-空气湍流燃烧和NO生成进行了研究。在燃料中加入少量氨模拟燃料氮。研究结果表明,随旋流数的增大（由0到1）,热NO排放先上升后下降,而总NO和燃料NO排放则先下降后上升。旋流数增大使湍流强度先下降然后稍有上升,使进口附近温度先上升然后稍有下降。热NO的生成受温度的影响更大而燃料NO的生成受湍流的影响更大。     

水泥回转窑内NO生成及其分布规律研究

本文在气固流动、煤粉燃烧和NO生成数学模型的基础上,对水泥回转窑内物料烧成过程的物理化学反应热效应采用分区段拟合的方法,建立了一套描述水泥回转窑窑内过程的数学模型。并对某3000吨／天生产能力的带四通道燃烧器的水泥回转窑进行了数值模拟,得到了回转窑内气体速度场、气体温度和组分浓度沿窑长的变化规律,对窑内NO生成进行了深入研究。研究结果表明:水泥回转窑内NO生成按机理可分为热力型NO和燃料型NO,由于窑内存在着高温、富氧环境,热力型NO为主要生成方式;热力型NO和燃料型NO生成过程存在着相互抑制作用。     

Prediction of soot and thermal radiation in a model gas turbine combustor burning kerosene fuel spray at different swirl levels

Combustion of kerosene fuel spray has been numerically simulated in a laboratory scale combustor geometry to predict soot and the effects of thermal radiation at different swirl levels of primary air flow. The two-phase motion in the combustor is simulated using an Eulerian–Lagragian formulation considering the stochastic separated flow model. The Favre-averaged governing equations are solved for the gas phase with the turbulent quantities simulated by realisable k–? model. The injection of the fuel is considered through a pressure swirl atomiser and the combustion is simulated by a laminar flamelet model with detailed kinetics of kerosene combustion. Soot formation in the flame is predicted using an empirical model with the model parameters adjusted for kerosene fuel. Contributions of gas phase and soot towards thermal radiation have been considered to predict the incident heat flux on the combustor wall and fuel injector. Swirl in the primary flow significantly influences the flow and flame structures in the combustor. The stronger recirculation at high swirl draws more air into the flame region, reduces the flame length and peak flame temperature and also brings the soot laden zone closer to the inlet plane. As a result, the radiative heat flux on the peripheral wall decreases at high swirl and also shifts closer to the inlet plane. However, increased swirl increases the combustor wall temperature due to radial spreading of the flame. The high incident radiative heat flux and the high surface temperature make the fuel injector a critical item in the combustor. The injector peak temperature increases with the increase in swirl flow mainly because the flame is located closer to the inlet plane. On the other hand, a more uniform temperature distribution in the exhaust gas can be attained at the combustor exit at high swirl condition.     

旋流扩散燃烧中旋流数对热NO生成的影响

本文对旋流扩散燃烧进行了数值模拟，研究旋流数对热ＮＯ生成的影响，其中对湍流采用Ｒｅｙｎｏｌｄｓ应力方程模型，对燃烧采用ＥＢＵ－Ａｒｒｈｅｎｉｕｓ模型，对热ＮＯ生成采用设定ＰＤＦ的模型。预报了不同旋流数下轴向和切向的平均和脉动速度、温度和ＮＯ浓度，指出随着旋流数的增大，计算得到的出口平均ＮＯ浓度首先升高然后下降。这一趋势和本文作者最近的实验结果的趋势一致．随着旋流数的增大，湍流脉动首先下阵然后升高，而进口附近的温度上升，二者综合效果造成上述趋势、因此在实际燃烧器中，完全靠改变族流数来控制ＮＯ生成是不现实的。应该采取其他方法来降低ＮＯ的生成。     

Conditional moment closure modelling of turbulent jet diffusion flames of helium-diluted hydrogen

First-order conditional moment closure (CMC) modelling of NO in non-premixed flames has met with limited success due to the need to consider turbulence influences on the conditional production rate of chemical species. This paper presents results obtained using a second-order approach where such effects are incorporated through solution of a transport equation for the conditional variance. In contrast to earlier work, second-order chemistry is implemented using a more robust numerical technique, with predictions obtained using a Reynolds stress turbulence model. First-order CMC and k–? turbulence model predictions are presented for comparison purposes. For the hydrogen flames examined, results demonstrate small differences between first- and second-order calculations of major species and temperature, although second-order corrections reduce NO and OH levels. Additionally, variations occur between results for these species derived using the two turbulence models due to differences in conditional variance predictions. This and the numerical solution method employed are responsible for deviations with earlier results. It is concluded that while the higher-order CMC model does not significantly improve NO predictions, agreement with OH data is superior. Physical space predictions are sufficiently accurate for assessing flame characteristics, with the Reynolds stress model providing superior results.     

Modeling of a self-excited pulse combustor and stability analysis

The major bottleneck for popularization and utilization of the conventional mechanical valve pulse combustors is the self-priming mode of gas supply. An aerodynamic valve (as against mechanical valve) self-excited pulse combustor of the Helmholtz-type with continuous supply of gas and air was designed and a mathematical model was established in this paper. The theoretical model employed well-stirred reactor model and a single step Arrhenius chemistry, and took those factors which might affect the combustion stability into account. The factors include the variation of the mass rate of the reactants affected by the pressure in the combustion chamber, the convective and radiation heat loss in the combustion chamber, and the heat transfer and wall friction in the tailpipe. The effect of wall temperature of combustion chamber, wall heat transfer coefficient, tailpipe length and friction coefficient on combustionstability were analyzed. The range of combustion oscillations can be predicted. It is theoretically and experimentally shown that combustion oscillations can be produced with a continuous supply of fuel and air without mechanical valves. The experimental data show qualitative agreement with predictions from the theoretical model. The theoretical model could be a tool for designing and optimizing the self-excited pulse combustor.     

Simulation of a round supersonic combustor using wall-modeled large eddy simulation and partially-stirred reactor models

Simulations are presented for a generic, round supersonic combustor. Turbulence is modeled in the combustor using a wall-modeled large eddy simulation approach. Combustion is modeled using a small quasi-global mechanism and a more detailed skeletal mechanism. Both mechanisms are used in conjunction with two variations of the partially-stirred reactor model for sub-grid turbulence chemistry interactions. Sensitivity of the solutions to grid resolution is investigated. It is found that in order to achieve reasonable grid convergence in the mean wall pressure, the model constant that appears in the partially-stirred reactor model must be a function of both the chemistry mechanism and the grid resolution. Most of the combinations of mechanism and turbulent combustion model tested can be tuned in order to predict the location of the pre-combustion shock train and the peak mean pressure in the combustor. It is found that while the different models are able to reproduce the mean wall pressure, there are significant differences in the mean temperature and heat release rate fields. The sensitivity of the different combinations of mechanisms and partially-stirred reactor formulation is quantified and some combinations are found to be more prone to blowout. Two of the tuned models were tested across several fuel equivalence ratios with a single value of the partially-stirred reactor model constant. One model produced reasonable predictions of shock location and peak mean pressure for each equivalence ratio. The second model captured the global trends in the mean wall pressure, but was unable to quantitatively predict the shock location and peak mean pressure for all equivalence ratios tested.     

Spherical turbulent flame propagation of pulverized coal particle clouds in an O2/N2 atmosphere

The present study aims to clarify the effects of turbulence intensity and coal concentration on the spherical turbulent flame propagation of a pulverized coal particle cloud. A unique experimental apparatus was developed in which coal particles can be dispersed homogeneously in a turbulent flow field generated by two fans. Experiments on spherical turbulent flame propagation of pulverized coal particle clouds in a constant volume spherical chamber in various turbulence intensities and coal concentrations were conducted. A common bituminous coal was used in the present study. The flame propagation velocity was obtained from an analysis of flame propagation images taken using a high-speed camera. It was found that the flame propagation velocity increased with increasing flame radius. The flame propagation velocity increases as the turbulence intensity increases. Similar trends were observed in spherical flames using gaseous fuel. The coal concentration has a weak effect on the flame propagation velocity, which is unique to pulverized coal combustions in a turbulent field. These are the first reports of experimental results for the spherical turbulent flame propagation behavior of pulverized coal particle clouds. The results obtained in the present study are obviously different from those of previous pulverized coal combustion studies and any other results of gaseous fuel combustion research.     

Effect of cavity location on combustion flow field of integrated hypersonic vehicle in near space

Abstract  

The cavity has been widely employed as the flame holder to prolong the residence time of fuel in supersonic flows since it improves the combustion efficiency in the scramjet combustor, and also imposes additional drag on the engine. In this paper, the two-dimensional coupled implicit Reynolds Average Navier–Stokes equations, the RNG k–ε turbulence model and the finite-rate/eddy-dissipation reaction model have been employed to numerically simulate the combustion flow field of an integrated hypersonic vehicle. The effect of cavity location on the combustion flow field of the vehicle has been investigated, and the fuel, namely hydrogen, was injected upstream of the cavity on the walls of the first stage combustor. The obtained results show that the viscous lift force, drag force and pitching moment of the vehicle are nearly unchanged by varying the cavity location over the location range and designs considered in this article, namely the configurations with single cavity, double cavities in tandem and double cavities in parallel. The variation of the fuel injection strategy affects the separation of the boundary layer, and the viscous effect on the drag force of the vehicle is remarkable, but the viscous effects on the lift force and the pitching moment are both small and they can be neglected in the design process of hypersonic vehicles. In addition to varying the location of the cavities, three fuel injection configurations were considered. It was found that one particular case can restrict the inlet unstart for the scramjet engine.     

Simultaneous OH-PLIF and PIV measurements in a gas turbine model combustor

In highly turbulent environments, combustion is strongly influenced by the effects of turbulence chemistry interactions. Simultaneous measurement of the flow field and flame is, therefore, obligatory for a clear understanding of the underlying mechanisms. In the current studies simultaneous PIV and OH-PLIF measurements were conducted in an enclosed gas turbine model combustor for investigating the influence of turbulence on local flame characteristics. The swirling CH₄/air flame that was investigated had a thermal power of 10.3 kW with an overall equivalence ratio of ϕ=0.75 and exhibited strong thermoacoustic oscillations at a frequency of approximately 295 Hz. The measurements reveal the formation of reaction zones at regions where hot burned gas from the recirculation zones mixes with the fresh fuel/air mixture at the nozzle exit. However, this does not seem to be a steady phenomenon as there always exist regions where the mixture has failed to ignite, possibly due to the high local strain rates present, resulting in small residence time available for a successful kinetic runaway to take place. The time averaged PIV images showed flow fields typical of enclosed swirl burners, namely a big inner recirculation zone and a small outer recirculation zone. However, the instantaneous images show the existence of small vortical structures close to the shear layers. These small vortical structures are seen playing a vital role in the formation and destruction of reaction zone structures. One does not see a smooth laminar flame front in the instantaneous OH-PLIF images, instead isolated regions of ignition and extinction highlighting the strong interplay between turbulence and chemical reactions. PACS 33.20.-t; 33.50.-j; 47.27.-i; 47.32.Ef; 47.70.Pq; 82.33.Vx; 82.40.-g