一维湍流模型对He平面羽流和CH4/H2/N2射流火焰的数值模拟

本文采用一维湍流模型（ODT）对氦气平面羽流和CH₄/H₂/N₂射流火焰进行数值模拟,和前人的实验结果进行定量地对比。结果表明,ODT模型能够准确地预测平面羽流基本特征,湍流涡团的分布同流场拉伸率之间具有密切的关系,涡团强度的分布能够直观地表明当地的湍流强度。ODT模型埘CH₄/H₂/N₂瞬态火焰的模拟定性反应了火焰特性及其与湍流作用的规律,对温度-混合分数的预测值和实测值进行比较,发现甲烷火焰燃烧在富燃料侧并未达到平衡状态,因而基元反应对火焰特性的预测具有重要作用。     

稀释剂对H2/CH4预混火焰熄灭特性影响

掺氢天然气在稀释气体作用下的熄灭特性研究对实际燃烧设备的设计和优化具有重要的指导意义。本文利用对冲火焰法测量了掺氢天然气层流火焰在N₂和CO₂作用下的熄灭拉伸率,并采用数值模拟耦合详细化学反应机理对N₂,CO₂和He的稀释剂效应展开研究。结果表明,Li、GRI Mech 3.0和FFCM-1机理均能定性反映燃料熄灭拉伸率随当量比的变化规律,且FFCM-1机理综合预测精度最高。实验和模拟发现,不同稀释剂气体对掺氢天然气熄灭拉伸率降低幅度满足:He22。进一步研究发现,CO₂由于热容大,在反应体系中会降低火焰温度,同时增强了链终止反应强度,通过热效应和化学效应两方面对火焰熄灭特性起作用。He则能显著改变燃料混合物的平均摩尔质量,从而改变体系中重要反应物和自由基的扩散特性,从扩散效应方面影响火焰的熄灭特性。     

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

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

H2/air旋流预混燃烧的流场特性和火焰结构分析

本文采用完全可压缩的N-S方程,对当量比为1.0的H₂/air旋流预混火焰进行了直接数值模拟研究。氢气和空气的化学反应采用9种组分19步的详细机理。模拟结果表明,强旋流流场中存在回流区,碗形旋流火焰稳定在回流区的外围。在火焰面上沿火焰法向提取了局部火焰结构,将局部湍流火焰结构与层流预混火焰的火焰结构进行了比较,发现局部湍流火焰比层流预混火焰更薄,燃烧强度更高。     

CH4/空气简化动力学机理数值验证

选择绕圆柱预混燃烧算例，验证CH₄/空气三种简化动力学机理（16s41r、15s19r和53s325r）.考虑均匀来流，忽略湍流和湍流与燃烧相互作用以及燃料扩散效应，假设层流有限反应速率，采用保自由流5阶WENO格式求解多组分Euler方程组，得到CH₄/空气预混燃烧流场温度等值线、沿驻点线压力和温度及其CH₄、CO和CO₂质量百分数分布.结果表明：三种简化动力学机理给出的流场均出现弓形激波和火焰面，弓形激波和火焰驻点距离及其形状、诱导区宽度和简化动力学机理相关.当圆柱直径增大，弓形激波和火焰向圆柱上游移动，对应的驻点距离均增大，诱导区宽度变短，点火延时变小，但火焰和弓形激波位置次序未变化.53s325r模型要比16s41r模型和15s19r模型精度要高，点火延时覆盖的压力和温度范围也较宽，所有简化机理均未完全反应，在较大圆柱直径下游达到化学平衡.     

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

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

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

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

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

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

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

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

CS2燃烧火焰光谱辐射模型的构建与特征污染产物浓度反演研究

CS₂在当今化工等领域占据了重要地位,而CS₂火灾污染事故危害性极大。通过研究CS₂燃烧火焰光谱辐射以探究其火灾污染特性极为必要。搭建了CS₂燃烧火焰光谱测试平台,采用黑体辐射源对VSR仪器进行了标定,通过多用途傅里叶变换(VSR)红外光谱辐射仪测试了5,10和20 cm三种燃烧尺度下CS₂燃烧的火焰光谱,并通过热电偶测试了整个燃烧时间段内不同燃烧时刻下的火焰温度,以及在火焰上方安装了烟气分析仪对火焰中的燃烧产物浓度进行监测。测量了CS₂整个燃烧时间段内火焰温度,以及不同燃烧时间、不同燃烧尺度下的火焰光谱、燃烧产物组分信息。测试结果表明,CS₂火焰中主要含有高温SO₂,CO₂,CO气体和空气中卷入的H₂O分子,并获取了特征污染产物SO₂的浓度。由于现有光谱仪测量分辨率有限,室内实验测量的火焰尺度有限,为了能实现火灾在线监测需要建立一个火焰光谱辐射模型来反演CS₂火灾时的污染物浓度相关信息。基于HITRAN数据库可知在2.7 μm附近为高温水蒸气的发射峰,4.2 μm附近特征峰为高温CO₂气体的发射峰,4.7 μm附近有CO微弱的发射峰,在7.4 μm附近特征峰为高温SO₂气体的发射峰,并获得了CS₂燃烧时产生的SO₂,CO₂,CO和H₂O气体在火焰燃烧相同温度下的吸收系数,通过计算得到了CS₂燃烧时产生的SO₂,CO₂,CO和H₂O混合气体的透过率与发射率,并结合气体辐射传输方程、气体吸收系数等方程,创建了CS₂燃烧的火焰光谱辐射模型。利用该光谱辐射模型反演了不同燃烧时间下特征污染产物SO₂的浓度,并与实验测得的数据进行了对比分析。结果表明,该模型精度高,可用于燃烧产物浓度的定量化反演,SO₂分子含量在燃烧时间20,40,60和80 s时的反演精度分别是89.5%,82.5%,85.6%和86.5%。为遥感反演CS₂型大尺度火灾中燃烧产物的浓度奠定基础。     

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.     

Application and theoretical analysis of the flamelet model for supersonic turbulent combustion flows in the scramjet engine

In the framework of Reynolds-averaged Navier–Stokes simulation, supersonic turbulent combustion flows at the German Aerospace Centre (DLR) combustor and Japan Aerospace Exploration Agency (JAXA) integrated scramjet engine are numerically simulated using the flamelet model. Based on the DLR combustor case, theoretical analysis and numerical experiments conclude that: the finite rate model only implicitly considers the large-scale turbulent effect and, due to the lack of the small-scale non-equilibrium effect, it would overshoot the peak temperature compared to the flamelet model in general. Furthermore, high-Mach-number compressibility affects the flamelet model mainly through two ways: the spatial pressure variation and the static enthalpy variation due to the kinetic energy. In the flamelet library, the mass fractions of the intermediate species, e.g. OH, are more sensible to the above two effects than the main species such as H₂O. Additionally, in the combustion flowfield where the pressure is larger than the value adopted in the generation of the flamelet library or the conversion from the static enthalpy to the kinetic energy occurs, the temperature obtained by the flamelet model without taking compressibility effects into account would be undershot, and vice versa. The static enthalpy variation effect has only little influence on the temperature simulation of the flamelet model, while the effect of the spatial pressure variation may cause relatively large errors. From the JAXA case, it is found that the flamelet model cannot in general be used for an integrated scramjet engine. The existence of the inlet together with the transverse injection scheme could cause large spatial variations of pressure, so the pressure value adopted for the generation of a flamelet library should be fine-tuned according to a pre-simulation of pure mixing.     

Large eddy simulation of Cambridge bluff-body coal (CCB2) flames with a flamelet progress variable model

In the present study, we report the first large eddy simulation (LES) study of the Cambridge CCB2 coal flames, one of the target flames in the Workshop on Measurement and Simulation of Coal and Biomass Conversion, with an extended flamelet progress variable (FPV) model. The extended FPV model is based on two mixture fractions considering the volatiles and char off-gases. The normalized total enthalpy is used for the interphase heat transfer modelling. Turbulence-chemistry interaction is treated with an assumed probability density function approach. The results show that the present LES can generally capture the flow field and particle distribution, while there are considerable deviations in the OH prediction due to the boundary treatment of using a mixture of volatiles and carrier gas to replace the methane-containing mixtures in the primary and pilot flow. It indicates that for such gas-assisted coal flames, the pilot fuel stream needs to be rigorously considered in the flamelet tabulation that could be resolved by extending the two-mixture-fraction model into a three-mixture-fraction model. The instantaneous Lagrangian particles histories show that the increasing of coal load has a negligible effect on devolatilization, but delays the char conversion.     

带控制的超声速二维后台阶流动显示

在Mach数3.4的来流条件下,对二维后台阶流动精细结构开展了实验研究.实验分为后台阶上游无控制加粗糙带扰动及微涡流发生器（micro-vortex generator,MVG）扰动3种状态,采用基于纳米示踪的平面激光散射（nano-tracer based planar laser scattering,NPLS）方法获得了流向和展向切面内的高时空分辨率流动显示图像,并测量了模型表面静压分布.对大量NPLS图像取平均,研究了流场结构的时间平均规律,对比不同时刻的瞬态流场精细结构图像,发现不同状态下的湍流大尺度结构的特征时间.有粗糙带状态相对无粗糙带台阶下游回流区压力更低,而下游压力较高,台阶上游区别不大;受MVG控制后台阶下游附近区域压力突增;MVG对流动的控制改变能力较强,粗糙带能调整台阶上下游附近流动平稳过渡,流场壁面压力没有突变.      

Scalar dissipation rates in isothermal and reactive turbulent opposed-jets: 1-D-Raman/Rayleigh experiments supported by LES

Isothermal and reactive turbulent opposed flows are presented, which are appropriate to test the applicability and performance of models for turbulence, mixing, chemical reaction, and turbulence-chemistry interaction. Transient flow and scalar fields are measured using laser Doppler velocimetry and one-dimensionally resolved Raman/Rayleigh spectroscopy. Aside of statistical moments of temperature, mean species, and velocity components, scalar dissipation rate across the mixing and reaction layer is determined on a single-shot base. Using large eddy simulation in connection with a steady flamelet model, it is shown how numerical data can serve to estimate the influence of experimental noise upon a measured quantity, such as scalar dissipation. As a key result, it is shown that an increase in scalar rate of dissipation by chemical reactions is caused by a significant increase in the mixture fraction diffusivity, which outweighs the decrease in mixture fraction gradients. In mixture fraction space, local maxima of scalar dissipation rate are found on the rich side, which cannot be correctly reproduced by the steady flamelet model assuming equal species diffusivity. Furthermore, the impact of experimental noise on conditional probability density functions of scalar dissipation rate is shown (exemplary) to lead to erroneous conclusions from experimental data.     

Effects of Lewis number,density ratio and gravity on burning velocity and conditional statistics in stagnating turbulent premixed flames

DNS is performed to analyse the effects of Lewis number (Le), density ratio and gravity in stagnating turbulent premixed flames. The results show good agreement with those of Lee and Huh (Combustion and Flame, Vol. 159, 2012, pp. 1576–1591) with respect to the turbulent burning velocity, S_T, in terms of turbulent diffusivity, flamelet thickness, mean curvature and displacement speed at the leading edge. In all four stagnating flames studied, a mean tangential strain rate resulting in a mean flamelet thickness smaller than the unstretched laminar flame thickness leads to an increase in S_T. A flame cusp of positive curvature involves a superadiabatic burned gas temperature due to diffusive–thermal instability for an Le less than unity. Wrinkling tends to be suppressed at a larger density ratio, not enhanced by hydrodynamic instability, in the stagnating flow configuration. Turbulence is produced, resulting in highly anisotropic turbulence with heavier unburned gas accelerating through a flame brush by Rayleigh–Taylor instability. Results are also provided on brush thickness, flame surface density and conditional velocities in burned and unburned gas and on flame surfaces to represent the internal brush structures for all four test flames.