On laminar premixed flame propagating into autoigniting mixtures under engine-relevant conditions

Usually premixed flame propagation and laminar burning velocity are studied for mixtures at normal or elevated temperatures and pressures, under which the ignition delay time of the premixture is much larger than the flame resistance time. However, in spark-ignition engines and spark-assisted compression ignition engines, the end-gas in the front of premixed flame is at the state that autoignition might happen before the mixture is consumed by the premixed flame. In this study, laminar premixed flames propagating into an autoigniting dimethyl ether/air mixture are simulated considering detailed chemistry and transport. The emphasis is on the laminar burning velocity of autoigniting mixtures under engine-relevant conditions. Two types of premixed flames are considered: one is the premixed planar flame propagating into an autoigniting DME/air without confinement; and the other is premixed spherical flame propagating inside a closed chamber, for which four stages are identified. Due to the confinement, the unburned mixture is compressed to high temperature and pressure close to or under engine-relevant conditions. The laminar burning velocity is determined from the constant-volume propagating spherical flame method as well as PREMIX. The laminar burning velocities of autoigniting DME/air mixture at different temperatures, pressures, and autoignition progresses are obtained. It is shown that the first-stage and second-stage autoignition can significantly accelerate the flame propagation and thereby greatly increase the laminar burning velocity. When the first-stage autoignition occurs in the unburned mixture, the isentropic compression assumption does not hold and thereby the traditional method cannot be used to calculate the laminar burning velocity. A modified method without using the isentropic compression assumption is proposed. It is shown to work well for autoigniting mixtures. Besides, a power law correlation is obtained based on all the laminar burning velocity data. It works well for mixtures before autoignition while improvement is still needed for mixtures after autoignition.     

Effects of fuel stratification on ignition kernel development and minimum ignition energy of n-decane/air mixtures

Fuel-stratified combustion has broad application due to its promising advantages in extension of lean flammability limit, improvement of flame stabilization, enhancement of lean combustion, etc. In the literature, there are many studies on flame propagation in fuel-stratified mixtures. However, there is little attention on ignition in fuel-stratified mixtures. In this study, one-dimensional numerical simulation is conducted to investigate the ignition and spherical flame kernel propagation in fuel-stratified n-decane/air mixtures. The emphasis is placed on assessing the effects of fuel stratification on the ignition kernel propagation and critical ignition condition. First, ignition and flame kernel propagation in homogeneous n-decane/air mixture are studied and different flame regimes are identified. The minimum ignition energy (MIE) of the homogeneous n-decane/air mixture is obtained and it is found to be very sensitive to the equivalence ratio under fuel-lean conditions. Then, ignition and flame kernel propagation in fuel-stratified n-decane/air mixture are investigated. The inner equivalence ratio and stratification radius are found to have great impact on ignition kernel propagation. The MIEs at different fuel-stratification conditions are calculated. The results indicate that for fuel-lean n-decane/air mixture, fuel stratification can greatly promote ignition and reduce the MIE. Six distinct flame regimes are observed for successful ignition in fuel-stratified mixture. It is shown that the ignition kernel propagation can be induced by not only the ignition energy deposition but also the fuel-stratification. Moreover, it is found that to achieve effective ignition enhancement though fuel stratification, one needs properly choose the values of stratification radius and inner equivalence ratio.     

LES/CMC modelling of ignition and flame propagation in a non-premixed methane jet

The Large Eddy Simulation (LES) / Conditional Moment Closure (CMC) model with detailed chemistry is used for modelling spark ignition and flame propagation in a turbulent methane jet in ambient air. Two centerline and one off-axis ignition locations are simulated. We focus on predicting the flame kernel formation, flame edge propagation and stabilization. The current LES/CMC computations capture the three stages reasonably well compared to available experimental data. Regarding the formation of flame kernel, it is found that the convection dominates the propagation of its downstream edge. The simulated initial downstream and radial flame propagation compare well with OH-PLIF images from the experiment. Additionally, when the spark is deposited at off-centerline locations, the flame first propagates downstream and then back upstream from the other side of the stoichiometric iso-surface. At the leading edge location, the chemical source term is larger than others in magnitude, indicating its role in the flame propagation. The time evolution of flame edge position and the final lift-off height are compared with measurements and generally good agreement is observed. The conditional quantities at the stabilization point reflect a balance between chemistry and micro-mixing. This investigation, which focused on model validation for various stages of spark ignition of a turbulent lifted jet flame through comparison with measurements, demonstrates that turbulent edge flame propagation in non-premixed systems can be reasonably well captured by LES/CMC.     

Theoretical analysis on the forced ignition of a quiescent mixture by repetitive heating pulse

Recently, forced ignition by nanosecond-repetitive-pulsed-discharge (NRPD) has received great attention since it can greatly promote ignition. However, there is no theoretical analysis on ignition induced by multiple heating pulses such as NRPD. Therefore, this work attempts to provide a theoretical interpretation on the ignition of a quiescent, flammable mixture by multiple discharges and to assess the effects of repetitive pulse on ignition characteristics. Based on fully transient formulation, analytical expressions describing ignition kernel propagation induced by multiple pulses are derived. The key parameters of multiple pulse heating, such as energy distribution among individual pulse, intermittent duration between neighboring pulse and total pulse numbers are appropriately incorporated, and their effects on ignition characteristics are assessed. It is found that because of memory effect, the flame kernel continues to propagate after switching off the external heating source. Sequentially introducing identical heating pulses at appropriate intermittent duration repetitively exploits the memory effect during ignition kernel development and thus extends propagating distance of the flame front. The repetitive pulse heating can promote ignition capability due to the flame revitalization effect, i.e., ignition reinforcement to the flame front thanks to additional thermal energy supplied by subsequent pulse. This is consistent with the synergistic effect of NRPD observed in previous experimental studies. In particular, as the pulse number increases, the minimum ignition energy decreases and approaches to an asymptotic value. In the limit of large pulse number, it is found that the integration of the implicit expressions describing the ignition kernel evolution does not depend on the pulse number explicitly. This substantiates the invariance of minimum ignition energy with heating pulse number. The present theory explains the effects of multiple discharges on ignition and provides insights on ignition enhancement.     

LES-CMC of high-altitude relight in an RQL aeronautical combustor

Large-Eddy Simulations with the Conditional Moment Closure sub-grid combustion model and detailed chemistry for kerosene were performed for the ignition process in an Rich-Quench-Lean aviation gas turbine combustor at high-altitude conditions. The simulations used realistic boundary conditions for the flow inlet and spray droplet size distributions and velocity. Due to the large droplets, the Central Recirculation Zone (CRZ) is filled with fuel, mostly in liquid form. The first phase of the ignition process is critical and the results show that the spark kernel must provide enough energy to evaporate the spray and pyrolyse the fuel for the flame to grow and establish in the corner of the combustor. The second phase is characterised by the flame burning the mixture in the scorner and propagating around the Inner Shear Layer. This phase is also critical, as the flame needs the prevaporised fuel and smaller droplets in the corner to sufficiently increase the temperature and be able to propagate inside the CRZ, filled with liquid fuel and cold air. If this propagation inside the CRZ is achieved, phase three is accomplished and the burner is fully ignited. The simulations demonstrate the particular importance of detailed chemistry and proper boundary conditions for flame ignition simulations in high-altitude relight conditions.     

Diesel flame lift-off stabilization in the presence of laser-ignition: a numerical study

Diesel flame lift-off and stabilization in the presence of laser-ignition were numerically investigated with the method of Eulerian stochastic fields. The aim was to scrutinise the interaction between the lifted diesel flame and an ignition kernel upstream of the lifted flame. The numerical simulation was carried out in a constant-volume combustion vessel with n-heptane as fuel. The process was studied previously in an experiment employing Diesel #2 as the fuel in the same combustion vessel. In the experiment a lifted flame was first established at a position downstream of the nozzle. An ignition kernel was then initiated using a high-energy pulse laser at a position upstream of the natural lift-off position of the diesel flame. The laser-ignition kernel was modelled using a high-temperature (～2000 K) hot spot. In both experiment and simulations the upstream front of the ignition kernel was shown to remain around the initial laser ignition site for a substantially long period of time, while the downstream front of the ignition kernel propagates rapidly towards the natural lift-off position downstream of the laser ignition site. The lift-off position oscillated before the final stabilization at the natural lift-off position. The structures and the propagation speed of the reaction fronts in the laser-ignition kernel and the main flame were analysed. Two different stabilization mechanisms, the auto-ignition mechanism and the flame propagation mechanism, were identified for the naturally lifted flame and the laser-induced reaction front, respectively. A mechanism was proposed to explain the oscillation of the lift-off position.     

Structure of propagating edge diffusion flames in hydrocarbon fuel jets

Direct numerical simulations with a C₃-chemistry model have been performed to investigate the transient behavior and internal structure of flames propagating in an axisymmetric fuel jet of methane, ethane, ethylene, acetylene, or propane in normal earth gravity (1g) and zero gravity (0g). The fuel issued from a 3-mm-i.d. tube into quasi-quiescent air for a fixed mixing time of 0.3 s before it was ignited along the centerline where the fuel–air mixture was at stoichiometry. The edge of the flame formed a vigorously burning peak reactivity spot, i.e., reaction kernel, and propagated through a flammable mixture layer, leaving behind a trailing diffusion flame. The reaction kernel broadened laterally across the flammable mixture layer and possessed characteristics of premixed flames in the direction of propagation and unique flame structure in the transverse direction. The reaction kernel grew wings on both fuel and air sides to form a triple-flame-like structure, particularly for ethylene and acetylene, whereas for alkanes, the fuel-rich wing tended to merge with the main diffusion flame zone, particularly methane. The topology of edge diffusion flames depend on the properties of fuels, particularly the rich flammability limit, and the mechanistic oxidation pathways. The transit velocity of edge diffusion flames, determined from a time series of calculated temperature field, equaled to the measured laminar flame speed of the stoichiometric fuel–air mixtures, available in the literature, independent of the gravity level.     

Laminar burning velocity and ignition delay time for premixed isooctane–air flames with syngas addition

The effects of blending syngas in different proportions to isooctane on the laminar burning velocity and ignition delay time of the fuel–air mixture have been studied in SI engine relevant conditions. The syngas is assumed to be composed of 50% H₂ and 50% CO. Simulations have been carried out using a skeletal mechanism containing 143 species and 643 reaction steps. It has been found that the blending of syngas augments the laminar burning velocity of isooctane due to increase of the thermal diffusivity of the reactant mixture and alteration in the chemistry of the flame reactions. For the mixture of 30% isooctane/70% syngas, the laminar burning velocity and the ignition delay time values are very close to those corresponding to pure isooctane. Additionally, the effects of exhaust gas recirculation have been explored for the 30% isooctane/70% syngas–air flame. It is seen that the reduction in laminar burning velocity due to the dilution by the recirculated exhaust gas can be compensated by an increase in the unburnt gas temperature. The effect of the exhaust gas dilution on the ignition delay time of 30% isooctane/70% syngas–air mixture has been found to be negligible.     

Minimum ignition energy and propagation dynamics of laminar premixed cool flames

Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.     

Parameters influencing the burning rate of laminar flames propagating into a reacting mixture

The laminar flame speed is an important property of a reacting mixture and it is used extensively for the characterization of the combustion process in practical devices. However, under engine-relevant conditions, considerable reactivity may be present in the unburned mixture, introducing thus challenges due to couplings of auto-ignition and flame propagation phenomena. In this study, the propagation of transient, one-dimensional laminar flames into a reacting unburned mixture was investigated numerically in order to identify the parameters influencing the flame burning rate in the conduction-reaction controlled regime at constant pressure. It was found that the fuel chemical classification significantly influences the burning rate. More specifically, for hydrogen flames, the “evolution” of the burning rate does not depend on the initial unburned mixture temperature. On the other hand, for n-heptane flames that exhibit low temperature chemistry, the burning rate depends on the instantaneous temperature and composition of the unburned mixture in a coupled way. A new approach was developed allowing for the decoupling the flame chemistry from the ignition dynamics as well as for the decoupling of parameters influencing the burning rate, so that meaningful sensitivity analysis could be performed. It was determined that the burning rate is not directly affected by fuel specific reactions even in the presence of low temperature chemistry whose effect is indirect through the modification of the reactants composition entering the flame. The controlling parameters include but not limited to mixture conductivity, enthalpy, and the species composition evolution in the unburned mixture.     

Dynamics of triple-flames in ignition of turbulent dual fuel mixture: A direct numerical simulation study

Pilot-ignited dual fuel combustion involves a complex transition between the pilot fuel autoignition and the premixed-like phase of combustion, which is challenging for experimental measurement and numerical modelling, and not sufficiently explored. To further understand the fundamentals of the dual fuel ignition processes, the transient ignition and subsequent flame development in a turbulent dimethyl ether (DME)/methane-air mixing layer under diesel engine-relevant conditions are studied by direct numerical simulations (DNS). Results indicate that combustion is initiated by a two-stage autoignition that involves both low-temperature and high-temperature chemistry. The first stage autoignition is initiated at the stoichiometric mixture, and then the ignition front propagates against the mixture fraction gradient into rich mixtures and eventually forms a diffusively-supported cool flame. The second stage ignition kernels are spatially distributed around the most reactive mixture fraction with a low scalar dissipation rate. Multiple triple flames are established and propagate along the stoichiometric mixture, which is proven to play an essential role in the flame developing process. The edge flames gradually get close to each other with their branches eventually connected. It is the leading lean premixed branch that initiates the steady propagating methane-air flame. The time required for the initiation of steady flame is substantially shorter than the autoignition delay time of the methane-air mixture under the same thermochemical condition. Temporal evolution of the displacement speed at the flame front is also investigated to clarify the propagation characteristics of the combustion waves. Cool flame and propagation of triple flames are also identified in this study, which are novel features of the pilot-ignited dual fuel combustion.     

An assessment of large eddy simulations of premixed flames propagating past repeated obstacles

This paper presents an assessment of Large Eddy Simulations (LES) in calculating the structure of turbulent premixed flames propagating past solid obstacles. One objective of the present study is to evaluate the LES simulations and identify the drawbacks in accounting the chemical reaction rate. Another objective is to analyse the flame structure and to calculate flame speed, generated overpressure at different time intervals following ignition of a stoichiometric propane/air mixture. The combustion chamber has built-in repeated solid obstructions to enhance the turbulence level and hence increase the flame propagating speed. Various numerical tests have also been carried out to determine the regimes of combustion at different stages of the flame propagation. These have been identified from the calculated results for the flow and flame characteristic parameters. It is found that the flame lies within the ‘thin reaction zone’ regime which supports the use of the laminar flamelet approach for modelling turbulent premixed flames. A submodel to calculate the model coefficient in the algebraic flame surface density model is implemented and examined. It is found that the LES predictions are slightly improved owing to the calculation of model coefficient by using submodel. Results are presented and discussed in this paper are for the flame structure, position, speed, generated pressure and the regimes of combustion during all stages of flame propagation from ignition to venting. The calculated results are validated against available experimental data.     

Turbulent flame propagation limits of ammonia/methane/air premixed mixture in a constant volume vessel

Ammonia is one of promising energy carriers that can be directly used as carbon-neutral fuel for combustion applications. However, because of the low-burning velocity of ammonia, it is challenging to introduce ammonia to practical combustors those are designed for general hydrocarbon fuels. One of ways to enhance the combustibility of ammonia is by mixing it with other hydrocarbon fuels, such as methane, with a burning velocity is much higher than the burning velocity of ammonia. In this study, we conducted flame propagation experiments of ammonia/methane/air using a fan-stirred constant volume vessel to clarify the effect of methane addition to ammonia on the turbulent flame propagation limit. From experimental results, we constructed the flame propagation maps and clarified the flame propagation limits. The results show that the flame propagation limits were extended with an increase in mixing a fraction of methane to ammonia. Additionally, ammonia/methane/air mixtures with the equivalence ration of 0.9 can propagate at the highest turbulent intensity, even though the peak of the laminar burning velocity is the fuel-rich side because of the diffusional-thermal instability of the flame surface. Furthermore, the Markstein number of the mixture obtained in this research successfully expressed the strength of the diffusional-thermal instability effect on the flame propagation capability. The turbulence Karlovitz number at the flame propagation limit monotonically increases with the decreasing Markstein number.     

DNS and LES of spark ignition with an automotive coil

The cycle to cycle combustion variability which is observed in spark-ignition engines is often caused by fluctuations of the early flame development. LES can be exploited for a better understanding and mastering of their origins. For that purpose appropriate models taking into account energy deposition, mixture ignition and transition to propagation are necessary requirements. This paper presents first DNS and LES of spark ignition with a real automotive coil and simplified pin-pin electrodes. The electrical circuit characteristics are provided by ISSIM while the energy deposition is modelled by Lagrangian particles. The ignition model is first evaluated in terms of initial spark radius on a pin-pin ignition experiment in pure air performed at CORIA and EM2C laboratories, showing that it pilots the radius of the torus formed by the initial shock wave. DNS of a quiescent lean propane/air mixture are then performed with this ignition system and a two-step mechanism. The impact of the modelled transferred energy during glow phase as well as the initial arc radius on the minimum ignition energy (MIE) are examined and compared to experimental values. Replacing the two-step chemistry by an analytically reduced mechanism leads to similar MIE but shows a different ignition kernel shape. Finally, LES of turbulent ignition using a Lagrangian arc model show a realistic prediction of the arc shape and its important role on the energy transfer location and thus on the flame kernel shape.     

Flame propagation in gaseous nitrous oxide

The characteristics of the propagation of a nitrous oxide decomposition flame in a tube with an internal diameter of 70 mm were measured. It was demonstrated that the pattern of flame propagation and the extent of burnout are determined by the convective motion of the flame kernel because of a very slow burning of nitrous oxide. The laminar flame speed estimated from pressure oscillograms and calculated using thermal theory of flame propagation was found to be ～1 cm/s. The critical diameter of flame quenching in channels were measured to decrease from 10 to 4 mm as the pressure was increased from 15 to 20 atm. Because of the possibility of reignition of the fresh mixture behind the flame arrester by the outflowing combustion products, the channel should be significantly longer than 200 mm.     

Measurement of unstable burning velocities of iso-octane–air mixtures at high pressure and the derivation of laminar burning velocities

A new technique is reported for measuring burning velocities at high pressures in the final stages of two inwardly propagating flame kernels in an explosion bomb. The flames were initiated at diametrically opposite spark electrodes, close to the wall, in quiescent mixtures. Measurements of pressure and flame kernel propagation speeds by high-speed photography showed the burning velocities to be elevated above the corresponding laminar burning velocities as a result of the developing flame instabilities. The enhancement increased with increase in pressure and decreased with increase in Markstein number. When the Markstein number was negative, instabilities could be appreciable, as could the enhancement. For the iso-octane–air mixtures investigated, where the mixtures had well-characterised Markstein numbers or critical Peclet numbers at the relevant pressures and temperatures, it was possible to explain the enhancement quantitatively by the spherical explosion flame instability theory of Bechtold and Matalon, provided the critical Peclet number was that observed experimentally, and allowance was made for the changing pressure. With this theoretical procedure, it was possible to derive values of laminar burning velocity from the measured values of burning velocity over a wide range of equivalence ratios, pressures, and temperatures. The values became less reliable at the higher temperatures and pressures as the data on Markstein and critical Peclet numbers became less certain. It was found that with iso-octane as the fuel the laminar burning velocity decreased during isentropic compression.     

Effects of blending 2,5-dimethylfuran on the laminar burning velocity and ignition delay time of isooctane/air mixture

The prospects of 2,5-dimethylfuran (DMF) as a bio-derived fuel that can be blended with gasoline are believed to be impressive. However, the effects of blending DMF on the key combustion parameters like the laminar burning velocity and ignition delay time of gasoline/air mixture need to be studied extensively for the successful implementation of the fuel mixture in spark ignition engines. Therefore, a skeletal chemical kinetic mechanism, comprising of 999 reactions among 218 species, has been developed in the present work for this purpose. The proposed chemical kinetic model has been validated against a wide range of experimental data for the laminar burning velocity and ignition delay time of isooctane (representing gasoline), DMF and their blends. It has been found from the present study that the thermal diffusivity of the unburnt gas mixture changes by a very small amount from the corresponding value for the pure isooctane/air mixture when DMF is added. Unlike isooctane, the DMF molecule does not consume H radicals during its primary breakup. Therefore, the maximum laminar burning velocity increases marginally when 50% DMF is blended with isooctane due to the increased presence of H radicals in the flame. The negative temperature coefficient behaviour in the ignition delay time of the isooctane fuel vanishes when 30% DMF (v/v) is blended to it.     

Spatio-temporal evolution of cavity ignition in supersonic flow

Successful ignition in the recirculating flow of a scramjet flame holder can be highly dependent upon the location of energy deposition because of the spatial variation of fuel concentration and flow properties. The current work experimentally investigated ignition processes when energy was deposited (～100 mJ) via a spark discharge at four locations in the base of a cavity or by laser-induced breakdown in a Mach 2 flow with a stagnation temperature and pressure of 590 K and 483 kPa, respectively. The cavity was directly fueled with ethylene injection. The time dependent heat release was imaged at 40,000 frames per second and fuel concentration and distribution measurements were taken in the cavity prior to ignition. The average fuel concentration at the lean and rich ignition limits near the energy deposition locations measured 4.4–9.3% (Φ= 0.75 to 1.47). Energy deposition near the cavity step resulted in near immediate ignition kernel development and rapid achievement of self-sustained flame propagation in the front of the cavity, often faster than the bulk recirculation time of the cavity, leading to a spike in heat release. Energy deposition away from the cavity step region led to competition between local flow velocity, fuel concentration, and flame propagation rates. Ignition kernels formed along the floor of the cavity towards the closeout ramp and were rapidly advected towards the cavity step region before flame propagation could ensue. The fastest and most robust ignition events for all fueling cases showed rapid spanwise flame propagation near the cavity step.     

2D computations of FREI with cool flames for n-heptane/air mixture

Two-dimensional axisymmetric numerical simulation reproduced flames with repetitive extinction and ignition (FREI) in a micro flow reactor with a controlled temperature profile with a stoichiometric n-heptane/air mixture, which have been observed in the experiment. The ignition of hot flame occurred from consumption reactions of CO that was remained in the previous cycle of FREI. Between extinction and ignition locations of hot flames, several other heat release rate peaks related to cool and blue flames were observed for the first time. After the extinction of the hot flame, cool flame by the low-temperature oxidation of n-heptane appeared first and was stabilized in a low wall temperature region. In the downstream of the stable cool flame, a blue flame by the consumption reactions of cool flame products of CH₂O and H₂O₂ appeared. After that, the hot flame ignition occurred from the remaining CO in the downstream of the blue flame. Then after the next hot flame ignition, the blue flame was swept away by the propagating hot flame. Soon before the hot flame merged with the stable cool flame, the hot flame propagation was intensified by the cool flame. After the hot flame merged with the stable cool flame, the hot flame reacted with the incoming fresh mixture of n-C₇H₁₆ and O₂.