Tabulated chemistry approaches for laminar flames: Evaluation of flame-prolongation of ILDM and flamelet methods

The present study considers the performance of tabulation methods for numerical simulation of complex chemical kinetics in laminar combusting flows and compares their predictions to results obtained by direct calculation. Two tabulation methods are considered: the Flame Prolongation of Intrinsic low-dimensional manifold (FPI) method and Steady Laminar Flamelet Model (SLFM). The FPI method is of current interest as it is a potentially unifying approach capable of dealing with both premixed and non-premixed flames for gaseous fuels. SLFM tabulation methods are popular for non-premixed flames and form a good basis for comparing the performance of the FPI approach. The performance of each method is also evaluated by comparing the results to the direct simulation of the laminar flames using two chemical kinetic schemes: simplified chemistry involving five species and one reaction and detailed chemistry involving 53 species and 325 reaction steps. As part of the evaluation process, the computational cost of each method is also assessed. The laminar flames considered in this study include: freely propagating laminar premixed flames, a two-dimensional axisymmetric methane–air opposed-jet diffusion flame, and a two-dimensional axisymmetric methane–air co-flow diffusion flame. Both tabulation methods are implemented in a parallel adaptive mesh refinement (AMR) framework for solving the complete set of governing partial differential equations. These equations are solved using a fully-coupled finite-volume formulation on body-fitted multi-block quadrilateral mesh. Significant improvements in terms of reduced computational requirements, as measured by both storage and processing time, are demonstrated for the tabulated methods.     

Combined experimental and computational study of laminar, axisymmetric hydrogen–air diffusion flames

We investigate the structure of two-dimensional, axisymmetric, laminar hydrogen–air flames in which a cylindrical fuel stream is surrounded by coflowing air, using laser-diagnostic and computational methods. Spontaneous Raman scattering and coherent anti-Stokes Raman scattering (CARS) are used to measure the distributions of major species and temperature. Computationally, we solve the governing conservation equations for mass, momentum, energy, and species, using detailed chemistry and transport. The fuel is diluted with nitrogen (1:1) to reduce heat transfer to the burner, to match the zero temperature gradient at the fuel exit. Three average fuel exit velocities are studied: 18, 27, and 50 cm/s. Comparisons of the measured and computed results are performed for radial profiles at a number of axial positions, and along the axial centerline. Peak major species mole fractions and temperatures are quantitatively predicted by the computations, and the axial species profiles are predicted to within the experimental uncertainty. In the radial profiles studied, base-case computations excluding thermal diffusion of light species were in excellent agreement with the measurements. While the addition of thermal diffusion led to some discrepancy with the measured results, the magnitude of the differences was no more than 25%. The computations predicted the axial centerline profiles from the burner exit to the maximum temperature well, though the experimental temperatures in the downstream mixing region decreased somewhat faster than the computed profiles. Radiative losses are seen to be negligible in these flames, and changes in transport properties and variations in initial flow velocities generally led to only modest changes in the axial profiles. The results also show that the detailed axial profiles of major species and temperature at different fuel jet velocities scale quantitatively with the jet velocity.     

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.     

Numerical study of the effects of gravity on soot formation in laminar coflow methane/air diffusion flames under different air stream velocities

Numerical simulations of laminar coflow methane/air diffusion flames at atmospheric pressure and different gravity levels were conducted to gain a better understanding of the effects of gravity on soot formation by using relatively detailed gas-phase chemistry and complex thermal and transport properties coupled with a semi-empirical two-equation soot model. Thermal radiation was calculated using the discrete-ordinates method coupled with a non-grey model for the radiative properties of CO, CO₂, H₂O, and soot. Calculations were conducted for three coflow air velocities of 77.6, 30, and 5 cm/s to investigate how the coflowing air velocity affects the flame structure and soot formation at different levels of gravity. The coflow air velocity has a rather significant effect on the streamwise velocity and the fluid parcel residence time, especially at reduced gravity levels. The flame height and the visible flame height in general increase with decreasing the gravity level. The peak flame temperature decreases with decreasing either the coflow air stream velocity or the gravity level. The peak soot volume fraction of the flame at microgravity can either be greater or less than that of its normal gravity counterpart, depending on the coflow air velocity. At sufficiently high coflow air velocity, the peak soot volume fraction increases with decreasing the gravity level. When the coflow air velocity is low enough, soot formation is greatly suppressed at microgravity and extinguishment occurs in the upper portion of the flame with soot emission from the tip of the flame owing to incomplete oxidation. The numerical results provide further insights into the intimate coupling between flame size, residence time, thermal radiation, and soot formation at reduced gravity level. The importance of thermal radiation heat transfer and coflow air velocity to the flame structure and soot formation at microgravity is demonstrated for the first time.     

Head-on quenching of laminar premixed methane flames diluted with hot combustion products

Transient head-on quenching of laminar premixed methane flames diluted with hot combustion products is analyzed using full-chemistry 1D DNS. The impact of the dilution level, pressure and wall temperature on carbon monoxide (CO) emissions is investigated. Increasing dilution level and pressure reduce peak average near-wall CO concentrations, and reduce the near-wall CO reduction rate. However, the peak average near-wall CO and near-wall CO reduction rate increase with increasing wall temperature. Analysis of the species transport budget for CO near the wall before, during and after quenching indicates that there are conditions where diffusion is the dominant transport term. As a consequence, it may be possible to model the near-wall CO using only the integrated diffusion term within certain spatial regions. Dilution increases the size of these regions, whereas increasing pressure reduces this size.     

A computational framework for predicting laminar reactive flows with soot formation

Numerical modeling is an attractive option for cost-effective development of new high-efficiency, soot-free combustion devices. However, the inherent complexities of hydrocarbon combustion require that combustion models rely heavily on engineering approximations to remain computationally tractable. More efficient numerical algorithms for reacting flows are needed so that more realistic physics models can be used to provide quantitative soot predictions. A new, highly-scalable combustion modeling tool has been developed specifically for use on large multiprocessor computer architectures. The tool is capable of capturing complex processes such as detailed chemistry, molecular transport, radiation, and soot formation/destruction in laminar diffusion flames. The proposed algorithm represents the current state of the art in combustion modeling, making use of a second-order accurate finite-volume scheme and a parallel adaptive mesh refinement (AMR) algorithm on body-fitted, multiblock meshes. Radiation is modeled using the discrete ordinates method (DOM) to solve the radiative transfer equation and the statistical narrow-band correlated-k (SNBCK) method to quantify gas band absorption. At present, a semi-empirical model is used to predict the nucleation, growth, and oxidation of soot particles. The framework is applied to two laminar coflow diffusion flames which were previously studied numerically and experimentally. Both a weakly-sooting methane–air flame and a heavily-sooting ethylene–air flame are considered for validation purposes. Numerical predictions for these flames are verified with published experimental results and the parallel performance of the algorithm analyzed. The effects of grid resolution and gas-phase reaction mechanism on the overall flame solutions were also assessed. Reasonable agreement with experimental measurements was obtained for both flames for predictions of flame height, temperature and soot volume fraction. Overall, the algorithm displayed excellent strong scaling performance by achieving a parallel efficiency of 70% on 384 processors. The proposed algorithm proved to be a robust, highly-scalable solution method for sooting laminar flames.     

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.     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

Computational and experimental study of JP-8, a surrogate, and its components in counterflow diffusion flames

Non-sooting counterflow diffusion flames have been studied both computationally and experimentally, using either JP-8, or a six-component JP-8 surrogate mixture, or its individual components. The computational study employs a counterflow diffusion flame model, the solution of which is coupled with arc length continuation to examine a wide variety of inlet conditions and to calculate extinction limits. The surrogate model includes a semi-detailed kinetic mechanism composed of 221 gaseous species participating in 5032 reactions. Experimentally, counterflow diffusion flames are established, in which multicomponent fuel vaporization is achieved through the use of an ultrasonic nebulizer that introduces small fuel droplets into a heated nitrogen stream, fostering complete vaporization without fractional distillation. Temperature profiles and extinction limits are measured in all flames and compared with predictions using the semi-detailed mechanism. These measurements show good agreement with predictions in single-component n-dodecane, methylcyclohexane, and iso-octane flames. Good agreement also exists between predicted and measured variables in flames of the surrogate, and the agreement is even better between the experimental JP-8 flames and the surrogate predictions.     

A numerical study on extinction behaviour of laminar micro-diffusion flames

We conducted a numerical study on the fluid dynamic, thermal and chemical structures of laminar methane–air micro flames established under quiescent atmospheric conditions. The micro flame is defined as a flame on the order of one millimetre or less established at the exit of a vertically-aligned straight tube. The numerical model consists of convective–diffusive heat and mass transport with a one-step, irreversible, exothermic reaction with selected kinetics constants validated for near-extinction analyses. Calculations conducted under the burner rim temperature 300 K and the adiabatic burner wall showed that there is the minimum burner diameter for the micro flame to exist. The Damköhler number (the ratio of the diffusive transport time to the chemical time) was used to explain why a flame with a height of less than a few hundred microns is not able to exist under the adiabatic burner wall condition. We also conducted scaling analysis to explain the difference in extinction characteristics caused by different burner wall conditions. This study also discussed the difference in governing mechanisms between micro flames and microgravity flames, both of which exhibit similar spherical flame shape.     

A computational study of soot formation and flame structure of coflow laminar methane/air diffusion flames under microgravity and normal gravity

A numerical study is conducted of methane–air coflow diffusion flames at microgravity (μg) and normal gravity (1g), and comparisons are made with experimental data in the literature. The model employed uses a detailed gas phase chemical kinetic mechanism that includes PAH formation and growth, and is coupled to a sectional soot particle dynamics model. The model is able to accurately predict the trends observed experimentally with reduction of gravity without any tuning of the model for different flames. The microgravity sooting flames were found to have lower temperatures and higher volume fraction than their normal gravity counterparts. In the absence of gravity, the flame radii increase due to elimination of buoyance forces and reduction of flow velocity, which is consistent with experimental observations. Soot formation along the wings is seen to be surface growth dominated, while PAH condensation plays a more major role on centreline soot formation. Surface growth and PAH growth increase in microgravity primarily due to increases in the residence time inside the flame. The rate of increase of surface growth is more significant compared to PAH growth, which causes soot distribution to shift from the centreline of the flame to the wings in microgravity.     

Identification and analysis of instability in non-premixed swirling flames using LES

Large eddy simulations (LES) of turbulent non-premixed swirling flames based on the Sydney swirl burner experiments under different flame characteristics are used to uncover the underlying instability modes responsible for the centre jet precession and large scale recirculation zone. The selected flame series known as SMH flames have a fuel mixture of methane-hydrogen (50:50 by volume). The LES solves the governing equations on a structured Cartesian grid using a finite volume method, with turbulence and combustion modelling based on the localised dynamic Smagorinsky model and the steady laminar flamelet model respectively. The LES results are validated against experimental measurements and overall the LES yields good qualitative and quantitative agreement with the experimental observations. Analysis showed that the LES predicted two types of instability modes near fuel jet region and bluff body stabilised recirculation zone region. The mode I instability defined as cyclic precession of a centre jet is identified using the time periodicity of the centre jet in flames SMH1 and SMH2 and the mode II instability defined as cyclic expansion and collapse of the recirculation zone is identified using the time periodicity of the recirculation zone in flame SMH3. Finally frequency spectra obtained from the LES are found to be in good agreement with the experimentally observed precession frequencies.     

Laser-induced fluorescence measurements and modeling of absolute CH concentrations in strained laminar methane/air diffusion flames

We have applied linear laser-induced fluorescence to obtain spatially resolved profiles of CH radicals in laminar methane/air and methane/nitric oxide/air counterflow diffusion flames at atmospheric pressure. Excitation and detection of transitions in the A–X band and calibrating the optical detection efficiency via Rayleigh scattering allowed the determination of absolute radical concentrations. Flames at strain rates from 59 to 269 s⁻¹ were studied to characterize the strain rate dependence of the CH concentration. The work shows that CH concentrations increase with increasing strain rate. Comparisons have been made with predicted CH levels obtained using two different chemical kinetic mechanisms (Lindstedt et al. and GRI-Mech. 3.0). Computed concentrations are shown to be in good agreement with experimental data. It was furthermore found that the addition of up to 600 ppm NO to the fuel did not have a measurable effect on the CH radical concentration. This is also in agreement with predictions from both mechanisms. The current work has shown that measurements of absolute CH radical concentrations are possible in non-premixed flames without the need for spatial temperature or quenching corrections.     

Particle formation in opposed-flow diffusion flames of ethylene: An experimental and numerical study

An experimental and numerical study on particles inception and growth is performed in opposed-flow diffusion flames of ethylene and air characterized by different sooting tendencies. Spectrally resolved UV-visible laser induced fluorescence, laser induced incandescence and laser light scattering measurements are used to characterize different classes of combustion-generated compounds. A detailed kinetic model accounting for both gas-phase and particle formation is used. Comparison between experimental results and numerical predictions gives a qualitative view of the mechanism of particle formation in opposed-flow flames.Particle inception is the result of both chemical growth and coagulation of aromatic compounds. In the region close to the flame front where the temperature is relatively high and radicals are abundant, the particle inception is due to a chemical growth mechanism by which aromatic molecules add aromatic radicals leading to the formation of biphenyl-like structures. The growth process continues as high-molecular mass aromatics are moved away from the flame zone towards the stagnation plane by the addition of acetylene and other aromatics forming particles of increasing sizes. Graphitization of these particles and thermal annealing lead to the formation of soot particles. At relatively lower temperatures, found across the stagnation plane, particles growth still occurs and it is mainly due to a process of physical coagulation of PAHs.The experimental and numerical results obtained in this work demonstrate and explain the sensibility of inception and growth of particles to radical concentration and temperature in opposed-flow flame configurations.     

Cellular instabilities of expanding hydrogen/propane spherical flames at elevated pressures: theory and experiment

An experimental and theoretical investigation of the onset of cellular instabilities on spherically expanding flames in mixtures of hydrogen and propane in air at elevated pressures was conducted. Critical conditions for the onset of instability were measured and mapped out over a range of pressures and mixture compositions. An asymptotic theory of hydrodynamic and diffusional-thermal cell development on flames in mixtures comprised of two scarce fuels burning in air was also formulated. Predicted values of Peclet number, defined as the flame radius at the onset of instability normalized by the flame thickness, were shown to compare favorably with the experimentally measured values.     

Analytical predictions of shapes of laminar diffusion flames in microgravity and earth gravity

Flame shape is an important observed characteristic of flames that can be used to scale flame properties such as heat release rates and radiation. Flame shape is affected by fuel type, oxygen levels in the oxidiser, inverse burning and gravity. The objective of this study is to understand the effect of high oxygen concentrations, inverse burning, and gravity on the predictions of flame shapes. Flame shapes are obtained from recent analytical models and compared with experimental data for a number of inverse and normal ethane flame configurations with varying oxygen concentrations in the oxidiser and under earth gravity and microgravity conditions. The Roper flame shape model was extended to predict the complete flame shapes of laminar gas jet normal and inverse diffusion flames on round burners. The Spalding model was extended to inverse diffusion flames. The results show that the extended Roper model results in reasonable predictions for all microgravity and earth gravity flames except for enhanced oxygen normal diffusion flames under earth gravity conditions. The results also show trends towards cooler flames in microgravity that are in line with past experimental observations. Some key characteristics of the predicted flame shapes and parameters needed to describe the flame shape using the extended Roper model are discussed.     

Flame propagation speed and Markstein length of spherically expanding flames: Assessment of extrapolation and measurement techniques

Laminar burning velocities are of great importance in many combustion models as well as for validation and improvement of chemical kinetic schemes. Determining laminar burning velocities with high accuracy is quite challenging and different approaches exist. Hence, a comparison of existing methods measuring and evaluating laminar burning velocities is of interest. Here, two optical diagnostics, high speed tomography and Schlieren cinematography, are simultaneously set up to investigate methods for evaluating laminar flame speed in a spherical flame configuration. The hypothesis to obtain the same flame propagation radii over time with the two different techniques is addressed. Another important aspect is the estimation of flame properties, such as the unstretched flame propagation speed and Markstein length in the burnt gas phase and if these are estimated satisfactorily by common experimental approaches. Thorough evaluation of the data with several extrapolation techniques is undertaken. A systematic extrapolation approach is presented to give more confidence into results generated experimentally. The significance of the linear extrapolation routine is highlighted in this context. Measurements of spherically expanding flames are carried out in two high-pressure, high-temperature, constant-volume vessels at RWTH in Aachen, Germany and at ICARE in Orleans, France. For the discussion of the systematic extrapolation approach, flame speed measurements of methane / air mixtures with mixture Lewis numbers moderately away from unity are used. Conditions were varied from lean to rich mixtures, at temperatures of 298–373 K, and pressures of 1 atm and 5 bar.     

An experimental and numerical study on the adequacy of CH as a flame marker in premixed methane flames

The CH radical is frequently used as a flame marker because it is relatively short-lived and is present over a narrow region in flames. Discontinuities in the CH field are thus often interpreted as localized extinction of the flame. Recently, however, the adequacy of CH laser-induced fluorescence (LIF) as a flame marker was questioned by an experimental study of flame–vortex interactions in highly N₂-diluted premixed methane flames. We demonstrate both experimentally and numerically that anomalies in the transient response of CH in this earlier study were due to reactant composition variations in the vortex. In addition, we evaluate the adequacy of CH LIF as a flame marker over a much broader range of conditions. Previous numerical studies showed that heat release rate correlates reasonably well with peak [HCO] and the concentration product [OH][CH₂O], but poorly with [CH], in highly N₂-diluted premixed methane flames. Here, the correlation between heat release rate and CH is investigated both experimentally, by performing simultaneous measurements of CH, OH, and CH₂O LIF, and numerically. We consider undiluted and N₂-diluted premixed methane flames over a range of strain rates and stoichiometries. Results are reported for flames subjected to unsteady stretch and reactant composition variations. For all N₂-dilution levels considered, the peak CH LIF signal correlates poorly with heat release rate when the stoichiometry of the reactant mixture changes from rich to lean. However, when flames are subjected to stretch, the correlation between CH and heat release rate improves as the N₂-dilution level decreases. The correlation is reasonably good for undiluted flames with equivalence ratios of 0.8 < Φ < 1.2. This result is particularly encouraging, given the relevance of undiluted flames to practical applications, and it motivates further investigation of the parameter space for which difficulties may exist in using CH as a flame marker.     

Experimental and kinetic study of 2,4,4-trimethyl-1-pentene and iso-octane in laminar flames

Laminar flame speeds of 2,4,4-trimethyl-1-pentene are investigated at equivalence ratios of 0.7–1.6, initial temperatures of 298–453 K and initial pressures of 0.1–0.5 MPa. The comparison between 2,4,4-trimethyl-1-pentene and iso-octane is also performed. Results show that 2,4,4-trimethyl-1-pentene has faster laminar flame speed than iso-octane. Chemical kinetic models (Metcalfe model, Modified model I) were tested against the present experimental data. The laminar flame speeds are apparently over-estimated by the Metcalfe model and under-predicted by the Modified model I. Therefore, high-level quantum mechanical calculations were used to revise the Modified model I to obtain Modified model II and it can give fairly good prediction at various conditions on laminar flame speeds. In addition, the chemical kinetic analysis was conducted. The analysis indicates both thermal and kinetic effects result in the discrepancy of laminar flame speeds between 2,4,4-trimethyl-1-pentene and iso-octane. Furthermore, IC₄H₈ plays a dominant role in laminar flame speeds of 2,4,4-trimethyl-1-pentene and iso-octane.