Simulations and experiments on the ignition probability in turbulent premixed bluff-body flames

The ignition characteristics of a premixed bluff-body burner under lean conditions were investigated experimentally and numerically with a physical model focusing on ignition probability. Visualisation of the flame with a 5 kHz OH* chemiluminescence camera confirmed that successful ignitions were those associated with the movement of the kernel upstream, consistent with previous work on non-premixed systems. Performing many separate ignition trials at the same spark position and flow conditions resulted in a quantification of the ignition probability P_ign, which was found to decrease with increasing distance downstream of the bluff body and a decrease in equivalence ratio. Flows corresponding to flames close to the blow-off limit could not be ignited, although such flames were stable if reached from a richer already ignited condition. A detailed comparison with the local Karlovitz number and the mean velocity showed that regions of high P_ign are associated with low Ka and negative bulk velocity (i.e. towards the bluff body), although a direct correlation was not possible. A modelling effort that takes convection and localised flame quenching into account by tracking stochastic virtual flame particles, previously validated for non-premixed and spray ignition, was used to estimate the ignition probability. The applicability of this approach to premixed flows was first evaluated by investigating the model's flame propagation mechanism in a uniform turbulence field, which showed that the model reproduces the bending behaviour of the S_T-versus-u′ curve. Then ignition simulations of the bluff-body burner were carried out. The ignition probability map was computed and it was found that the model reproduces all main trends found in the experimental study.     

Numerical study of the ignition behavior of a post-discharge kernel in a turbulent stratified crossflow

Ensuring robust ignition is critical for the operability of aeronautical gas-turbine combustors. For ignition to be successful, an important aspect is the ability of the hot gas generated by the spark discharge to initiate combustion reactions, leading to the formation of a self-sustained ignition kernel. This study focuses on this phenomena by performing simulations of kernel ignition in a crossflow configuration that was characterized experimentally. First, inert simulations are performed to identify numerical parameters correctly reproducing the kernel ejection from the ignition cavity, which is here modeled as a pulsed jet. In particular, the kernel diameter and the transit time of the kernel to the reacting mixture are matched with measurements. Considering stochastic perturbations of the ejection velocity of the ignition kernel, the variability of the kernel transit time is also reproduced by the simulations. Subsequently, simulations of a series of ignition sequences are performed with varying equivalence ratio of the fuel-air mixture in the crossflow. The numerical results are shown to reproduce the ignition failure that occurs for the leanest equivalence ratio ($?=0.6$). For higher equivalence ratios, the simulations are shown to capture the sensitivity of the ignition to the equivalence ratio, and the kernel successfully transitions into a propagating flame. Significant stochastic dispersion of the ignition strength is observed, which relates to the variability of the transit time of the kernel to the reactive mixture. An analysis of the structure of the ignition kernel also highlights the transition towards a self-propagating flame for successful ignition conditions.     

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

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

Simultaneous OH-PLIF and PIV measurements in a gas turbine model combustor

In highly turbulent environments, combustion is strongly influenced by the effects of turbulence chemistry interactions. Simultaneous measurement of the flow field and flame is, therefore, obligatory for a clear understanding of the underlying mechanisms. In the current studies simultaneous PIV and OH-PLIF measurements were conducted in an enclosed gas turbine model combustor for investigating the influence of turbulence on local flame characteristics. The swirling CH₄/air flame that was investigated had a thermal power of 10.3 kW with an overall equivalence ratio of ϕ=0.75 and exhibited strong thermoacoustic oscillations at a frequency of approximately 295 Hz. The measurements reveal the formation of reaction zones at regions where hot burned gas from the recirculation zones mixes with the fresh fuel/air mixture at the nozzle exit. However, this does not seem to be a steady phenomenon as there always exist regions where the mixture has failed to ignite, possibly due to the high local strain rates present, resulting in small residence time available for a successful kinetic runaway to take place. The time averaged PIV images showed flow fields typical of enclosed swirl burners, namely a big inner recirculation zone and a small outer recirculation zone. However, the instantaneous images show the existence of small vortical structures close to the shear layers. These small vortical structures are seen playing a vital role in the formation and destruction of reaction zone structures. One does not see a smooth laminar flame front in the instantaneous OH-PLIF images, instead isolated regions of ignition and extinction highlighting the strong interplay between turbulence and chemical reactions. PACS 33.20.-t; 33.50.-j; 47.27.-i; 47.32.Ef; 47.70.Pq; 82.33.Vx; 82.40.-g     

Catalytic ignition and autothermal combustion of JP-8 and its surrogates over a Pt/γ-Al2O3 catalyst

With the aim of utilizing JP-8 fuel for small scale portable power generation systems, catalytic combustion of JP-8 is studied. The surface ignition, extinction and autothermal combustion of JP-8, of a six-component surrogate fuel mixture, and the individual components of the surrogate fuel over a Pt/γ-Al₂O₃ catalyst are experimentally investigated in a packed bed flow reactor. The surrogate mixture exhibits similar ignition–extinction behavior and autothermal temperatures compared to JP-8 suggesting the possibility of using this surrogate mixture for detailed kinetics of catalytic combustion of JP-8. It is shown that JP-8 ignites at low temperatures in the presence of catalyst. Upon ignition, catalytic combustion of JP-8 and the surrogate mixture is self-sustained and robust combustion is observed under fuel lean as well as fuel rich conditions. It is shown that the ignition temperature of the hydrocarbon fuels increases with increasing equivalence ratio. Extinction is observed under fuel lean conditions, whereas sustained combustion was also observed for fuel rich conditions. The effect of dilution in the air flow on the catalytic ignition and autothermal temperatures of the fuel mixture is also investigated by adding helium to the air stream while keeping the flow rate and the equivalence ratio constant. The autothermal temperature decreases linearly as the amount of dilution in the flow is increased, whereas the ignition temperature shows no dependence on the dilution level under the range of our conditions, showing that ignition is dependent only on the type and relative concentration of the active species.     

涡流控制阀对四气门汽油机缸内空气运动影响的研究

在一台单缸四气门汽油机上,用激光多普勒测速仪实验研究了涡流控制阀打开和关闭后,缸内平均速度、湍流强度和湍流结构的变化．结果表明：当涡流控制阀打开时,在压缩后期缸内流动主要是滚流的畸变和破碎过程,而当涡流控制阀关闭后,在压缩后期缸内流动主要是涡流运动．两种情况下的湍流强度和湍流结构也有较大差别．     

Towards the development of an efficient low-NOx ammonia combustor for a micro gas turbine

Recent studies have demonstrated stable generation of power from pure ammonia combustion in a micro gas turbine (MGT) with a high combustion efficiency, thus overcoming some of the challenges that discouraged such applications of ammonia in the past. However, achievement of low NOx emission from ammonia combustors remains an important challenge. In this study, combustion techniques and combustor design for efficient combustion and low NOx emission from an ammonia MGT swirl combustor are proposed. The effects of fuel injection angle, combustor inlet temperature, equivalence ratio, and ambient pressure on flame stabilization and emissions were investigated in a laboratory high pressure combustion chamber. An FTIR gas analyser was employed in analysing the exhaust gases. Numerical modeling using OpenFOAM was done to better understand the dependence of NO emissions on the equivalence ratio. The result show that inclined fuel injection as opposed to vertical injection along the combustor central axis resulted to improved flame stability, and lower NH₃ and NOx emissions. Numerical and experimental results showed that a control of the equivalence ratio upstream of the combustor is critical for low NOx emission in a rich-lean ammonia combustor. NO emission had a minimum value at an upstream equivalence ratio of 1.10 in the experiments. Furthermore, NO emission was found to decrease with ambient pressure, especially for premixed combustion. For the rich-lean combustion strategy employed in this study, lower NOx emission was recorded in premixed combustion than in non-premixed combustion indicating the importance of mixture uniformity for low NOx emission from ammonia combustion. A prototype liner developed to enhance the control and uniformity of the equivalence ratio upstream of the combustor further improved ammonia combustion. With the proposed liner design, NOx emission of 42?ppmv and ammonia combustion efficiency of 99.5% were achieved at 0.3?MPa for fuel input power of 31.44?kW.     

Investigation on the flameholding mechanisms in supersonic flows: backward-facing step and cavity flameholder

Abstract  

As effective devices to extend the fuel residence time in supersonic flow and prolong the duration time for hypersonic vehicles cruising in the near-space with power, the backward-facing step and the cavity are widely employed in hypersonic airbreathing propulsive systems as flameholders. The two-dimensional coupled implicit RANS equations, the standard k-ε turbulence model, and the finite-rate/eddy-dissipation reaction model have been used to generate the flow field structures in the scramjet combustors with the backward-facing step and the cavity flameholders. The flameholding mechanism in the combustor has been investigated by comparing the flow field in the corner region of the backward-facing step with that around the cavity flameholder. The obtained results show that the numerical simulation results are in good agreement with the experimental data, and the different grid scales make only a slight difference to the numerical results. The vortices formed in the corner region of the backward-facing step, in the cavity and upstream of the fuel injector make a large difference to the enhancement of the mixing between the fuel and the free airstream, and they can prolong the residence time of the mixture and improve the combustion efficiency in the supersonic flow. The size of the recirculation zone in the scramjet combustor partially depends on the distance between the injection and the leading edge of the cavity. Further, the shock waves in the scramjet combustor with the cavity flameholder are much stronger than those that occur in the scramjet combustor with the backward-facing step, and this causes a large increase in the static pressure along the walls of the combustor.     

Simulation of spray development and turbulent combustion processes in low and high speed diesel engines by the CMC-ISR model

Simulation is performed to analyse the characteristics of turbulent spray combustion in conventional low and high speed diesel engine conditions. Turbulence–chemistry interaction is resolved by the Conditional Moment Closure (CMC) model in the spatially integrated form of an Incompletely Stirred Reactor (ISR). After validation against measured pressure traces, characteristic length and time scales and dimensionless numbers are estimated at the locations of sequentially injected fuel groups. Conditional flame structures are calculated for sequentially evaporated fuel groups to consider different available periods for ignition chemistry. Injection overlaps the combustion period in the high rpm engine, while most combustion occurs after injection and evaporation are complete in the low rpm engine. Ignition occurs in rich premixture with the initial peak temperature at the equivalence ratio around 2–4 as observed in Dec [2]. It corresponds to the most reactive mixture fraction of the minimum ignition delay for the given mixture states. Combustion proceeds to lean and rich sides in the mixture fraction space as a diffusion process by turbulence. The mean scalar dissipation rates (SDRs) are lower than the extinction limit to show stability of diffusion flames throughout the combustion period.     

Analysis of first stage ignition delay times of dimethyl ether in a laminar flow reactor

The combustion chemistry of the first stage ignition and chemistry/flow interactions are studied for dimethyl ether (DME) with a mathematical analysis of two systems: a plug flow reactor study is used to reduce the reaction chemistry systematically. A skeletal reaction mechanism for the low temperature chemistry of DME until the onset of ignition is derived on the basis of the detailed DME mechanism of the Lawrence Livermore National Laboratory – see Curran, Fischer and Dryer, Int. J. Chem. Kinetics, Vol. 32 (2000). It is shown that reasonably good results for ignition delay times can be reached using a simple system of three ordinary differential equations and that the resulting analytical solution depends only on two reaction rates and the initial fuel concentration. The stepwise reduction of the system based on assumptions yields an understanding on why these reactions are so important. Furthermore, the validation of the assumptions yields insight into the influence of the fuel and the oxygen concentration on the temperature during the induction phase. To investigate the influence of chemistry/flow interactions, a 2D model with a laminar Hagen–Poiseuille flow and 2D-polynomial profiles for the radial species concentration is considered. For the 2D model, it is found that only the diffusion coefficients and the reactor radius need to be taken into consideration additionally to describe the system sufficiently. Also, the coupling of flow and chemistry is clarified in the mathematical analysis. The insight obtained from the comparison of the 2D model and the plug flow model is used to establish an average velocity for the conversion of ignition locations to ignition delay times in a laminar flow reactor. Finally, the 2D analytical solution is compared against new experimental data, obtained in such a laminar flow reactor for an undiluted DME/air mixture with an equivalence ratio of φ = 0.835 and a temperature range of 555 to 585 K at atmospheric pressure.     

燃油分级多点喷射低污染燃烧室的化学反应网络模型分析

本文采用基于详细化学反应机理的化学反应网络模型分析了航空发动机燃油径向分级多点喷射低污染燃烧室的NO_x排放特性。该分级燃烧室不同于传统燃烧室,头部由值班区和主燃区两个不同的燃烧区域,根据CFD得到的流场特性和当量比的分布特性对燃烧室进行分区构建化学反应器网络模型,研究了值班级当量比以及值班级和主燃级两级供油比例对排放的影响。同时,还分析了空气进口温度对NO_x排放的影响。得到了较为合理的变化趋势,为低污染燃烧室的初步设计提供了有益的指导。     

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.     

A parametric study of ignition dynamics at ECN Spray A thermochemical conditions using 2D DNS

The ignition process in diesel engines is highly complex and incompletely understood. In the present study, two-dimensional direct numerical simulations are performed to investigate the ignition dynamics and their sensitivity to thermochemical and mixing parameters. The thermochemical and mixing conditions are matched to the benchmark Spray A experiment from the Engine Combustion Network. The results reveal a complex ignition process with overlapping stages of: low-temperature ignition (cool flames), rich premixed ignition, and nonpremixed ignition, which are qualitatively consistent with prior experimental and numerical investigations, however, this is the first time that fully-resolved simulations have been reported at the actual Spray A thermochemical condition. Parametric variations are then performed for the Damköhler number Da, oxidiser temperature, oxygen concentration, and peak mixture fraction (a measure of premixedness), to study their effect on the ignition dynamics. It is observed that with both increasing oxidiser temperature and decreasing oxygen concentration, that the cool flame moves to richer mixtures, the overlap in the ignition stages decreases, and the (nondimensional) time taken to reach a fully burning state increases. With increasing Da, the cool-flame speed is decreased due to lower mean mixing rates, which causes a delayed onset of high-temperature ignition. With increasing peak mixture fraction, the onset of each stage of ignition is not affected, but the overall duration of the ignition increases leading to a longer burn duration. Overall, the results suggest that turbulence–chemistry interactions play a significant role in determining the timing and location in composition space of the entire ignition process.     

Ammonia–methane interaction in jet-stirred and flow reactors: An experimental and kinetic modeling study

The influence of the addition of ammonia on the oxidation of methane was investigated both experimentally and numerically. Experiments were carried out at atmospheric pressure, using a fused silica jet-stirred reactor, and a recrystallized alumina tubular reactor designed on purpose to reach temperatures as high as ～2000 K. A temperature range of 600–1200 K was investigated in the jet-stirred reactor at a residence time of 1.5 s, while experiments in the flow reactor were carried out between 1200 and 2000 K, for a fixed residence time of about 25 ms in the reactive zone. A methane/ammonia mixture, diluted in helium, was used in both reactors with equivalence ratios varied between 0.5 and 2 in the first reactor, while stoichiometric conditions were investigated in the second one. The measurements indicate that CH₄ reactivity was promoted by NH₃ addition below 1200 K, but not so much influenced above. These results were interpreted and explained using a comprehensive kinetic model, previously validated over a wider range of operating conditions. The mechanism allowed to shed light on the underlying causes of the anticipated methane reactivity at low temperature, and of the major role played by NO_x in it. This effect was shown to become less significant at higher temperatures, where the reactivity is mainly governed by H-abstractions on both fuels.     

Generalized entropy production analysis for mechanism reduction

The paper introduces a generalized formulation for the computation of the relative contribution of each elementary reaction to the total entropy production, which has been proposed as a measure of the importance of elementary reactions and used for the reduction of detailed chemical reaction mechanisms. The reduction method is extended for the cases where the principle of detailed balance does not hold or apply, namely in the case of irreversible reactions or when the reverse rate constants are not computed via the thermodynamic equilibrium constants. Using a mechanism for n-butane consisting exclusively of reversible reactions, the new formulation is compared to the original one, and then applied for the construction of a skeletal mechanism for n-dodecane starting from a detailed mechanism which includes predominantly irreversible reactions. The skeletal scheme is found to accurately capture the ignition delay times over an extended range of pressure, initial temperature and equivalence ratio, the steady-state temperature as function of the residence time in a non-isothermal adiabatic perfectly stirred reactor, and the laminar flame speed of atmospheric flames at different unburned mixture temperatures and equivalence ratios.     

Adjoint-based sensitivity and ignition threshold mapping in a turbulent mixing layer

Successful ignition in non-premixed turbulent flows remains a fundamental challenge in combustion systems. Current design strategies typically rely on iterative testing to map the spatial distribution of ignition probability. We propose to accelerate this by formulating the adjoint of the perturbed and linearised governing equations in such a way that sensitivity of an ignition indicator can be obtained with a cost comparable to the flow solution. A space–time discrete adjoint method for multi-component chemically reacting flows is developed, and the gradient formed via the corresponding adjoint solution is used to identify regions favourable to ignition in a direct numerical simulation of non-premixed turbulent free shear flow. This approach requires a specific definition of an ignition metric, although this can be problematic because ignition either succeeds or fails after some period and thus gradients for some metrics become ill-defined near the ignition threshold. To this end, a quantity of interest is designed to provide short-time sensitivity in conjunction with an indicator function over a long-time period that informs whether successful ignition occurred. The gradients are used in a line-search algorithm to map the ignition boundary under specific constraints. Finally, parametric sensitivity is evaluated at different flow realisations to analyse factors governing local sensitivity in unsteady chemically reacting flows.     

Simulation of a round supersonic combustor using wall-modeled large eddy simulation and partially-stirred reactor models

Simulations are presented for a generic, round supersonic combustor. Turbulence is modeled in the combustor using a wall-modeled large eddy simulation approach. Combustion is modeled using a small quasi-global mechanism and a more detailed skeletal mechanism. Both mechanisms are used in conjunction with two variations of the partially-stirred reactor model for sub-grid turbulence chemistry interactions. Sensitivity of the solutions to grid resolution is investigated. It is found that in order to achieve reasonable grid convergence in the mean wall pressure, the model constant that appears in the partially-stirred reactor model must be a function of both the chemistry mechanism and the grid resolution. Most of the combinations of mechanism and turbulent combustion model tested can be tuned in order to predict the location of the pre-combustion shock train and the peak mean pressure in the combustor. It is found that while the different models are able to reproduce the mean wall pressure, there are significant differences in the mean temperature and heat release rate fields. The sensitivity of the different combinations of mechanisms and partially-stirred reactor formulation is quantified and some combinations are found to be more prone to blowout. Two of the tuned models were tested across several fuel equivalence ratios with a single value of the partially-stirred reactor model constant. One model produced reasonable predictions of shock location and peak mean pressure for each equivalence ratio. The second model captured the global trends in the mean wall pressure, but was unable to quantitatively predict the shock location and peak mean pressure for all equivalence ratios tested.     

Application of an aerosol shock tube to the measurement of diesel ignition delay times

Shock tube ignition delay times were measured for DF-2 diesel/21% O₂/argon mixtures at pressures from 2.3 to 8.0 atm, equivalence ratios from 0.3 to 1.35, and temperatures from 900 to 1300 K using a new experimental flow facility, an aerosol shock tube. The aerosol shock tube combines conventional shock tube methodology with aerosol loading of fuel-oxidizer mixtures. Significant efforts have been made to ensure that the aerosol mixtures were spatially uniform, that the incident shock wave was well-behaved, and that the post-shock conditions and mixture fractions were accurately determined. The nebulizer-generated, narrow, micron-sized aerosol size distribution permitted rapid evaporation of the fuel mixture and enabled separation of the diesel fuel evaporation and diffusion processes that occurred behind the incident shock wave from the chemical ignition processes that occurred behind the higher temperature and pressure reflected shock wave. This rapid evaporation technique enables the study of a wide range of low-vapor-pressure practical fuels and fuel surrogates without the complication of fuel cracking that can occur with heated experimental facilities. These diesel ignition delay measurements extend the temperature and pressure range of earlier flow reactor studies, provide evidence for NTC behavior in diesel fuel ignition delay times at lower temperatures, and provide an accurate data base for the development and comparison of kinetic mechanisms for diesel fuel and surrogate mixtures. Representative comparisons with several single-component diesel surrogate models are also given.     

Direct numerical simulations of turbulent flame expansion in fine sprays

Direct Numerical Simulations of expanding flame kernels following localized ignition in decaying turbulence with the fuel in the form of a fine mist have been performed to identify the effects of the spray parameters on the possibility of self-sustained combustion. Simulations show that the flame kernel may quench due to fuel starvation in the gaseous phase if the droplets are large or if their number is insufficient to result in significant heat release to allow for self-sustained flame propagation for the given turbulent environment. The reaction proceeds in a large range of equivalence ratios due to the random location of the droplets relative to the igniter location that causes a wide range of mixture fractions to develop through pre-evaporation in the unreacted gas and through evaporation in the preheat zone of the propagating flame. The resulting flame exhibits both premixed and non-premixed characteristics.