Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

Capturing multi-regime combustion in turbulent flames with a virtual chemistry approach

Multiple flame regimes are encountered in industrial combustion chambers, where premixed, stratified and non-premixed flame regions may coexist. To obtain a predictive tool for pollutant formation predictions, chemical flame modeling must take into account the influence of such complex flame structure. The objective of this article is to apply and compare two reduced chemistry models on both laminar and turbulent multi-regime flame configurations in order to analyze their capabilities in predicting flame structure and CO formation. The challenged approaches are (i) a premixed flamelet-based tabulated chemistry method, whose thermochemical variables are parameterized by a mixture fraction and a progress variable, and (ii) a virtual chemical scheme which has been optimized to retrieve the properties of canonical premixed and non-premixed 1-D laminar flames. The methods are first applied to compute a series of laminar partially-premixed methane-air counterflow flames. Results are compared to detailed chemistry simulations. Both approaches reproduced the thermal flame structure but only the virtual chemistry captures the CO formation in all ranges of equivalence ratio from stoichiometry premixed flame to pure non-premixed flame. Finally, the two chemical models combined with the Thickened Flame model for LES are challenged on a piloted turbulent jet flame with inhomogeneous inlet, the Sydney inhomogeneous burner. Mean and RMS of temperature and CO mass fraction radial profiles are compared to available experimental data. Scatter data in mixture fraction space and Wasserstein metric of numerical and experimental data are also studied. The analyses confirm again that the virtual chemistry approach is able to account for the impact of multi-regime turbulent combustion on the CO formation.     

Nonlinear hydrodynamic and thermoacoustic oscillations of a bluff-body stabilised turbulent premixed flame

Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.     

Turbulent flame propagation limits of polymethylmethacrylate particle cloud–ammonia–air co-combustion

Ammonia is a highly promising energy carrier for achieving a carbon-neutral society. The co-combustion of solid particle clouds–ammonia, in particular, is considered an efficient and feasible method of reducing carbon dioxide emissions. Understanding turbulent flame stabilization and extinguishment processes during the two-phase hybrid-mixture co-combustion of solid particle clouds–ammonia is essential for the co-combustion technology to be used in combustors. To the best of our knowledge, this is the first study to describe the turbulent flame propagation limits and associated mechanism on the co-combustion of solid particle clouds–ammonia–air. Turbulent flame propagation experiments on silica particle clouds–ammonia–air mixing combustion and polymethylmethacrylate (PMMA) particle cloud–ammonia–air co-combustion were conducted in this work using a novel fan-stirred constant-volume vessel to clarify the turbulent flame propagation limits and associate mechanism of solid particle cloud–ammonia–air co-combustion. Results showed that adding inert silica particles contracted the turbulent flame propagation limits of premixed ammonia–air mixtures. However, adding PMMA particles expanded and then contracted the turbulent flame propagation limits of a premixed ammonia–air mixture as the ammonia equivalence ratio increased from lean to rich. In the solid particle cloud–ammonia–air co-combustion, reactive particles induce two types of effects on the turbulent flame propagation limits of premixed ammonia–air mixtures: The local equivalence ratio increment effect is caused by adding volatile matter from preheated particles in the preheat zone of the flame front, and the heat sink negative effect is induced by the unburned particles.     

An analytical model for the prediction of the dynamic response of premixed flames stabilized on a heat-conducting perforated plate

The dynamic response of a premixed flame stabilized on a heat-conducting perforated plate depends critically on their coupled thermal interaction. The objective of this paper is to develop an analytical model to capture this coupling. The model predicts the mean flame base standoff distance; the flame base area, curvature and speed; and the burner plate temperature given the operating conditions; the mean velocity, temperature and equivalence ratio of the reactants; thermal conductivity and the perforation ratio of the burner. This coupled model is combined with our flame transfer function (FTF) model to predict the dynamic response of the flame to velocity perturbations. We show that modeling the thermal coupling between the flame and the burner, while accounting for the two-dimensionality of the former, is critical to predicting the dynamic response characteristics such as the overshoot in the gain curve (resonant condition) and the phase delay. Good agreement with the numerical and experimental results is demonstrated over a range of conditions.     

Large eddy simulation of multi-regime combustion with a two-progress variable approach for carbon monoxide

Simulations of two cases in a novel multi-regime burner configuration are undertaken using a presumed joint probability density function (PDF) approach with tabulated chemistry. The flame conditions are varied by changing the central jet equivalence ratio, which produces different multi-regime combustion modes in the non-premixed inner flame. An outer premixed flame and recirculation zone behind a bluff body are present to supply heat and combustion products to stabilise the inner flame. A two-progress variable approach is tested to improve predictions of carbon monoxide (CO) in the post-flame regions, where CO oxidation occurs. The large eddy simulation set-up and sub-grid combustion model are assessed through comparisons with time-averaged measurements for radial profiles at different streamwise locations. The jet break-up length, the shear layers and the mixture fraction distribution are well captured in both cases. The temperature distribution is well captured for the inner flame in each case but the temperature and mixture fraction are over predicted in the downstream regions of the outer premixed flame, which is due to increased dilatation that suppresses air entrainment. Improved predictions of the CO mass fraction are obtained for the outer premixed flames with the two-progress variable approach. Over predictions are seen in the upstream regions of the inner flame when the CO mass fraction is obtained from a look-up table, suggesting that the CO mass fraction should be transported to include the convection/diffusion balance in regions where there is no flame. Furthermore, transporting the CO mass fraction with a one-progress variable approach produces over predictions in the burnt regions, suggesting a two-progress variable model is needed to capture the consumption region of CO. The multi-regime combustion characteristics are observed to be stronger in flame MRB26b, where non-premixed and rich premixed combustion is present. For flame MRB18b, the non-premixed contribution is smaller and weak stratified combustion is observed.     

Measurements of three-dimensional mean flame surface area ratio in turbulent premixed Bunsen flames

A data processing scheme with particular emphasis on proper flame contour smoothing is developed and applied to measure the three-dimensional mean flame surface area ratio in turbulent premixed flames. The scheme is based on the two-sheet imaging technique such that the mean flame surface area ratio is an average within a window covering a finite section of the turbulent flame brush. This is in contrast to the crossed-plane tomograph technique which applies only to a line. Two sets of Bunsen flames have been investigated in this work with the turbulent Reynolds number up to 4000 and the Damköhler number ranging from less than unity to close to 10. The results show that three-dimensional effects are substantial. The measured three-dimensional mean flame surface area ratio correlates well with a formula similar to the Zimont model for turbulent burning velocity but with different model constants. Also, the mean flame surface area ratio displays a weak dependency on turbulence intensity but a strong positive dependency on the turbulence integral length scale.     

Flame front analysis of high-pressure turbulent lean premixed methane–air flames

An experimental study on lean turbulent premixed methane–air flames at high pressure is conducted by using a turbulent Bunsen flame configuration. A single equivalence ratio flame at Φ = 0.6 is explored for pressures ranging from atmospheric pressure to 0.9 MPa. LDA measurements of the cold flow indicate that turbulence intensities and the integral length scale are not sensitive to pressure. Due to the decreased kinematic viscosity with increasing pressure, the turbulent Reynolds numbers increase, and isotropic turbulence scaling relations indicate a large decrease of the smallest turbulence scales. Available experimental results and PREMIX code computations indicate a decrease in laminar flame propagation velocities with increasing pressure, essentially between the atmospheric pressure and 0.5 MPa. The u′/S_L ratio increases therefore accordingly. Instantaneous flame images are obtained by Mie scattering tomography. The images and their analysis show that pressure increase generates small scale flame structures. In an attempt to generalize these results, the variance of the flamelet curvatures, the standard deviation of the flamelet orientation angle, and the flamelet crossing lengths have been plotted against which is proportional to the ratio between the integral and Taylor length scales, and which increases with pressure. These three parameters vary linearly with the ratio between large and small turbulence scales and clearly indicate the strong effect of this parameter on premixed turbulent flame dynamics and structure. An obvious consequence is the increase in flame surface density and hence burning rate with pressure, as confirmed by its direct determination from 2D tomographic images.     

Theoretical and experimental studies on mesoscale flame propagation and extinction

Mesoscale flame propagation and extinction of premixed flames in channels are investigated theoretically and experimentally. Emphasis is placed on the effect of wall heat loss and the wall–flame interaction via heat recirculation. At first, an analytical solution of flame speed in mesoscale channels is obtained. The results showed that channel width, flow velocity, and wall thermal properties have dramatic effects on the flame propagation and lead to multiple flame regimes and extinction limits. With the decrease in channel width, there exist two distinct flame regimes, a fast burning regime and a slow burning regime. The existence of the new flame regime and its extended flammability limit render the classical quenching diameter inapplicable. Furthermore, the results showed that at optimum conditions of flow velocity and wall thermal properties, mesoscale flames can propagate faster than the adiabatic flame. Second, numerical simulation with detailed chemistry demonstrated the existence of multiple flame regimes. The results also showed that there is a non-linear dependence of the flame speed on equivalence ratio. Moreover, it is shown that the Nusselt number has a significant impact on this non-linear dependence. Finally, the non-linear dependence of flame speed on equivalence ratio for both flame regimes is measured using a C₃H₈–air mixture. The results are in good agreement with the theory and numerical simulation.     

Modelling of turbulent lifted jet flames using flamelets: a priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z–c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z–c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences.     

Temperature measurement of axisymmetric partially premixed methane/air flame in a co-annular burner using Mach–Zehnder interferometry

In this paper partially premixed laminar methane/air co-flow flame is studied experimentally. Methane–air flame is established on an axisymmetric co-annular burner. The fuel-air jet flows from the central tube while the secondary air flows from the region between the inner and the outer tube. The aim is to investigate the flame characteristics for methane/air axisymmetric partially premixed flame using Mach–Zehnder interferometry. Different equivalence ratios (φ=1.4–2.2) and Reynolds numbers (Re=100–1200) are considered in the study. Flame generic visible appearance and the corresponding fringe map structures are also investigated. It is seen that the fringe maps are poorly influenced by equivalence ratio variations at constant Reynolds number but are significantly affected by Reynolds number variations in constant equivalence ratio. Temperatures obtained from optical techniques are compared with those obtained from thermocouples and good agreement is observed. It is concluded that the effect of Reynolds number increment on maximum flame temperature is negligible while equivalence ratio reduction increases maximum flame temperature substantially.     

Extinction and flame bifurcations of stretched dimethyl ether premixed flames

Extinction limits and flame bifurcation of lean premixed dimethyl ether–air flames are numerically investigated using the counterflow flame with a reduced chemistry. Emphasis is paid to the combined effect of radiation and flame stretch on the extinction and flammability limits. A method based on the reaction front is presented to predict the Markstein length. The predicted positive Markstein length agrees well with the experimental data. The results show that flow stretch significantly reduces the flame speed and narrows the flammability limit of the stretched dimethyl ether–air flame. It is found that the combined effect of radiation and flow stretch results in a new flame bifurcation and multiple flame regimes. At an equivalence ratio slightly higher than the flammability limit of the planar flame, the distant flame regime appears at low stretch rates. With an increase in the equivalence ratio, in addition to the distant flame, a weak flame isola emerges at moderate stretch rates. With a further increase in the equivalence ratio, the distant flame and the weak flame branches merge together, resulting in the splitting of the weak flame branch into two weak flame branches, one at low stretch and the other at high stretch. Flame stability analysis demonstrates that the high stretch weak flame is also stable. Furthermore, a K-shaped flammability limit diagram showing various flame regimes and their extinction limits is obtained.     

Characterization of premixed swirling methane/air diffusion flame through filtered Rayleigh scattering

Characteristics of a premixed, swirl methane/air diffusion flame at atmospheric pressure are measured by filtered Rayleigh scattering (FRS). Three operating conditions are investigated with the equivalence ratios of the methane/air flame covering a range of 0.67—0.83. Under each condition, single-shot and averaged FRS images over a region measured 39.3×65.6 mm² at seven cross sections of the flame are collected to demonstrate the flame behavior. A gradient calculation algorithm is applied to identify reaction zone locations and structures in the instantaneous FRS measurements. Statistical analysis for the mean FRS measurements is performed by means of joint probability density functions. The experimental results indicate that thermochemical state of the swirl flame is strongly influenced by equivalence ratio, leading to varieties of flame structures and temperature distributions. The gradient of the instantaneous FRS images clearly illustrates the characteristics of the reaction zone. The results also demonstrate that FRS can provide detailed insights into the behavior of turbulent flames.     

Combustion of partially premixed spray jets

Gas turbines, liquid rocket motors, and oil-fired furnaces utilize the spray combustion of continuously injected liquid fuels. In most cases, the liquid spray is mixed with an oxidizer prior to combustion, and further oxidizer is supplied from the outside of the spray to complete diffusion combustion. This rich premixed spray is called “partially premixed spray.” Partially premixed sprays have not been studied systematically although they are of practical importance. In the present study, the burning behavior of partially premixed sprays was experimentally studied with a newly developed spray burner. A fuel spray and an oxidizer, diluted with nitrogen, was injected into the air. The overall equivalence ratio of the spray jet was set larger than unity to establish partially premixed spray combustion. In the present burner, the mean droplet diameter of the atomized liquid fuel could be varied without varying the overall equivalence ratio of the spray jet. Two combustion modes with and without an internal flame were observed. As the mean droplet diameter was increased or the overall equivalence ratio of the spray jet was decreased, the transition from spray combustion only with an external group flame to that with the internal premixed flame occurred. The results suggest that the internal flame was supported by flammable mixture through the vaporization of fine droplets, and the passage of droplet clusters deformed the internal flame and caused internal flame oscillation. The existence of the internal premixed flame enhanced the vaporization of droplets in the post-premixed-flame zone within the external diffusion flame.     

预混火焰传递函数的测量与分析

火焰传递函数是理解和控制振荡燃烧的理论基础.本文通过自发化学荧光法测量放热率,双传声器技术获得燃烧器出口的压力和速度脉动,研究了不同雷诺数、当量比和扰动频率下的预混火焰传递函数.结果显示,随着脉动幅值的增加,火焰传递函数将出现由线性到非线性的变化过程.当量比、雷诺数和扰动频率的改变,都会影响火焰传递函数的幅值和相位特性...     

Laser-induced plasma in methane and dimethyl ether for flame ignition and combustion diagnostics

In this paper we report the investigation of the laser-induced breakdown and ignition behaviour of methane/air and dimethyl ether (DME)/air mixtures. Moreover, the optical emission from the induced plasma is utilized for determining the mixture composition quantitatively by means of laser-induced breakdown spectroscopy (LIBS). To the best of the authors’ knowledge, LIBS and laser ignition of DME have not been reported in literature before. The technique under investigation is finally employed for combustion diagnostics in laminar as well as turbulent flames. In the laminar premixed and non-premixed flames the LIBS spectra allow spatially resolved measurements of the equivalence ratio and enable studying the mixing of gases provided through the burner with the surrounding room air. In addition, the breakdown threshold of the applied laser pulse energy yields an estimate for the local temperature. In the turbulent cases single-shot LIBS spectra are recorded at fixed position allowing the derivation of local statistical fluctuations of the equivalence ratio in partially premixed jet flames. The results show that laser-induced breakdowns have a strong potential for flame diagnostics and, under suitable conditions, for the ignition of combustible mixtures.     

Temperature measurement of a premixed radially symmetric methane flame jet using the Mach-Zehnder Interferometry

The temperature field of a premixed methane symmetric laminar flame jet is visualized by studying the interferograms of the flame, using the Mach-Zehnder Interferometry. Two kinds of oxidizers are chosen for combustion: industrially pure oxygen and oxygen-enriched air. The flame is chosen to be both lean, and rich. For the lean oxygen-enriched combustion (OEC), the equivalence ratio was held constant at 0.5, and the oxygen enrichment was adjusted to 0.5 and 0.6, and for rich OEC, equivalence ratio is chosen to be 1.2 while the oxygen enrichment was 0.7 and 0.8. For methane/oxygen combustion, the equivalence ratio varied from 0.35 to 0.55 for the lean flame, and 1.3 and 1.7 for the rich flame. Attempt was made to keep the Reynolds number unchanged at 500, for OEC, and 1000, for methane/oxygen flame. In the present study a non-contact method is successfully developed to measure the temperature field of a premixed radially symmetric laminar methane flame jet. The effect of oxygen enrichment and equivalence ratio on temperature field is also investigated and depicted.     

Direct numerical simulations of rich premixed turbulent n-dodecane/air flames at diesel engine conditions

Rich premixed turbulent n-dodecane/air flames at diesel engine conditions are analyzed using direct numerical simulations. The conditions correspond to a parametric variation of the Engine Combustion Network Spray A (pressure 60 atm; oxidizer oxygen level and temperature 21% and 900 K, respectively; fuel temperature 363 K). Three simulations with equivalence ratios of 3, 5, and 7 are performed with a Karlovitz number (Ka, based on flame time) of order 100 to match the estimated Ka of the rich premixed combustion region in Spray A. At these conditions, the reference laminar flames exhibit a complex structure which involves both low-temperature chemistry (LTC) and high-temperature chemistry over a wide range of length scales. In the presence of turbulence, the flame structure is strongly affected in physical space and the reaction zone exhibits a very complex structure in which broken, distributed, and thin regions co-exist, especially for the leanest case. However, the contribution of the LTC pathway is only weakly affected by turbulence. In progress variable space, the mean flame structure, including the chemical source terms, is found to match remarkably well that of the corresponding unity Lewis number laminar flame, particularly for the $?=$ 3 and 5 cases. This behavior is attributed to the strong turbulent mixing occurring throughout the flames/reaction zones, which suppresses differential diffusion effects. Nevertheless, large conditional fluctuations around the mean chemical source terms are identified. These are found to correlate very well with radical species mass fractions such as OH. In addition, a similar functional dependence is obtained from counterflow laminar flames. As such, it appears from these results that laminar flame models have a potential to be used to represent the thermochemical state of rich premixed turbulent flames under diesel engine conditions.     

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.