Studies of the dynamics of autoignition assisted outwardly propagating spherical cool and double flames under shock-tube conditions

The initiation, propagation, and transition of the autoignition assisted spherical cool flame and double flame are studied numerically and experimentally using n-heptane/air/He mixtures under shock-tube experimental conditions over a wide range of temperatures. The primary goal of the current study is to understand the effects of the ignition Damkohler number, ignition energy, flame curvature, and autoignition-induced flow compression on the propagation of spherical flames to ensure the proper interpretation of shock-tube flame speed measurements at engine-relevant conditions. The results show that at high ignition Damkohler number, there are three different flame regimes, cool flame, double flame, and hot flame. The cool flame speed accelerates dramatically with the increase of ignition Damkohler number. In addition, it is found that the change of flame regime, low-temperature autoignition, flame stretch, and autoignition-induced flow compression result in a complicated non-linear dependence of flame speed on stretch. The results also reveal that the spherical cool flame has much lower Markstein length compared to the hot flame at T > 600 K. Moreover, it is found that both the autoignition assisted cool flame and the trailing hot flame front in the double flame can propagate much faster that the hot flame alone at the same mixture conditions, leading to a nonlinear dependence of flame speed on the mixture initial temperature. The simulated flame trajectories and the flame speed dependence on temperature agree qualitatively well with the shock-tube experiments. A quantitative criterion to ensure the accurate speed measurement of the cool and hot flame is proposed. The present study provides important physical insight and guidance for the flame speed measurement using a shock-tube at engine relevant conditions.     

Understanding cool flames and warm flames

Low-temperature flames such as cool flames, warm flames, double flames, and auto-ignition assisted flames play a critical role in the performance of advanced engines and fuel design. In this paper, an overview of the recent progresses in understanding low-temperature flames and dynamics as well as their impacts on combustion, advanced engines, and fuel development will be presented. Specifically, at first, a brief review of the history of cool flames is made. Then, the recent experimental studies and computational modeling of the flame structures, dynamics, and burning limits of non-premixed and premixed cool flames, warm flames, and double flames are presented. The flammability limit diagram and the temperature-dependent chain-branching reaction pathways, respectively, for hot, warm, and cool flames at elevated temperature and pressure will be discussed and analyzed. After that, the effect of low temperature auto-ignition of auto-igniting mixtures at high ignition Damköhler numbers at engine conditions on the propagation of cool flames, warm flames, and double flames as well as turbulent flames will be discussed. Finally, a new platform using low temperature flames for the development and validation of chemical kinetic models of alternative fuels will be presented. Discussions of future research of the dynamics and control of low temperature flames under engine conditions will be made.     

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.     

Effects of pressure rise rate on laminar flame speed under normal and engine-relevant conditions

Laminar flame speed (LFS) is one of the most important physicochemical properties of a combustible mixture. At normal and elevated temperatures and pressures, LFS can be measured using propagating spherical flames in a closed chamber. LFS is also used in certain turbulent premixed flame modelling for combustion in spark ignition engines. Inside the closed chamber or engine, transient pressure rise occurs during the premixed flame propagation. The effects of pressure rise rate (PRR) on LFS are examined numerically in this study. One-dimensional simulations are conducted for spherical flame propagation in a closed chamber. Detailed chemistry and transport are considered. Different values of PRR at the same temperature and pressure are achieved through changing the spherical chamber size. It is found that the effect of PRR on LFS is negligible under the normal and engine-relevant conditions considered in this study. This observation is then explained through the comparison between the unsteady and convection terms in the energy equation for a premixed flame.     

Numerical simulation of a mixed-mode reaction front in a PPC engine

The ignition process, mode of combustion and reaction front propagation in a partially premixed combustion (PPC) engine running with a primary reference fuel (87% iso-octane, 13% n-heptane by volume) is studied numerically in a large eddy simulation. Different combustion modes, ignition front propagation, premixed flame and non-premixed flame, are observed simultaneously. Displacement speed of CO iso-surface propagation describes the transition of premixed auto-ignition to non-premixed flame. High temporal resolution optical data of CH₂O and chemiluminescence are compared with simulated results. A high speed ignition front is seen to expand through fuel-rich mixture and stabilize around stoichiometry in a non-premixed flame while lean premixed combustion occurs in the spray wake at a much slower pace. A good qualitative agreement of the distribution of chemiluminescence and CH₂O formation and destruction shows that the simulation approach sufficiently captures the driving physics of mixed-mode combustion in PPC engines. The study shows that the transition from auto-ignition to flame occurs over a period of several crank angles and the reaction front propagation can be captured using the described model.     

Dynamics and burning limits of near-limit hot,warm, and cool diffusion flames of dimethyl ether at elevated pressures

The near-limit diffusion flame regimes and extinction limits of dimethyl ether at elevated pressures and temperatures are examined numerically in the counterflow geometry with and without radiation at different oxygen concentrations. It is found that there are three different flame regimes—hot flame, warm flame, and cool flame—which exist, respectively, at high, intermediate, and low temperatures. Furthermore, they are governed by three distinct chain-branching reaction pathways. The results demonstrate that the warm flame has a double reaction zone structure and plays a critical role in the transition between cool and hot flames. It is also shown that the cool flame can be formed in several different ways: by either radiative extinction or stretch extinction of a hot flame or by stretch extinction of a warm flame. A warm flame can also be formed by radiative extinction of a hot flame or ignition of a cool flame. A general €-shaped flammability diagram showing the burning limits of all three flame regimes at different oxygen mole fractions is obtained. The results show that thermal radiation, reactant concentration, temperature, and pressure all have significant impacts on the flammable regions of the three flame regimes. Increases in oxidizer temperature, oxygen concentration, and pressure shift the cool flame regime to higher stretch rates and cause the warm flame to have two extinction limits. At elevated temperatures, it is found that there is a direct transition between the hot flame and warm flame at low stretch rates. The results also show that, unlike the hot flame, the cool flame structure cannot be scaled by using pressure-weighted stretch rates due to the its significant reactant leakage and strong dependence of reactivity on pressure. The present results advance the understanding of near-limit flame dynamics and provide guidance for experimental observation of different flame regimes.     

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₂.     

Ignition of dimethyl ether/air mixtures by hot particles: Impact of low temperature chemical reactions

Understanding and characterizing ignition of flammable mixtures by hot particles is important for assessing and reducing the risk of accidental ignition and explosion in industry and aviation. Recently, many studies have been conducted for ignition of gaseous mixtures by hot particles. However, the effects of low-temperature chemistry (LTC) on ignition by hot particles received little attention. LTC plays an important role in the ignition of most hydrocarbon fuels and may induce cool flames. The present study aims to numerically assess the effects of LTC on ignition by the hot particles. We consider the transient ignition processes induced by a hot spherical particle in quiescent and flowing stoichiometric dimethyl ether/air mixtures. 1D and 2D simulations, respectively, are conducted for the ignition process by hot-particles in quiescent and flowing mixtures. A detailed kinetic model including both LTC and high-temperature chemistry (HTC) is used in simulations. The results exhibit a premixed cool flame to be first initiated by the hot particle. Then a double-flame structure with both premixed cool and hot flames is observed at certain conditions. At zero or low inlet flow velocities, the hot flame catches up and merges with the leading cool flame. At high inlet flow velocities, the hot flame cannot be initiated due to the short residence time and large convective loss of heat and radicals. Comparing the results with and without considering LTC confirms that LTC accelerates substantially ignition via HTC in a certain range of hot particle temperatures. The mechanism of ignition promotion by LTC is interpreted by analyzing the radical pool produced by the LTC and HTC surrounding the hot particle. Moreover, the influence of inlet flow velocity on ignition by hot particles is assessed. Non-monotonic change of ignition delay time with flow velocity is observed and discussed.     

超声速预混可燃气流的点火与燃烧

在激波风洞一激波管组合设备上开展了碳氢燃料超声速预混可燃气流的点火与燃烧实验研究。实验结果表明:利用激波对燃料进行预热,并以高温燃气作为引导火焰,可以有效缩短汽油空气超声速可燃混气的点火延迟时间,使之缩短到 0.2 ms以下。利用纹影照片对超声速燃烧流场结构作出了分析;研究了超声速预混可燃气流的温度以及当量比对超声速燃烧流场结构、点火与火焰传播特性的影响。     

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.     

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.     

Numerical simulation and validation of SI-CAI hybrid combustion in a CAI/HCCI gasoline engine

SI-CAI hybrid combustion, also known as spark-assisted compression ignition (SACI), is a promising concept to extend the operating range of CAI (Controlled Auto-Ignition) and achieve the smooth transition between spark ignition (SI) and CAI in the gasoline engine. In this study, a SI-CAI hybrid combustion model (HCM) has been constructed on the basis of the 3-Zones Extended Coherent Flame Model (ECFM3Z). An ignition model is included to initiate the ECFM3Z calculation and induce the flame propagation. In order to precisely depict the subsequent auto-ignition process of the unburned fuel and air mixture independently after the initiation of flame propagation, the tabulated chemistry concept is adopted to describe the auto-ignition chemistry. The methodology for extracting tabulated parameters from the chemical kinetics calculations is developed so that both cool flame reactions and main auto-ignition combustion can be well captured under a wider range of thermodynamic conditions. The SI-CAI hybrid combustion model (HCM) is then applied in the three-dimensional computational fluid dynamics (3-D CFD) engine simulation. The simulation results are compared with the experimental data obtained from a single cylinder VVA engine. The detailed analysis of the simulations demonstrates that the SI-CAI hybrid combustion process is characterised with the early flame propagation and subsequent multi-site auto-ignition around the main flame front, which is consistent with the optical results reported by other researchers. Besides, the systematic study of the in-cylinder condition reveals the influence mechanism of the early flame propagation on the subsequent auto-ignition.     

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.     

Measurements of the critical initiation radius and unsteady propagation of n-decane/air premixed flames

Unsteady flame propagation, the critical radius for flame initiation, and multiple flame regimes of n-decane/air mixtures are studied experimentally and computationally using outwardly propagating spherical flames at various equivalence ratios and pressures. The transient flame speeds, trajectories, and critical radius are measured. The experimental results are compared with direct numerical simulations using detailed high temperature kinetic models. Both experimental and numerical results show that there exist multiple flame regimes in the unsteady spherical flame initiation process. The transition between the flame regimes depends strongly on the mixture equivalence ratio (or Lewis number). It is found that there is a critical flame radius and that it increases dramatically as the mixture equivalence ratio and pressure decrease. The large increase of critical flame radius leads to a dramatic increase of the minimum ignition energy. Furthermore, the flame thickness and the radical pool concentration change significantly during the transition from the ignition flame regime to the self-sustained propagating flame regime. For the same steady state flame speeds, the predicted unsteady flame speeds and the critical flame radius differ significantly from the experimental results. Moreover, different chemical kinetic mechanisms predict different unsteady flame speeds. The existence of multiple flame regimes and the large critical radius of lean liquid fuel mixtures make the ignition of lean mixtures at low pressure and the development of a validated kinetic model more challenging. The unsteady flame regimes, speeds, and critical flame radius should be included as targets of future kinetic model development for turbulent combustion modeling.     

Transient interactions between a premixed double flame and a vortex

The interaction between a laminar flame and a vortex is an important study for understanding the fundamentals of turbulent combustion. In the past, however, flame-vortex interactions have been investigated only for high-temperature flames. In this study, the impact of a vortex on a premixed double flame, which consists of a coupled cool flame and a hot flame, is examined experimentally and computationally using dimethyl ether/oxygen/ozone mixtures. The double flame is first shown to occur near the extinction limit of the hot flame. The differences between steady-state cool flames, double flames, and hot flames are explored in a one-dimensional counterflow configuration. The transient interactions between double flames and impinging vortices are then investigated experimentally using a micro-jet and numerically in two-dimensional transient modeling. It is seen that the vortex can extinguish the near-limit hot flame locally, resulting in a lone cool flame. At higher vortex intensities, the cool flame may also be extinguished after the extinction of the hot flame. It is found that there can be three different transient flame structures coexisting at the same time: an extinguished flame hole, a cool flame, and a double flame. Moreover, flame curvature is shown to play an important role in determining whether the vortex weakens or strengthens the cool flame and double flame.     

Premixed flames for arbitrary combinations of strain and curvature

Many modeling strategies for combustion rely on laminar flamelet concepts to determine structure and properties of multi-dimensional and turbulent flames. Using flamelet tabulation strategies, the user anticipates certain aspects of the combustion process prior to the simulation and selects a flamelet model which mimics local flame conditions in the more complex configuration. Flame stretch, which can be decomposed into contributions from strain and curvature, is one of the conditions influencing a flame’s properties, structure, and stability. The objective of this work is to study premixed flame structures in the strain-curvature space using a recently published composition space model (CSM) and three physical space models for canonical flame configurations (stagnation flame, spherical expanding flame and inwardly propagating flame). Flames with effective Lewis numbers both smaller and larger than unity are considered. For canonical laminar flames, the stretch components are inherently determined through boundary conditions and their specific flame configuration. Therefore, canonical flames can only represent a certain sub-set of stretch effects experienced by multi-dimensional and turbulent flames. On the contrary, the CSM allows arbitrary combinations of strain and curvature to be prescribed for premixed flames exceeding the conditions attainable with the canonical flame setups. Thereby, also influences of negative strain effects and large curvatures can be studied. A parameter variation with the CSM shows that flame structures still significantly change outside the region of the canonical flame configurations. Furthermore, limits in the strain-curvature space are discussed. The present paper highlights advantages of composition space modeling which is achieved by detaching the representation of the flame structure from a specific canonical flame configuration in physical space.     

Kinetic effects of NO addition on n-dodecane cool and warm diffusion flames

The kinetic effects of NO addition on the flame dynamics and burning limits of n-dodecane cool and warm diffusion flames are investigated experimentally and computationally using a counterflow system. The results show that NO plays different roles in cool and warm flames due to their different reaction pathway sensitivities to the flame temperature and interactions with NO. We observe that NO addition decreases the cool flame extinction limit, delays the extinction transition from warm flame to cool flame, and promotes the ignition transition from warm flame to hot flame. In addition, jet-stirred reactor (JSR) experiments of n-dodecane oxidation with and without NO addition are also performed to develop and validate a n-dodecane/NO_x kinetic model. Reaction pathway and sensitivity analyses reveal that, for cool flames, NO addition inhibits the low-temperature oxidation of n-dodecane and reduces the flame temperature due to the consumption of RO₂ via NO+RO₂?NO₂+RO, which competes with the isomerization reaction that continues the peroxy radical branching sequence. The model prediction captures well the experimental trend of the inhibiting effect of NO on the cool flame extinction limit. For warm flames, two different kinds of warm flame transitions, the warm flame extinction transition to cool flame and the warm flame reignition transition to hot flame, were observed. The results suggest that warm extinction transition to cool flame is suppressed by NO addition while the warm flame reignition transition to hot flame is promoted. The kinetic model developed captures well the experimentally observed warm flame transitions to cool flame but fails to predict the warm flame reignition to hot flame at similar experimental conditions.     

Studies of autoignition-assisted nonpremixed cool flames

Autoignition-assisted nonpremixed cool flames of diethyl ether (DEE) are investigated in both laminar counterflow and turbulent jet flame configurations. First, the ignition and extinction limits of laminar nonpremixed cool flames of diluted DEE are measured and simulated using detailed kinetic models. The laminar flame measurements are used to validate the kinetic models and guide the turbulent flame measurements. The results show that, below a critical mixture condition, for elevated temperature and dilute mixtures, the cool flame extinction limit and the low-temperature ignition limit merge, leading to autoignition-assisted cool flame stabilization without hysteresis. Based on the findings from the laminar flame experiments, autoignition-assisted turbulent lifted cool flames are established using a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner. The lift-off heights of the turbulent cool flames are quantified using formaldehyde planar laser-induced fluorescence. Based on an analogy with autoignition-assisted lifted hot flames, a correlation is proposed such that the autoignition-assisted cool flame lift-off height scales with the product of the flow velocity and the square of the first-stage ignition delay time. Using this scaling, we demonstrate that the kinetic mechanism that most accurately predicts the laminar flame ignition and extinction limits also best predicts the turbulent cool flame lift-off height.