One-dimensional turbulence model simulations of autoignition of hydrogen/carbon monoxide fuel mixtures in a turbulent jet

The autoignition of hydrogen/carbon monoxide in a turbulent jet with preheated co-flow air is studied using the one-dimensional turbulence (ODT) model. The simulations are performed at atmospheric pressure based on varying the jet Reynolds number and the oxidizer preheat temperature for two compositions corresponding to varying the ratios of H₂ and CO in the fuel stream. Moreover, simulations for homogeneous autoignition are implemented for similar mixture conditions for comparison with the turbulent jet results. The results identify the key effects of differential diffusion and turbulence on the onset and eventual progress of autoignition in the turbulent jets. The differential diffusion of hydrogen fuels results in a reduction of the ignition delay relative to similar conditions of homogeneous autoignition. Turbulence may play an important role in delaying ignition at high-turbulence conditions, a process countered by the differential diffusion of hydrogen relative to carbon monoxide; however, when ignition is established, turbulence enhances the overall rates of combustion of the non-premixed flame downstream of the ignition point.     

Understanding ignition processes in spray-guided gasoline engines using high-speed imaging and the extended spark-ignition model SparkCIMM: Part B: Importance of molecular fuel properties in early flame front propagation

Recent high-speed imaging of ignition processes in spray-guided gasoline engines has motivated the development of the physically-based spark channel ignition monitoring model SparkCIMM, which bridges the gap between a detailed spray/vaporization model and a model for fully developed turbulent flame front propagation. Previously, both SparkCIMM and high-speed optical imaging data have shown that, in spray-guided engines, the spark plasma channel is stretched and wrinkled by the local turbulence, excessive stretching results in spark re-strikes, large variations occur in turbulence intensity and local equivalence ratio along the spark channel, and ignition occurs in localized regions along the spark channel (based upon a Karlovitz-number criteria).In this paper, SparkCIMM is enhanced by: (1) an extended flamelet model to predict localized ignition spots along the spark plasma channel, (2) a detailed chemical mechanism for gasoline surrogate oxidation, and (3) a formulation of early flame kernel propagation based on the G-equation theory that includes detailed chemistry and a local enthalpy flamelet model to consider turbulent enthalpy fluctuations. In agreement with new experimental data from broadband spark and hot soot luminosity imaging, the model establishes that ignition prefers to occur in fuel-rich regions along the spark channel. In this highly-turbulent highly-stratified environment, these ignition spots burn as quasi-laminar flame kernels. In this paper, the laminar burning velocities and flame thicknesses of these kernels are calculated along the mean turbulent flame front, using tabulated detailed chemistry flamelets over a wide range of stoichiometry and exhaust gas dilution. The criteria for flame propagation include chemical (cross-over temperature based) and turbulence (Karlovitz-number based) effects. Numerical simulations using ignition models of different physical complexity demonstrate the significance of turbulent mixture fraction and enthalpy fluctuations in the prediction of early flame front propagation. A third paper on SparkCIMM (companion paper to this one) focuses on the importance of molecular fuel properties and flame curvature on early flame propagation and compares computed flame propagation with high speed combustion imaging and computed heat release rates with cylinder pressure analysis.The goals of SparkCIMM development are to (a) enhance our fundamental understanding of ignition and combustion processes in highly-turbulent highly-stratified engine conditions, (b) incorporate that understanding into a physically-based submodel for RANS engine calculations that can be reliably used without modification for a wide range of conditions (i.e., homogeneous or stratified, low or high turbulence, low or high dilution), and (c) provide a submodel that can be incorporated into a future LES model for physically-based modeling of cycle-to-cycle variability in engines.     

一种可变湍流参数的定容燃烧弹试验装置

设计了一种可变湍流参数定容燃烧弹试验装置,用于预混湍流燃烧的基础研究。该试验装置利用带有通孔的板在容弹内的快速平动来产生可控的湍流,它能方便地改变容弹内充量的湍流参数（包括湍流强度和湍流标尺）、热力参数（包括燃空比、残余废气系数和压力）以及点火参数（火花塞位置、电极间隙与点火能量）。介绍了该试验装置的定容燃烧弹、湍流发生系统、混合气配制系统、点火系统、燃烧压力测量系统、纹影与高速摄影系统以及时序控制系统等子系统的设计原理、结果及其特点。     

用紊流卷吸燃烧模型计算研究点火时刻紊流参数对燃烧循环的影响

利用准维紊流卷吸燃烧模型计算研究了点火时刻气缸内紊流参数的变化对汽油机燃烧循环变动的影响。结果表明，点火时刻紊流强度u′及紊流长度积分标尺L的变动对汽油机燃烧循环变动均有很大影响，增大u′或L均有利于加快燃烧速度及火焰传播速度，缩短火焰发展期；u′或L的变动加大，燃烧循环变动也随之加大，其中，u′对燃烧循环变动的影响尤其大。     

初始进气充量成分对末端气体自燃化学特性的影响

本文介绍了在一台单缸二冲程爆震实验发动机上，研究进气充量中不同的残余气体成分对焰前反应及自燃的影响的实验结果。结果表明，初始进气充量中含有一定量的残余废气时，抑制了自燃；在残余废气中含有一定量的部分氧化产物时，促进了低温化学反应和自燃。低温化学反应所释放的热量以及活性中间产物的积累导致了终燃混合气的自燃。ＲＯＮ５３汽油的自燃过程表现为两阶段着火特征。     

Advancing lean combustion of hydrogen-air mixtures by laser-induced spark ignition

We investigate how ignition through laser-induced plasma can improve the application of lean combustion, in particular in environmental conditions relevant to hydrogen internal combustion engines (ICE). Major design goals when developing combustion engines are increasing thermal efficiency and decreasing combustion emissions. High compression ratios, lean combustion and precise ignition timings are contributing factors in ICE optimization. In our studies, several gains from laser spark ignition are investigated. The high energy content of laser-induced ignition kernels are shown to speed up the development of the early flame kernels. These extended ignition kernels transfer into self propagating flames even in lean mixtures. Leaner mixtures are ignited in our experiments using laser spark ignition in comparison to conventional electrical spark plugs. Precise ignition timing is realized. Multi-point ignitions are synchronized on the timescale of microseconds to enhance the progress of combustion. We modified the locus of ignition in a mixture flow to decrease the temporal extent of flame contact with the wall. Therefore, burning duration and heat loss can be reduced.     

Hydrogen autoignition in a turbulent jet with preheated co-flow air

The autoignition of hydrogen in a turbulent jet with preheated air is studied computationally using the stand-alone one-dimensional turbulence (ODT) model. The simulations are based on varying the jet Reynolds number and the mixture pressure. Also, computations are carried out for homogeneous autoignition at different mixture fractions and the same two pressure conditions considered for the jet simulations. The simulations show that autoignition is delayed in the jet configuration relative to the earliest autoignition events in homogeneous mixtures. This delay is primarily due to the presence of scalar dissipation associated with the scalar mixing layer in the jet configuration as well as with the presence of turbulent stirring. Turbulence plays additional roles in the subsequent stages of the autoignition process. Pressure effects also are present during the autoignition process and the subsequent high-temperature combustion stages. These effects may be attributed primarily to the sensitivity of the autoignition delay time to the mixture conditions and the role of pressure and air preheating on molecular transport properties. The overall trends are such that turbulence increases autoignition delay times and accordingly the ignition length and pressure further contribute to this delay.     

火花点火发动机准维湍流卷吸燃烧模型的适用性研究（1）

本文在分析火花点火发动机湍流涡结构及缸内湍流特性参数的基础上，提出了适用于火花点火发动机燃烧计算的准维湍流卷吸模型，通过建立相应的子模型及求解方程，实现了燃烧过程的计算；对压缩比为１０的紧凑型燃烧室，在改变发动机转速、负荷、空燃比以及点火正时的情况下，计算得到的压力示功图、质量燃烧率等与实测值一致，从而证实了该模型的合理性。     

The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels.     

DNS of spark ignition and edge flame propagation in turbulent droplet-laden mixing layers

A parametric study of forced ignition at the mixing layer between air and air carrying fine monosized fuel droplets is done through one-step chemistry direct numerical simulations to determine the influence of the size and volatility of the droplets, the spark location, the droplet-air mixing layer initial thickness and the turbulence intensity on the ignition success and the subsequent flame propagation. The propagation is analyzed in terms of edge flame displacement speed, which has not been studied before for turbulent edge spray flames. Spark ignition successfully resulted in a tribrachial flame if enough fuel vapour was available at the spark location, which occurred when the local droplet number density was high. Ignition was achieved even when the spark was offset from the spray, on the air side, due to the diffusion of heat from the spark, provided droplets evaporated rapidly. Large kernels were obtained by sparking close to the spray, since fuel was more readily available. At long times after the spark, for all flames studied, the probability density function of the displacement speed was wide, with a mean value in the range 0.55-0.75S_L, with S_L the laminar burning velocity of a stoichiometric gaseous premixed flame. This value is close to the mean displacement speed in turbulent edge flames with gaseous fuel. The displacement speed was negatively correlated with curvature. The detrimental effect of curvature was attenuated with a large initial kernel and by increasing the thickness of the mixing layer. The mixing layer was thicker when evaporation was slow and the turbulence intensity higher. However, high turbulence intensity also distorted the kernel which could lead to high values of curvature. The edge flame reaction component increased when the maximum temperature coincided with the stoichiometric contour. The results are consistent with the limited available experimental evidence and provide insights into the processes associated with ignition of practical spray flames.     

Ignition of turbulent swirling n-heptane spray flames using single and multiple sparks

This paper examines ignition processes of an n-heptane spray in a flow typical of a liquid-fuelled burner. The spray is created by a hollow-cone pressure atomiser placed in the centre of a bluff body, around which swirling air induces a strong recirculation zone. Ignition was achieved by single small sparks of short duration (2 mm; 0.5 ms), located at various places inside the flow so as to identify the most ignitable regions, or larger sparks of longer duration (5 mm; 8 ms) repeated at 100 Hz, located close to the combustion chamber enclosure so as to mimic the placement and characteristics of a gas turbine combustor surface igniter. The air and droplet velocities, the droplet diameter, and the total (i.e. liquid plus vapour) equivalence ratio were measured in inert flow by phase Doppler anemometry and sampling respectively. Fast camera imaging suggested that successful ignition events were associated with flamelets that propagated back towards the spray nozzle. Measurements of ignition probability with the single spark showed that localised ignition inside the spray is more likely to result in successful flame establishment when the spark is located in a region of negative velocity, relatively small droplet Sauter mean diameter, and mean equivalence ratio within the flammability limits. Ignition with the single spark was not possible at the location where the multiple spark experiments were performed. For those, the multiple spark sequence lasted approximately 1 to 5 s. It was found that a long spark sequence increases the ignition efficiency, which reached a maximum of 100% at the axial distance where the recirculation zone had maximum width. Ignition was not feasible with the spark downstream of about two burner diameters. Visualisation showed that small flame kernels emanate very often from the spark, which can be stretched as far as 20 mm from the electrodes by the turbulent velocity fluctuations. These kernels survive very little time. Successful overall ignition occurs at a random time from the spark initiation and, as in the case of the single spark, success is associated with kernels that move without getting extinguished towards the bluff body. The results demonstrate that the energy deposited by multiple sparks and spark stretching in a turbulent flow can have a spatially far-reaching effect to initiate combustion.     

Spark ignition transitions in premixed turbulent combustion

Recent discoveries and developments on the dynamic process of premixed turbulent spark ignition are reviewed. The focus here is on the variation of turbulent minimum ignition energies (MIE_T) against laminar MIE (MIE_L) over a wide range of r.m.s. turbulence fluctuation velocity (uʹ) alongside effects of the spark gap between electrodes, Lewis number, and some other parameters on MIE. Two distinguishable spark ignition transitions are discussed. (1) A monotonic MIE transition, where MIE_L sets the lower bound, marks a critical uʹ_c between linear and exponential increase in MIE_T with uʹ increased. (2) A non-monotonic MIE transition, where the lower bound is to be set by a MIE_T at some uʹ_c, stems from a great influence of Lewis number and spark gap despite turbulence. At sufficiently large Lewis number >> 1 and small spark gap (typically less than 1 mm), turbulence facilitated ignition (TFI), where MIE_T < MIE_L, occurs; then MIE_T increases rapidly at larger uʹ > uʹ_c because turbulence re-asserts its dominating role. Both phenomena are explained by the coupling effects of differential diffusion, heat losses to electrodes, and turbulence on the spark kernel. In particular, the ratio of small-scale turbulence diffusivity to reaction zone thermal diffusivity, a reaction zone Péclet number, captures the similarity of monotonic MIE transition, regardless of different ignition sources (conventional electrodes versus laser), turbulent flows, pressure, and fuel types. Furthermore, TFI does and/or does not occur when conventional spark is replaced by nanosecond-repetitively-pulsed-discharge and/or laser spark. The latter is attributed to the third lobe formation of laser kernel with some negative curvature segments that enhance reaction rate through differential diffusion, where MIE_L < MIE_T (no TFI). Finally, the implications of MIE transitions relevant to lean-burn spark ignition engines are briefly mentioned, and future studies are suggested.     

柴油引燃天然气发动机燃烧过程的研究

对柴油引燃天然气发动机的着火过程、燃烧过程进行研究,提出了相应的数学模型,并进行了试验验证。着火模型的基本框架仍然采用Shell模型,但考虑到CH4对柴油着火过程的影响,修正了Shell模型。对CH4的均质燃烧过程,以简单的一步表观反应动力学概念的阿伦纽斯(Arrhenius)公式为基础,综合湍流脉动对化学动力学的影响,提出了一个新的燃烧模型。模拟计算借助了KIVA-3软件,修正的Shell模型和燃烧模型作为独立子程序与KIVA-3进行了耦合。研究表明,讨论的模型能够较好模拟柴油引燃天然气发动机的着火和燃烧过程。     

Simulations of spray autoignition and flame establishment with two-dimensional CMC

The unsteady two-dimensional conditional moment closure (CMC) model with first-order closure of the chemistry and supplied with standard models for the conditional convection and turbulent diffusion terms has been interfaced with a commercial engine CFD code and analyzed with two numerical methods, an “exact” calculation with the method of lines and a faster fractional-step method. The aim was to examine the sensitivity of the predictions to the operator splitting errors and to identify the extent to which spatial transport terms are important for spray autoignition problems. Despite the underlying simplifications, solution of the full CMC equations allows a single model to be used for the autoignition, flame propagation (“premixed mode”), and diffusion flame mode of diesel combustion, which makes CMC a good candidate model for practical engine calculations. It was found that (i) the conditional averages have significant spatial gradients before ignition and during the premixed mode and (ii) that the inclusion of physical-space transport affects the calculation of the autoignition delay time, both of which suggest that volume-averaged CMC approaches may be inappropriate for diesel-like problems. A balance of terms in the CMC equation before and after autoignition shows the relative magnitude of spatial transport and allows conjectures on the structure of the premixed phase of diesel combustion. Very good agreement with available experimental data is found concerning ignition delays and the effect of background air turbulence on them.     

Chemical structure of autoignition in a turbulent lifted H2/N2 jet flame issuing into a vitiated coflow

In the present paper autoignition is studied as the main stabilization mechanism in turbulent lifted H₂/N₂ jet flames issuing into a vitiated hot coflow. The numerical study is performed using the joint scalar PDF approach with detailed chemistry in a two dimensional axisymmetric domain. The SSG Reynolds stress model is used as a turbulence model in the simulation. Chemical structure and characteristics of autoignition are investigated using various methods and parameters. Reaction rate analysis is made to analyze the ignition process at the flame base. The results show the occurrence of a chain branching reaction preceding thermal runaway, which boosts the chain branching process in the flame. This demonstrates the large impact of autoignition at the flame base on the stabilization of the lifted turbulent flame. Further investigation using the scatter-plots of scalars reveals the characteristics of the ignition. The relation between the behavior of temperature and of key intermediate species demonstrates the formation of OH through consumption of HO₂ at nearly isothermal conditions in a very lean-fuel mixture at the flame base. Flux analyses in the conservation equations of species are used to explore the impacts of mass transport on ignition process. Ignition is found to be mainly controlled by chemical features rather than the mixing processes near the flame base. Characteristics of autoignition are also investigated in terms of Damköhler number and progress variable.     

Effect of the turbulence intensity on knocking tendency in a SI engine with high compression ratio using biogas and blends with natural gas,propane and hydrogen

This research presents the test results carried out in a diesel engine converted to spark ignition (SI) using gaseous fuels, applying a geometry change of the pistons combustion chamber (GCPCC) to increase the turbulence intensity during the combustion process; with similar compression ratio (CR) of the original diesel engine; the increase in turbulence intensity was planned to rise turbulent flame speed of biogas, to compensate its low laminar flame speed. The research present the test to evaluate the effect of increase turbulence intensity on knocking tendency; using fuel blends of biogas with natural gas, propane and hydrogen; for each fuel blend the maximum output power was measured just into the knocking threshold before and after GCPCC; spark timing (ST) was adjusted for optimum generating efficiency at the knocking threshold. Turbulence intensity with GCPCC was estimated using Fluent 13, with 3D Combustion Fluid Dynamics (CFD) numerical simulations; 12 combustion chamber geometries were simulated in motoring conditions; the selected geometry had the greatest simulated turbulent kinetic energy (TKE) and Reynolds number (Re) during combustion. The increased turbulence intensity was measured indirectly through the periods of combustion duration to mass fraction burn 0–5%, 0–50% and 0–90%; for almost all the fuel blends the increased turbulence intensity of the engine, increased the knocking tendency requiring to reduce the maximum output power to keep engine operation just into the knocking threshold. Biogas was the only fuel without power derating by the conditions of higher pressure and higher turbulence during combustion by GCPCC and improve its generating efficiency. Peak pressure, heat release rate, mean effective pressure and exhaust temperature were lower after GCPCC. Tests results indicated that knocking tendency was increased because of the higher turbulent flame speed; fuel blends with high laminar flame speed and low methane number (MN) had higher knocking tendency and lower output power.     

Numerical investigation of a large bore,direct injection,spark ignition,hydrogen-fuelled engine

This paper presents Unsteady Reynolds Averaged Navier Stokes (URANS) simulations of a large bore, hydrogen-fuelled direct injection spark ignition (DISI) engine at different spark and start of injection (SOI) timings. Six cases are simulated, including three with various spark timings at a low boost level and three with advanced to late injection timings at a higher boost level. The numerical simulations are validated with experimental data for four out of six cases, while the other two are considered blind computational fluid dynamics (CFD) simulations. It is shown that the autoignition occurs with advanced spark timing due to high in-cylinder pressure and unburnt temperature. For different SOIs, it is demonstrated that flame propagation involves a spark-initiated flame combined with an autoignition generated flame. The case with the late injection timing features poor mixing and slower combustion due to the presence of lean mixtures near the spark plug. As a result, this case features the lowest thermal efficiency when SOI is varied. In all cases, both mixture and temperature stratification are present. Simulations of zero-dimensional chemical reactors demonstrate that this stratification must be correctly captured for accurate prediction of autoignition timing.     

Stochastic modeling of autoignition in turbulent non-homogeneous hydrogen–air mixtures

The autoignition of turbulent non-homogeneous hydrogen–air mixtures is studied using the linear-eddy model (LEM). The initial solution consists of fully segregated regions of fuel and oxidizer mixtures. The relative size of these regions represents a measure of mixture heterogeneity, while the specified turbulence conditions determine the subsequent evolution of the dissipation rate field. Chemistry and transport are described accurately using a detailed mechanism for hydrogen–air chemistry and the CHEMKIN libraries. The simulations are implemented for a range of pressures and initial mixing conditions to identify the effects of mixing on the dominant autoignition chemistry. The simulations show that some of the salient features of the coupling between autoignition chemistry and mixing may adequately be captured by the LEM. This coupling includes the competing roles of mixing on this chemistry. Mixing can increase the volumetric rate of reaction and heat release by increasing the interface between ignition fronts and unburnt gases; it also contributes to homogenizing the mixture.     

火花点火发动机准维湍流卷吸燃烧模型的适用性研究（二）

本文利用作者提出的火花点火发动机准维湍流卷吸燃烧模型，对压缩比为１０和１２的火球形燃烧室以及压缩比为１０的碗形燃烧室变工况进行了计算，将计算得到的示功图、质量燃烧率等与实验值进行了对比对分析。结果表明，合理选取与燃烧室结构相对应的四个经验常数，准维湍流卷吸燃烧模型完全适用于火花点火发动机变工况及不同燃烧室结构工作过程的计算，能够正确反映火花点火发动机结构参数和运转参数对燃烧过程的影响。