碳氢燃料超声速燃烧分区实验研究

高保真度空天发动机数值模拟通常基于快速化学反应火焰面假设,即超声速燃烧反应的特征尺度小于湍流Kolmgorov尺度,该模型方法对于氢气燃料仿真计算结果较好,但对于乙烯等碳氢燃料仍需进一步研究.受限于极端环境特种非接触测量技术,目前尚未见超声速燃烧火焰分区判别的实验研究,导致目前超声速燃烧火焰面模型适用性以及分区燃烧物理模型认识不清,进而也制约了数值发动机技术发展.本工作基于自主研发的MHz发动机内窥光纤传感器,针对单边扩张双模态冲压发动机超声速燃烧火焰分区开展实验研究,通过化学自发光信号的最小香农熵定义超声速燃烧的特征时间τsc,根据理论方法和来流工况估算了超声速燃烧的流动特征时间,结合分区燃烧理论分析了双模态超燃冲压发动机内碳氢燃料燃烧的分区情况.通过燃烧分区情况以及与泰勒尺度的比较结果,验证了碳氢燃料超燃冲压发动机典型飞行条件下燃烧室内超声速燃烧处于旋涡小火焰区域(Re?50 000; Da∈1.80～2.60, B区),多尺度湍流涡结构发挥重要作用,并随着相对于泰勒尺度的不同大小,分别对应了不同尺度的涡结构主导该过程.同时给出了当量比、通量比以及来流马赫数对燃烧特征时间的影响规律...     

高温激波管化学非平衡流动数值模拟研究

为了预测氢氧定容燃烧驱动的高温激波管性能,需要准确分析激波管非定常化学非平衡流动过程.本文在破膜前的驱动段定容燃烧以及破膜后的化学非平衡流动数值模拟中,引入双时间步长方法,发展高温激波管化学非平衡流动数值模拟方法,该方法在时间上具有二阶精度.计算结果与目前存在的激波管流动解析解以及零维化学反应系统的数值解进行了比较,吻合较好.对于典型高温激波管状态,采用有限体积方法离散准一维流动Euler控制方程,并通过将流动过程和化学反应动力学过程耦合求解,获得了激波管内部的化学非平衡流动特征.     

钝体射流火焰及其点火过程的大涡模拟

钝体燃烧器广泛应用于航空发动机、燃气轮机、锅炉等设备的燃烧室中.对其点火过程的了解和控制直接关系到设备的安全运行和污染物排放等重要问题.我们采用基于火焰面/过程变量燃烧模型的大涡模拟方法对湍流非预混钝体射流火焰及其点火过程进行了详细的数值模拟.以Sydney钝体燃烧器的无反应射流和有反应甲烷/氢气(CH_4/H_2)火焰为研究对象,首先通过统计平均的数值结果与实验测量及文献数据的对比,全面检验了所用数值方法和燃烧模型;随后,详细展示了钝体燃烧器点火和火焰发展的瞬态过程;最后,对钝体射流的点火过程进行了细致的分析和表征.根据温度峰值、羟基(OH)质量分数和甲醛(CH_2O)质量分数峰值随时间的变化表征了强制点火过程的4个阶段:点火源衰减、点火触发、点火核生成和点火成功.其中,点火核驻留的空间位置位于钝体燃烧器冷态流场外侧涡的尾部回流区域附近.     

Recent progress and challenges in fundamental combustion research（燃烧基础研究的进展和挑战）

超过80%的世界的能源转换是由燃烧方法来实现的. 发展可利用替代燃料的清洁和高效的新型发动机是解决可持续能源发展的关键之一. 在燃烧研究领域,实现这一目标的挑战是要揭示从燃料分子到发动机的多尺度燃烧过程中化学反应和火焰动力学机理,发展高效,定量的数值模拟方法和开发新的燃烧技术. 本文从7个方面综述最近几年燃烧领域的基础燃烧研究的进展和挑战. 它们包括低温清洁燃烧的发动机技术,极限条件下的燃烧机理和现象,替代燃料和混合燃料模型,多尺度化学反应模拟方法,高压燃烧反应动力学,基础燃烧的实验方法,和先进测量技术. 本文首先介绍均值充量压缩点火（HCCI）,反应控制压缩点火（RCCI）以及增压燃烧等新型发动机的概念,评述燃料特性和低温燃烧反应过程对湍流燃烧和发动机的影响,讨论发展基础燃烧研究的必要性. 第二,综述燃料浓度分层燃烧,稀薄燃烧,冷炎燃烧,以及等离子体助燃等极限燃烧条件下的新的燃烧现象和火焰机制. 第三,以航空煤油和生物柴油为例来讨论建立模拟真实燃料和替代燃料的混合燃料模型的方法. 介绍活性基指数和输运加权的反应焓的概念并用来比较燃料的高温反应特性和评价燃料的分子结构对燃烧特性的影响. 第四,评述详细化学反应机理简化的方法. 介绍多时间尺度（MTS）的化学反应的模拟和动态关联性自适应机理简化（CO-DAC）的方法来提高详细化学反应机理的计算效率. 第五,讨论高压燃烧的火焰传播速度的实验测量结果以及高压燃烧化学反应机理所存在的问题,并分析高压燃烧的关键组分和反应路径. 第六,评述测量火焰速度和组分等基础燃烧实验方法和模型中的问题和误差来源. 介绍一些改进测量方法和提高测量精度的方法. 最后,介绍测量低温燃烧中的关键组分和自由基的测量方法和最新进展.      

部分预混超声速燃烧火焰面模式研究综述

采用部分预混火焰面模式建模与仿真已经成为超声速燃烧数值研究的主要趋势之一.系统回顾了超声速燃烧火焰面模式的发展历程,针对其应用的两个基本问题进行总结:一是火焰面模式在超声速燃烧中的物理存在问题,二是超声速燃烧火焰面模式的建模问题.重点分析了火焰面模式应用于超声速燃烧的难点,提出了超声速燃烧火焰面模式建模应兼顾的问题.     

激波诱导火焰变形、混合和燃烧的数值研究

为了深入研究激波诱导的火焰变形以及由此带来的混合和燃烧变化特性,采用带单步化学反应的Navier-Stokes方程和高网格分辨率,对平面入射激波及其反射激波诱导球形火焰变形的现象进行了二维数值研究,计算结果与实验结果较好吻合。研究结果显示,在反射激波作用前,火焰的变形和皱褶主要受入射激波诱导等物理过程影响;而在反射激波与失稳火焰作用后,燃烧放热率、火焰有效面积和界面长度均迅速增加,控制火焰变形的机制逐渐向化学反应(燃烧)过程过渡;在失稳火焰发展的后期,增强的燃烧过程能够削弱火焰界面的皱褶,进而抑制未燃气和可燃气的混合。由此可以得出结论,激波诱导的火焰界面通过变形可促进界面两侧未燃气与可燃气的混合,进而强化燃烧过程,但燃烧的增强却反过来能抑制混合;认识两者之间的关系有助于利用或控制激波 火焰相互作用过程。     

湍流燃烧中的概率密度函数方法

用求解速度和化学热力学参数联合概率密度函数(pdf)输运方程的方法,计算湍流燃烧问题时,湍流输运和化学反应等过程可以精确计算,无须模拟。它还可以提供比统计矩模型方法更多的信息,因此是一个很有潜力的方法。本文给出了湍流燃烧概率密度函数的输运方程,扼要地介绍了目前应用较多的随机过程模型,以及求解概率密度函数方程的Monte Carlo算法。最后引用两个例子说明概率密度函数方法的优越性。      

煤油简化动力学和ISAT在超燃计算中的应用

利用基于``准稳态'方法建立的模型简化软件包SPARCK, 从建立的详细模型出发得到了一个包含22组分18步总包反应的煤油简化动力学模型. 简化模型计算得到的点火延迟时间与文献计算结果和实验结果相一致, 验证了模型的有效性. 采用简化模型和当地自适应建表(ISAT)方法, 对超燃冲压发动机进行了二维并行数值模拟, 计算得到的壁面压力分布与试验结果吻合较好, 表明简化模型能够很好地用来模拟煤油燃料超燃发动机内部的复杂燃烧过程. 在并行计算环境下, 和直接积分方法相比, 该方法将化学反应项的计算速度提高了3.73倍, 大大提高了计算效率.      

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性；热通量辐射模型估算辐射换热，分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解，液相采用Lagrange处理，并用PSIC算法对其进行求解．通过实验结果和计算结果的比较，表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟，3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程，而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

基于动态分区概念的高超声速燃烧大涡模拟

本文基于动态分区概念开展了亿级网格的高马赫数全尺寸超燃冲压发动机内外流耦合一体化改进延迟分离涡(IDDES)模拟研究. 研究建立了包括动态分区火焰面湍流燃烧模型(DZFM)、分区自适应化学(Z-DAC)和分区并行自适应建表(Z-ISAT)的完整动态分区燃烧模拟框架, 并通过1.15亿网格的马赫数12 REST标准高超声速燃烧室模型初步验证了分区模拟框架的保真性. DZFM通过分区解耦的思想既准确表征了当地湍流化学交互作用关系, 又有效提升了整场湍流燃烧的计算效率. Z-DAC和Z-ISAT通过在分区框架内对化学反应机理进行动态实时简化和建表查询, 可进一步提升当前分区内化学反应的求解效率. 基于1.25和1.4亿网格动态分区框架对比分析了马赫数10条件下中心支板(strut)和壁面撑挡型(pylon)两类构型氢气高超声速燃烧室特性. 支板或撑挡结构均诱发了明显的边界层分离和头部回流区, 由此两种燃烧室均出现了较长区域的喷注点前部燃烧现象. 基于Borghi图的数值分析表明当前氢气高超声速燃烧室中广泛存在扩散控制为主的火焰面模式, 效率提升的瓶颈在于高效增混. 壁面撑挡燃烧室具有较高的穿透深度和近场混合效率, 因而燃烧效率高于净推力准则80%, 相应的比冲1234 s也远高于中心支板燃烧室的437 s. 分区自适应化学方法在将近一半的计算域上降低了反应求解计算代价, 特别是在无燃料区反应机理的简化幅度更加明显. 相比与传统的有限速率PaSR模型, DZFM模型实现了高达11倍的加速比.      

Large Eddy Simulation of a Turbulent Lifted Flame using Conditional Moment Closure and Rate-Controlled Constrained Equilibrium

The objective of this paper is the simulation of a turbulent flame by employing the Rate-Controlled Constrained Equilibrium (RCCE) approach for the chemistry reduction, and Large-Eddy Simulation (LES) coupled with Conditional Moment Closure (CMC) for the turbulence-chemistry interaction modelling. RCCE is a systematic method for mechanism reduction, based on the concept that certain species characterized by faster time scales are in a constrained equilibrium state, determined by the concentration of the species controlled by the chemical kinetics. A general system of differential equations can be derived, independent on the selection of the fast and slow species (which appears as a parameter). The RCCE system is used to compute the conditional source term in the CMC equation. The flame simulated here is a methane flame issuing into a vitiated co-flow formed by hot combustion products, the ??Cabra?? flame, which is controlled by auto-ignition and is therefore sensitive to the chemical mechanism. The results show an influence of the chosen chemistry in the ignition length.     

Rapid detonation initiation by sparks in a short duct: a numerical study

Rapid onset of detonation can efficiently increase the working frequency of a pulse detonation engine (PDE). In the present study, computations of detonation initiation in a duct are conducted to investigate the mechanisms of detonation initiation. The governing equations are the Euler equations and the chemical kinetic model consists of 19 elementary reactions and nine species. Different techniques of initiation have been studied for the purpose of accelerating detonation onset with a relatively weak ignition energy. It is found that detonation ignition induced by means of multiple sparks is applicable to auto-ignition for a PDE. The interaction among shock waves, flame fronts and the strip of pre-compressed fresh (unburned) mixture plays an important role in rapid onset of detonation.     

Radical species in the cool-flame regime of diesel combustion: a comparative numerical and experimental study

Multidimensional modelling and experimental measurements are performed to study the early stages of diesel combustion. Numerical simulation is realised by means of a customised version of the KIVA 3 code, including the Shell model for auto-ignition. Experimentally, a spectroscopic analysis of the burning mixture is carried out under real operating conditions on a diesel engine equipped with an optically accessible combustion chamber. Changing the fuel injection law makes for auto-ignition to occur in environments characterised by different values of mixture pressure and temperature. Dependence of the ignition delay time upon this last variable is shown to follow a law with a negative temperature coefficient in the middle range of values. By means of natural chemiluminescence spectra, OH, CH and HCO radicals are detected as products of the reactions of thermal decomposition of the hydrocarbon molecules preceding auto-ignition. Distribution of the radicals emission intensity within the combustion chamber permits the localisation of auto-ignition sites. These are found to be in good agreement with the points of high energetic chemical activity, individuated numerically, under all the considered operating conditions. Experimentally identified radicals and fictitious species entering the reduced kinetic scheme employed within the numerical simulation are shown to exhibit an analogous behaviour regarding the trend with respect to time of the total amount of concentration, and, in a spatial sense, their distribution within the combustion chamber at the time of auto-ignition.     

Flame Stabilization Mechanisms in Lifted Flames

Flame stabilization and the mechanisms that govern the dynamics at the flame base of lifted flames have been subject to numerous studies in recent years. A combined Large Eddy Simulation-Conditional Moment Closure (LES-CMC) approach has been successful in predicting flame ignition and stabilization by auto-ignition, but accurate modelling of the competition between turbulent quenching and laminar and turbulent flame propagation at the anchor point had not been demonstrated. This paper will consolidate LES-CMC results by analysing a wide range of lifted flame geometries with different prevailing stabilization mechanisms. The simulations allow a clear distinction of these mechanisms. It is corroborated that LES-CMC accurately predicts the competition between turbulence and chemistry during the auto-ignition process, the dynamics of turbulent flame propagation can be captured, however, the dynamics of the extinction process are not approximated well under certain conditions. The averaging process inherent in the CMC methods does not allow for an instant response of the transported conditionally averaged reactive species to the changes in the flow conditions and any response of the scalars will therefore be delayed. The dimensionality of the CMC implementation affects the solution and higher dimensionality does no necessarily improve results. Stationary or quasi-stationary conditions, however, can be well predicted for all flame configurations.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.     

低温环境下甲烷爆炸流场特性模拟

在低浓度煤层气含氧液化工艺过程中,甲烷浓度会处于爆炸极限范围内,存在爆炸危险。采用流场模拟平台,对密闭容器内低温环境条件下的甲烷爆炸过程进行了数值模拟。通过研究得出:在反应体系体积及初始环境压力不变的情况下,环境温度越低,最大爆炸压力越大,到达最大爆炸压力所需时间越长;爆炸流场以化学反应区为阵面分别建立正负流动区,并不断向壁面推进,火焰传播过程受化学反应区正反馈机制的影响,在密闭容器内出现点火、加速传播、衰减传播和猝灭4个阶段;随着环境温度的降低,火焰传播速度明显降低,火焰持续时间延长。该结论可为认清低温条件下的甲烷爆炸机理及预防低浓度煤层气含氧液化工艺爆炸事故提供依据。     

圆球诱发斜爆轰波的数值研究

斜爆轰发动机是飞行器在高马赫数飞行条件下的一种新型发动机,具有结构简单、成本低和比冲高等优点.但是斜爆轰发动机的来流马赫数范围广,来流条件复杂,为实现斜爆轰波的迅速、可靠引发,采用钝头体来诱发.利用Euler方程和氢氧基元反应模型,对超声速氢气/空气混合气体中圆球诱导的斜爆轰流场进行了数值研究.不同于楔面诱发的斜爆轰波,球体首先会在驻点附近诱发正激波/爆轰波,然后在稀疏波作用下发展为斜激波/爆轰波.模拟结果显示,经过钝头体压缩的预混气体达到自燃温度后,会出现两种流场:当马赫数较低时,由于稀疏波的影响,燃烧熄灭,钝头体下游不会出现燃烧情况;而当马赫数较高时,燃烧阵面能传到下游.分析表明,当钝头体的尺度较小时,驻点附近的能量不足以诱发爆轰波,只会形成明显的燃烧带与激波非耦合结构;当钝头体的尺度较大时,流场中不会出现燃烧带与激波的非耦合现象,且这一特征与马赫数无关.通过调整球体直径,获得了激波和燃烧带部分耦合的燃烧流场结构,这一流场结构在楔面诱发的斜爆轰波中并不存在,说明稀疏波与爆轰波面的相互作用是决定圆球诱发斜爆轰波的关键.     

Numerical Study of n-Heptane Auto-ignition Using LES-PDF Methods

This work deals with pre-vaporized n-heptane auto-ignition in turbulent flows. The paper applies the Eulerian Stochastic Field method to the solution of the sub-grid joint-scalar probability density function (PDF) of the reacting scalars in the Large Eddy Simulations (LES) context. A reduced mechanism of 22 species and 18 reactions is used to model n-heptane chemical kinetics. The method is able to reproduce the different regimes observed experimentally and studies the influence of inflow temperature fluctuations.     

The Structure of the Auto-Ignition Region of Turbulent Dilute Methanol Sprays Issuing in a Vitiated Co-flow

This paper presents planar imaging of laser induced fluorescence (LIF) from key reactive species in the auto-ignition region of dilute turbulent spray flames of methanol. High-speed (5?kHz) LIF-OH imaging as well as low speed (10?Hz) imaging of joint LIF-OH-CH₂O is performed. The product of the OH and CH₂O signals is used as a qualitative indicator of local heat release. The burner is kept intentionally simple to facilitate computations and the spray is formed upstream of the jet exit plane and carried with air or nitrogen into a hot co-flowing stream of vitiated combustion products. The studied flames are all lifted but differ in the shape of their leading edge and heat release zones. Similarities with auto-ignition of gaseous fuels, as well as differences, are noted here. Formaldehyde is detected earlier than OH implying that the former is a key precursor in the initiation of auto-ignition. Growing kernels of OH that are advected from upstream, close in on the jet centreline and ignite the main flame. The existence of double reaction zones in some flames may be due to ignitable mixtures formed subsequent to local evaporation of droplets and subsequent mixing. When air is used as spray carrier, reaction zones broaden with distance, possibly due to increased partial premixing and regions of intense heat release occur near the flame centreline further downstream. With nitrogen as carrier, the flame maintains a nominal diffusion-like structure with reaction zones of uniform width and substantially less concentration of heat release on the flame centreline.     

Efficient Thermo-Chemistry Tabulation for Non-Premixed Combustion at High-Pressure Conditions

Propellant injection and turbulent combustion in high-pressure engines is often dominated by real-gas effects. However, previous studies suggested that the departure of the fluid properties from an ideal gas behavior has only a limited effect on the laminar flame structure. This is due to the fact that chemical reactions take place in the flame zone where the temperature is sufficiently high and molecular interactions are negligible, i.e., the ideal gas assumption is valid. On the other hand, various experimental and numerical studies of injection processes at high-pressure conditions demonstrated that real-gas effects can have a strong impact on the turbulent flow. Mixing is influenced by the rapid change of fluid properties. In this work, we exploit the gap in the fidelity of the thermodynamics model needed to describe the laminar flame structure and that needed to describe the turbulent flow field. We then propose a new real-gas flamelet model with increased numerical performance. The computational cost of the new formulation is not significantly higher than that of an ideal gas simulation. The performance of the method is analyzed and the error that is introduced by our assumptions is assessed by comparison to more complete modeling. Finally, the method is used to simulate a turbulent jet flame emanating from a coaxial injector at supercritical pressure and cryogenic oxidizer temperature. The results are compared with experimental OH^? images giving evidence of the suitability of the present method.