Numerical study of transient evolution of lifted jet flames: partially premixed flame propagation and influence of physical dimensions

Three-dimensional (3D) unsteady Reynolds-averaged Navier–Stokes simulations of a spark-ignited turbulent methane/air jet flame evolving from ignition to stabilisation are conducted for different jet velocities. A partially premixed combustion model is used involving a correlated joint probability density function and both premixed and non-premixed combustion mode contributions. The 3D simulation results for the temporal evolution of the flame's leading edge are compared with previous two-dimensional (2D) results and experimental data. The comparison shows that the final stabilised flame lift-off height is well predicted by both 2D and 3D computations. However, the transient evolution of the flame's leading edge computed from 3D simulation agrees reasonably well with experiment, whereas evident discrepancies were found in the previous 2D study. This difference suggests that the third physical dimension plays an important role during the flame transient evolution process. The flame brush's leading edge displacement speed resulting from reaction, normal and tangential diffusion processes are studied at different typical stages after ignition in order to understand the effect of the third physical dimension further. Substantial differences are found for the reaction and normal diffusion components between 2D and 3D simulations especially in the initial propagation stage. The evolution of reaction progress variable scalar gradients and its interaction with the flow and mixing field in the 3D physical space have an important effect on the flame's leading edge propagation.     

A comparative study of H2-air premixed flame in micro combustors with different physical and boundary conditions

A numerical study of H₂-air premixed combustion in the micro channels with a detailed chemical reaction mechanism is performed by solving the two-dimensional fully elliptic governing equations of continuity, momentum, energy and species, coupled with the energy equation in the solid wall. A reference case is defined as the combustion in a cylindrical tube with 0.8 mm inner diameter and 8 mm length with a non-slip wall and a uniform velocity profile at the inlet plane. Different physical and boundary conditions have been applied in order to investigate their respective effects on the flame temperature. The conditions studied in the current paper include the combustor size and geometry, inlet velocity profile, axial heat conduction in the solid wall and slip-wall and temperature jump at the gas–solid interface. It is noted that effects of Knudsen number (slip-wall and temperature jump) on the thermal and fluid field are not very significant in a d = 0.4 mm micro combustor. Furthermore, the qualitative effects of Knudsen number on the flame temperature are analysed. The results of this paper indicate that these various boundary and physical conditions have effects on the flame temperature to different extent and should be carefully monitored when applied for different applications.     

A theoretical study of premixed turbulent flame development

Flame development in a statistically stationary and uniform, planar, one-dimensional turbulent flow is theoretically studied. A generalized balance equation for the mean combustion progress variable, which includes turbulent diffusion and pressure-driven transport terms, as well as the mean rate of product creation, is introduced and analyzed by invoking the sole assumption of a self-similar flame structure, well-supported by numerous experiments. The assumption offers the opportunity to simplify the problem by splitting the aforementioned partial differential equation into two ordinary differential equations, which separately model spatial variations of the progress variable and time variations of flame speed and thickness. The self-similar profile of the progress variable, obtained in numerous experiments, is theoretically predicted. Closures of the normalized pressure-driven transport term and mean rate of product creation are obtained. The closed balance equation shows that turbulent diffusion dominates during the initial stage of flame development, followed by the transition to counter-gradient transport in a sufficiently developed flame. A criterion of the transition is derived. The transition is promoted by the heat release and pressure-driven transport. Fully developed mean flame brush thickness and speed are shown to decrease when either density ratio or pressure-driven transport increases. Solutions for the development of the thickness are obtained. The development is accelerated by the pressure-driven transport and heat release.     

Experimental and numerical study of a conical turbulent partially premixed flame

The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH₂O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH₂O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH₂O structure is somewhat thinner than the experiments.     

突扩燃烧室湍流预混反应流的数值模拟

对突扩燃烧室这一典型工程燃烧装置内的湍流预混反应流进行了数值模拟。时平均控制方程组的封闭采用k-ε湍流输运模型和EBU-Arhenius湍流反应模型。模拟结果给出了突扩燃烧室内湍流预混反应流的气体时均流场、组分浓度场与温度场的分布。通过数值模拟结果与实验的比较对EBU-Arhenius模型进行了讨论与评价。     

An experimental and kinetic modeling study of a premixed nitromethane flame at low pressure

A premixed nitromethane/oxygen/argon flame at low pressure (4.67 kPa) has been investigated using tunable vacuum ultraviolet (VUV) photoionization and molecular-beam mass spectrometry. About 30 flame species including hydrocarbons, oxygenated and nitrogenous intermediates have been identified by measurements of photoionization efficiency spectra. Mole fraction profiles of the flame species have been determined by scanning burner position at some selected photon energies. The results indicate that N₂ and NO are the major nitrogenous products in the nitromethane flame. Compared with previous studies on nitromethane combustion, a number of unreported intermediates, including C₃H₄, C₄H₆, C₄H₈, C₂H₂O, C₂H₄O, CH₃CN, H₂CNHO, C₃H₃N and C₃H₇N, are observed in this work. Based on our experimental results and previous modeling studies, a detailed oxidation mechanism including 69 species and 314 reactions has been developed to simulate the flame structure. Despite some small discrepancies, the predictions by the modeling study are in reasonable agreement with the experimental results.     

Hydrodynamic instability of premixed flame propagating in narrow planar channel in the presence of gas flow

The paper analyses the hydrodynamic instability of a flame propagating in the space between two parallel plates in the presence of gas flow. The linear analysis was performed in the framework of a two-dimensional model that describes the averaged gas flow in the space between the plates and the perturbations development of two-dimensional combustion wave. The model includes the parametric dependences of the flame front propagation velocity on its local curvature and on the combustible gas velocity averaged along the height of the channel. It is assumed that the viscous gas flow changes the surface area of the flame front and thereby affects the propagation velocity of the two-dimensional combustion wave. In the absence of the influence of the channel walls on the gas flow, the model transforms into the Darrieus–Landau model of flame hydrodynamic instability. The dependences of the instability growth rate on the wave vector of disturbances, the velocity of the unperturbed gas flow, the viscous friction coefficients and other parameters of the problem are obtained. It is shown that the viscous gas flow in the channel can lead, in some cases, to a significant increase in instability compared with a flame propagating in free space. In particular, the instability increment depends on the direction of the gas flow with respect direction of the flame propagation. In the case when the gas flow moves in the opposite direction to the direction of the flame propagation, the pulsating instability can appear.     

Concept and combustion characteristics of ultra-micro combustors with premixed flame

Under micro-scale combustion influenced by quenching distance, high heat loss, shortened diffusion characteristic time, and flow laminarization, we clarified the most important issues for the combustor of ultra-micro gas turbines (UMGT), such as high space heating rate, low pressure loss, and premixed combustion. The stability behavior of single flames stabilized on top of micro tubes was examined using premixtures of air with hydrogen, methane, and propane to understand the basic combustion behavior of micro premixed flames. When micro tube inner diameters were smaller than 0.4 mm, all of the fuels exhibited critical equivalence ratios in fuel-rich regions, below which no flame formed, and above which the two stability limits of blow-off and extinction appeared at a certain equivalence ratio. The extinction limit for very fuel-rich premixtures was due to heat loss to the surrounding air and the tube. The extinction limit for more diluted fuel-rich premixtures was due to leakage of unburned fuel under the flame base. This clarification and the results of micro flame analysis led to a flat-flame burning method. For hydrogen, a prototype of a flat-flame ultra-micro combustor with a volume of 0.067 cm³ was made and tested. The flame stability region satisfied the optimum operation region of the UMGT with a 16 W output. The temperatures in the combustion chamber were sufficiently high, and the combustion efficiency achieved was more than 99.2%. For methane, the effects on flame stability of an upper wall in the combustion chamber were examined. The results can be explained by the heat loss and flame stretch.     

Influence of water mist on propagation and suppression of laminar premixed flame

The combustion of premixed gas mixtures containing micro droplets of water was studied using one-dimensional approximation. The dependencies of the burning velocity and flammability limits on the initial conditions and on the properties of liquid droplets were analyzed. Effects of droplet size and concentration of added liquid were studied. It was demonstrated that the droplets with smaller diameters are more effective in reducing the flame velocity. For droplets vaporizing in the reaction zone, the burning velocity is independent of droplet size, and it depends only on the concentration of added liquid. With further increase of the droplet diameter the droplets are passing through the reaction zone with completion of vaporization in the combustion products. It was demonstrated that for droplets above a certain size there are two stable stationary modes of flame propagation with transition of hysteresis type. The critical conditions of the transition are due to the appearance of the temperature maximum at the flame front and the temperature gradient with heat losses from the reaction zone to the products, as a result of droplet vaporization passing through the reaction zone. The critical conditions are similar to the critical conditions of the classical flammability limits of flame with the thermal mechanism of flame propagation. The maximum decrease in the burning velocity and decrease in the combustion temperature at the critical turning point corresponds to predictions of the classical theories of flammability limits of Zel'dovich and Spalding. The stability analysis of stationary modes of flame propagation in the presence of water mist showed the lack of oscillatory processes in the frames of the assumed model.     

The unstable behavior of cellular premixed flames induced by intrinsic instability

The unstable behavior of cellular premixed flames induced by intrinsic instability is studied by two-dimensional unsteady calculations of reactive flows. In the present numerical simulation, the compressible Navier–Stokes equation including a one-step irreversible chemical reaction is employed. We consider two basic types of phenomena to account for the intrinsic instability of premixed flames, i.e., hydrodynamic and diffusive-thermal effects. The hydrodynamic effect is caused by the thermal expansion through the flame front; the diffusive-thermal effect is caused by the preferential diffusion of mass versus heat. A disturbance with several wavelength components is superimposed on a planar flame, and the formation of a cellular flame induced by hydrodynamic and diffusive-thermal effects is numerically simulated. After the cellular-flame formation, the combination and division of cells are observed. The behavior of cellular-flame fronts becomes more unstable when the Lewis number is lower than unity, since the diffusive-thermal effect has a great influence on the unstable behavior. The cell size changes with time, and its average is greater than the critical wavelength and becomes smaller by decreasing the Lewis number. The flame velocity of cellular flames depends strongly on the length of computational domain in the direction tangential to the flame front. As the length of computational domain increases, the flame velocity becomes larger. This is because the long-wavelength components of disturbances play an important role in the shape of cellular flames, i.e., in the flame-surface area.     

The role of duct thickness on the quenching process of premixed flame propagation

The propagation of premixed laminar flame in ducts of circular cross-section considering a thermal-diffusive model is investigated numerically. Heat losses by conduction to the channels walls are taken into account using the thermally thin wall regime. The effects and the relationship between thickness and diameter of the tube with the flame speed propagation are studied and the quenching condition is obtained as a function of the heat-loss parameter. The mathematical model employs the axisymmetric energy and species equations. The calculations are based on a two-step chemistry, with an Arrhenius, energetically neutral, radical production reaction followed by an exothermic radical recombination reaction. For large values of the heat-loss parameter, the wall temperature is close to the free stream temperature and all the heat losses through the wall are convected away. No heat feedback occurs. On the other hand, for small values of the heat-loss parameter, a feedback mechanism occurs by transferring heat from the burned gas to the fresh mixture along the tube wall. For values of the heat-loss parameter of order unity, the heat feedback mechanism is able to sustain the flame propagation and the quenching condition disappears, producing an almost planar flame front as the propagation velocity reduces. For this two-step reaction mechanism, the radical species behaviour at the duct walls seems to have negligible effect on the quenching process.     

A multi-zone chemistry mapping approach for direct numerical simulation of auto-ignition and flame propagation in a constant volume enclosure

A direct numerical simulation (DNS) coupling with multi-zone chemistry mapping (MZCM) is presented to simulate flame propagation and auto-ignition in premixed fuel/air mixtures. In the MZCM approach, the physical domain is mapped into a low-dimensional phase domain with a few thermodynamic variables as the independent variables. The approach is based on the fractional step method, in which the flow and transport are solved in the flow time steps whereas the integration of the chemical reaction rates and heat release rate is performed in much finer time steps to accommodate the small time scales in the chemical reactions. It is shown that for premixed mixtures, two independent variables can be sufficient to construct the phase space to achieve a satisfactory mapping. The two variables can be the temperature of the mixture and the specific element mass ratio of H atom for fuels containing hydrogen atoms. An aliasing error in the MZCM is investigated. It is shown that if the element mass ratio is based on the element involved in the most diffusive molecules, the aliasing error of the model can approach zero when the grid in the phase space is refined. The results of DNS coupled with MZCM (DNS-MZCM) are compared with full DNS that integrates the chemical reaction rates and heat release rate directly in physical space. Application of the MZCM to different mixtures of fuel and air is presented to demonstrate the performance of the method for combustion processes with different complexity in the chemical kinetics, transport and flame–turbulence interaction. Good agreement between the results from DNS and DNS-MZCM is obtained for different fuel/air mixtures, including H₂/air, CO/H₂/air and methane/air, while the computational time is reduced by nearly 70%. It is shown that the MZCM model can properly address important phenomena such as differential diffusion, local extinction and re-ignition in premixed combustion.     

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.     

An investigation of flame zones and burning velocities of laminar unconfined methane-oxygen premixed flames

Numerical and experimental investigations of unconfined methane-oxygen laminar premixed flames are presented. In a lab-scale burner, premixed flame experiments have been conducted using pure methane and pure oxygen mixtures having different equivalence ratios. Digital photographs of the flames have been captured and the radial temperature profiles at different axial locations have been measured using a thermocouple. Numerical simulations have been carried out with a C2 chemical mechanism having 25 species and 121 reactions and with an optically thin radiation sub-model. The numerical results are validated against the experimental and numerical results for methane-air premixed flames reported in literature. Further, the numerical results are validated against the results from the present methane-oxygen flame experiments. Visible regions in digital flame photographs have been compared with OH isopleths predicted by the numerical model. Parametric studies have been carried out for a range of equivalence ratios, varying from 0.24 to 1.55. The contours of OH, temperature and mass fractions of product species such as CO, CO₂ and H₂O, are presented and discussed for various cases. By using the net methane consumption rate, an estimate of the laminar flame speed has been obtained as a function of equivalence ratio.     

Investigation of flame structure and burning intensity of partially premixed methane enrichment of syngas using OH-PLIF and kinetic simulation

Various experiments were conducted to study the combustion characteristics of partially premixed methane enrichment of syngas by using the OH-PLIF technique. Experiments were conducted on a co-flow burner, and the methane concentration (X_CH4 = CH₄/(H₂+CO+CH₄)) was varied from 0 to 20%, the overall equivalence ratio was varied from 0.4 to 1.2 and the inner equivalence ratio was varied from 1.5 to 3.5. Kinetic simulation was conducted by using OPPDIF module of CHEMKIN-Pro software. Results show that an increase in X_CH4 and ?_overall weakens the OH signal intensity. Adding methane into the fuel greatly increases the height of the inner flame front, and the increase of methane concentration has a negative effect on flame propagation speed. Meanwhile, simulation results remain consistent with the experiments. The main OH radical production reaction changes from R46: H+HO₂ = 2OH to R38: H+O₂ = O+OH when methane concentration contained in the fuel mixture increases. Sensitivity analysis also indicates that reaction which plays a dominant effect on temperature changes with the increase of methane concentration.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

Thermal expansion effect on the propagation of premixed flames in narrow channels of circular cross-section: Multiplicity of solutions,axisymmetry and non-axisymmetry

The present study examines, in presence of thermal expansion effects, the existence of the multiplicity of solutions previously reported within the context of diffusive-thermal modeling in [15], for lean premixed flames with low Lewis number (Le?<?1) propagating in narrow circular adiabatic channels subject to a Poiseuille flow. For this, direct numerical simulations have been carried out within the framework of variable-density Navier–Stokes equations and zero-Mach-number approximation. The simulations, conducted for both axisymmetric and three-dimensional cylindrical geometries, confirm the coexistence of multiple steady flame structures for a given flow rate. They show that axisymmetric flames concave towards the upstream are more unstable to three-dimensional perturbations than convex (toward the upstream) flames. This result evinces earlier findings obtained from stability analysis. The non-axisymmetry property of the flame is also found to push back the critical flashback limits at larger flow rate when compared to those predicted under the assumption of flame axisymmetry.     

湍流分层燃烧的直接数值模拟和小火焰模型研究

湍流分层燃烧广泛应用于工业燃烧装置,但是目前还比较缺乏适用于湍流分层燃烧的高精度数值模型。本文利用直接数值模拟数据库,对高Karlovitz数分层射流火焰的小火焰模型表现进行了先验性评估。考虑了两种小火焰模型,一种是基于自由传播层流预混火焰的小火焰模型M1,另一种是基于分层对冲小火焰的小火焰模型M2。研究发现M1和M2在c-Z空间的结果与直接数值模拟在定性上是一致的。在物理空间,M2对过程变量反应速率脉动值的预测结果要优于M1.     

Mechanism of magnetic field effect on OH density distribution in a methane–air premixed jet flame

The mechanism of magnetic field effect on OH density distribution in a methane-air premixed jet flame was investigated by means of PLIF measurement and numerical simulation. In the experiment, magnetically induced change in the OH density profile in the flame in a N₂ atmosphere was much smaller than that in air (mixture of 80% N₂ and 20% O₂), and such a phenomenon was qualitatively reproduced by solving the equations for reactive gas dynamics and magnetism in the numerical simulation. Here, N₂ is diamagnetic and little affected by the magnetic field, while O₂ is paramagnetic and influenced due to the magnetic field. This provided the experimental and numerical verification for the mechanism of the magnetic field effect suggested in our previous study. That is, the magnetic force does not directly and selectively induce the change in the diffusion velocity of OH itself. Alternatively, the magnetic force acting on O₂ in the surrounding air, whose mass density and magnetic susceptibility are much larger than that of other chemical species in the flame, causes the change in the convection velocity of the gas mixture and displaces the OH density distribution indirectly and passively. In other words, the most important cause of the OH density change is not the diffusion of OH, but the convection of air containing a large amount of O₂. Furthermore, by careful examination of the magnetic field effect on the flame in the N₂ atmosphere, it was found out that the magnetic force does not only act on O₂ in the surrounding air, but also on O₂ in the premixed gas injected from the burner.     

Computational study of lean premixed turbulent flames using RANSPDF and LESPDF methods

A computational study is performed on a series of four piloted, lean, premixed turbulent jet flames. These flames use the Sydney Piloted Premixed Jet Burner (PPJB), and with jet velocities of 50, 100, 150 and 200 m/s are denoted PM150, PM1100, PM1150 and PM1200, respectively. Calculations are performed using the RANSPDF and LESPDF methodologies, with different treatments of molecular diffusion, with detailed chemistry and flamelet-based chemistry modelling, and using different imposed boundary conditions. The sensitivities of the calculations to these different aspects of the modelling are compared and discussed. Comparisons are made to experimental data and to previously-performed calculations. It is found that, given suitable boundary conditions and treatment of molecular diffusion, excellent agreement between the calculations and experimental measurements of the mean and variance fields can be achieved for PM150 and PM1100. The application of a recently developed implementation of molecular diffusion results in a large improvement in the computed variance fields in the LESPDF calculations. The inclusion of differential diffusion in the LESPDF calculations provides insight on the behaviour in the near-field region of the jet, but its effects are found to be confined to this region and to the species CO, OH and H₂. A major discrepancy observed in many previous calculations of these flames is an overprediction of reaction progress in PM1150 and PM1200, and this discrepancy is also observed in the LESPDF calculations; however, a parametric study of the LESPDF mixing model reveals that, with a sufficiently large mixing frequency, calculations of these two flames are capable of yielding improved reaction progress in good qualitative agreement with the mean and RMS scalar measurements up to an x/D of 30. Lastly, the merits of each computational methodology are discussed in light of their computational costs.