Theoretical analysis of the evolution from ignition kernel to flame ball and planar flame

Dynamics of flame kernel evolution with and without external energy addition has been investigated analytically and numerically. Emphasis is placed on the effects of radiation heat loss, ignition power and Lewis number on the correlation and transition between the initial flame kernel, the self-extinguishing flame, the flame ball, the outwardly propagating spherical flame and the propagating planar flame. The present study extends previous results by bridging the theories of the non-adiabatic stationary flame balls and travelling flames and allowing rigorous consideration of radiation heat losses. The results show that the effects of radiation heat loss play an important role in flame regimes and flame transition and result in a new isolated self-extinguishing flame. Furthermore, it is found that radiation heat losses significantly increase the critical ignition radius and result in three different dependences of the minimum ignition power on the Lewis number. Comparisons between the results from the transient numerical simulation and those from the quasi-steady state analysis show a good agreement. The results suggest that prediction of flame initiation without appropriate consideration of radiation is not acceptable.     

A numerical study on the effect of CO addition on extinction limits and NO x formation in lean counterflow CH4/air premixed flames

The effect of CO addition on extinction and NO _x formation in lean premixed counterflow CH₄/air flames was investigated by numerical simulation. Detailed chemistry and complex thermal and transport properties were employed. A method that gradually switched off the initial reactions of NO formation from different routes was used to analyse the variation of NO formation mechanism. The results indicate that the addition of certain amount of CO increases the strain extinction limits and reduces the radiation extinction limits. As a result, the lean flammability limit of CH₄/air premixed flame is extended to leaner side by the addition of CO. The formation of NO in a flame is increased with the addition of CO at a constant equivalence ratio. For an ultra-lean flame, the increase in the formation of NO is mainly because of the increase in the contribution from the NNH intermediate route, while for a near stoichiometric flame, this increase is mainly attributable to the rise in the contribution from the thermal route. With the fraction of added CO being gradually increased, the formation of NO₂ in a flame first decreases and then increases at a given equivalence ratio. The addition of CO reduces the formation N₂O in an ultra-lean flame, while affects little on the formation of N₂O in a near stoichiometric flame.     

A DNS study on the flame structures and flame stabilization mechanism of a laboratory-scale lean premixed jet flame in crossflow

In the present work, direct numerical simulation (DNS) of a laboratory-scale lean premixed reacting jet flame in crossflow was performed to understand the flame structures and the flame stabilization mechanism. In the DNS, an ethylene-air jet with an equivalence ratio of 0.6 was injected into a hot vitiated crossflow. The jet Reynolds number reaches 6161. The DNS results were compared with those of the experiment with a good agreement. It was found that the windward and leeward branches of the flame show significantly different behaviors. The windward flame branch, appearing lifted and discontinuous, is located in the shear layer regions with high temperature, low vorticity and low scalar dissipation rate. The location of the peak heat release rate shifts to a higher mixture fraction with increasing distance from the jet exit. The leeward branch of the flame anchors in the shear layer near the jet exit. The recirculation zone in the wake of the jet facilitates the stabilization of the leeward flame. The chemical explosive mode analysis (CEMA) and species budget analysis were employed to characterize the local combustion mode. Auto-ignition plays a key role in the stabilization of the windward flame where a large range of extinction is also found due to the high strain rate. In contrast, premixed flame propagation is dominant on the leeward side.     

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.     

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.     

Computational study of influences of a seam line of a ball for baseball on flows

Flows around a ball used in baseball games are calculated using third-order upwinddifference method with various seam positions determined by two rotation angles. Those are four-seam rotation with an angle: a and two-seam rotation with an angle: b. The computed results of the four-seam rotation are compared with experimental data measured in a wind tunnel and computed drag coefficients qualitatively agree well with experiments. However, lift coefficients do not agree well. The computed results and geometrical symmetry suggest that a supporting rod in the wind tunnel would have strong influence on the accuracy of the measurement. Flow changes in two-seam rotation are also simulated. It is found that the lowest drag force is observed atb=90 and that the value is less than half of the largest drag force ata=30 and 60 degrees. The largest lift force is observed atb=20 degree. In this case, a projection of the seam line on the top causes a large separation while smooth surface without the seam at the bottom dose not separate the flow. A pair of longitudinal vortices are found in the wake, which make wake slant and generate large lift force.     

Experimental and kinetic modeling study of 1-hexanol combustion in an opposed-flow diffusion flame

Biofuels are of particular interest as they have the potential to reduce our dependence on petroleum-derived fuels and the levels of greenhouse gas emissions from transportation. 1-Hexanol is a promising renewable long chain alcohol that can be used in conventional fuel blends or as a cosolvent for biodiesel mixtures. However, the fundamental combustion properties of 1-hexanol have not been fully characterized in the literature. Thus, new experimental results, consisting of temperature and concentration profiles of stable species were obtained for the oxidation of 1-hexanol generated in an opposed-flow diffusion flame at 0.101 MPa. This experimental data were compared to the predicted values of a detailed chemical kinetic mechanism proposed in the literature to study the combustion of 1-hexanol.This mechanism consists of 361 chemical species and 2687 chemical reactions (most of them reversible). Reaction pathway and sensitivity analyses were performed to interpret the results. In addition, improvements were investigated to optimize the proposed mechanism.     

Microgravity droplet combustion: effect of tethering fiber on burning rate and flame structure

Droplets tethering on fibers has become a well established technique for conducting droplet combustion experiments in microgravity conditions. The effects of these supporting fibers are frequently assumed to be negligible and are not considered in the experimental analysis or in numerical simulations. In this work, the effect of supporting fibers on the characteristics of microgravity droplet combustion has been investigated numerically; a priori predictions have then been compared with published experimental data. The simulations were conducted using a transient one-dimensional spherosymmetric droplet combustion model, where the effect of the supporting fiber was implicitly taken into account. The model applied staggered convective flux finite volume method combined with high-order implicit time integration. Thermal radiation was evaluated using a statistical narrow band radiation model. Chemical kinetics and thermophysical properties were represented in rigorous detail. Tether fiber diameter, droplet diameter, ambient pressure and oxygen concentration were varied over a range for n-decane droplets in the simulations. The results of the simulations were compared to previously published experiments conducted in the Japan Microgravity Center (JAMIC) 10 second drop tower and the NASA Glenn Research Center (GRC) 5.2 second drop tower. The model reproduces closely nearly all aspects of tethered n-decane droplet burning phenomena, which included droplet burning history, transient and average burning rate, and flame standoff ratio. The predictions show that the presence of the tethering fiber significantly influences the observed burning rate, standoff ratio, and extinction.     

Spatiotemporal measurements of flame stretch and propagation rates for lean and rich CH4/air premixed flames interacting with a turbulent-wake

A 1.5 m long turbulent-wake combustion vessel with a 0.15 m × 0.15 m cross-sectional area is proposed for spatiotemporal measurements of curvature, strain, dilatation and burning rates along a freely downward-propagating premixed flame interacting with a parallel row of staggered vortex pairs having both compression (negative) and extension (positive) strains simultaneously. The wanted wake is generated by rapidly withdrawing an electrically-controlled, horizontally-oriented sliding plate of 5 mm thickness for flame–wake interactions. Both rich and lean CH₄/air flames at the equivalence ratios  = 1.4 and  = 0.7 with nearly the same laminar burning velocity are studied, where flame–wake interactions and their time-dependent velocity fields are obtained by high-speed, high-resolution DPIV and laser-tomography. Correlations among curvature, strain, stretch, and dilatation rates along wrinkled flame fronts at different times are measured and thus their influences on front propagation rates can be analyzed. It is found that strain-related effects have significant influence on front propagation rates of rich CH₄/air (diffusionally stable) flames even when the curvature weights more in the total stretch than the strain rate does. The local propagation rates along the wrinkled flame front are more intense at negative strain rates corresponding to positive peak dilatation rates but the global propagation rate averaged along the rich flame front remains constant during all period of flame–wake interaction. For lean CH₄/air (diffusionally unstable) flames, the curvature becomes a dominant parameter influencing the structure and propagation of the wrinkled flame front, where both local and global propagation rates increase significantly with time, showing unsteady flame propagation. These experimental results suggest that the theory of laminar flame stretch can be applicable to a more complex flame–wake interaction involving unsteadiness and multitudinous interactions between vortices.     

Control of oscillating combustion and noise based on local flame structure

To control combustion oscillations, the characteristics of an oscillating swirl injection premixed flame have been investigated, and control of oscillating combustion and noise based on local flame structure has been conducted. The r.m.s. value of pressure fluctuations and noise level show significantly large values between = 0.8 and 1.1. The beating of pressure fluctuations is observed for the large oscillating flame conditions in this combustor. Relationship between beating of pressure fluctuations and local flame structure was observed by the simultaneous measurement of CH/OH planar laser induced fluorescence and pressure fluctuations. The local flame structure and beating of pressure fluctuations are related and the most complicated flame is formed in the middle pressure fluctuating region of beating. The beating of pressure fluctuations, which plays important roles in noise generation and nitric oxide emission in this combustor, could be controlled by injecting secondary fuel into the recirculating region of oscillating flames. Injecting secondary fuel prevented lean blowout, and low NO_x combustion was also achieved even for the case of pure methane injection as a secondary fuel. By injecting secondary fuel into the recirculating region near the swirl injector, the flame lifted from the swirl injector and its reaction region became uniform and widespread, hence resulting in low nitric oxide emission. Secondary mixture injection, fuel diluted with air, is not effective for control of combustion oscillations suppression and lean blowout prevention.     

Effects of turbulence-chemistry interactions on auto-ignition and flame structure for n-dodecane spray combustion

The Engine Combustion Network (ECN) spray A under diesel engine conditions is investigated with a non-adiabatic 5D Flamelet Generated Manifolds (FGM) model with the consideration of detailed chemical kinetic mechanisms. The enthalpy deficit due to droplet vapourisation is considered by employing an additional controlling parameter in the FGM library. In this FGM model, β-PDF is used for the PDF integration over the control variable space. Validation results in non-reacting conditions indicate relatively good agreement between the predicted and experimental data in terms of liquid and vapour penetrations and mixture fraction spatial distribution. In reacting conditions, the effects of variance of mixture fraction and progress variable were examined. The ignition delay time and the quasi-steady flame structure are both affected by the variances. The variance of mixture fraction delays the ignition process and the variance of progress variable accelerates it. For mixture fraction, the ignition process is quicker at any stage in the case of neglecting variance. While things are more complex for progress variable, the ignition process is advanced in the case of neglecting variance at early times, but surpassed by the case of β-PDF later and until auto-ignition. When variance of mixture fraction is considered, the OH mass fraction shows a wide spatial distribution. While if not, a very thin flame is observed with a higher peak in OH, and a very large lift-off length. The variance of progress variable has little impact on the global flame structure, but makes the flame lift-off length much shorter. This study confirms the general observation, that the variance of mixture fraction is of higher importance in high temperature non-premixed combustion, however, we found that the variance of progress variable is far from negligible.     

Experimental and numerical study of premixed, lean ethylene flames

Ethylene is a key intermediate in the combustion mechanisms of most practical fuels. It plays also an important role in the formation of aromatic hydrocarbons and soot particules. The latter has motivated many experimental and numerical studies carried out on rich ethylene-air mixtures. Less studies have been devoted to lean mixtures, and the development of strategies based on lean, premixed flames to reduce soot and NO_x production requires additional experimental data in lean conditions. In this work, the chemical structure of lean premixed ethylene-oxygen-nitrogen flames stabilized on a flat-flame burner at atmospheric pressure was determined experimentally. The species mole fraction profiles were also computed by the Premix code (Chemkin II version) and four detailed reaction mechanisms. A very good agreement was observed for the main flame properties: reactants consumption, final products (CO₂, H₂O) and the main intermediates: CO and H₂. Marked differences occurred in the prediction of active intermediate species present in small concentrations. Pathways analyses were performed to identify the origins of these discrepancies. It was shown that the same reactions were involved in the four mechanisms to describe the consumption of ethylene, but with marked differences in their relative importance. C₂H₃ and CH₂HCO are the main radicals formed in this first step and their consumption increases the differences between the mechanisms either by the use of different kinetic data for common reactions or by differences in the nature of the consumption reactions.     

Effects of subgrid scale and combustion modelling on flame structure of a turbulent premixed flame within LES and tabulated chemistry framework

The effects of combustion and SubGrid Scale (SGS) modelling on the overall flame characteristics of a turbulent premixed flame are investigated. This is achieved in terms of mean flow statistics, variances and flame surfaces. In particular, the chemical flame structure is analysed and compared. The Artificially Thickened Flame (ATF) approach coupled with the Flamelet Generated Manifolds (FGMs) and Filtered TAbulated Chemistry for LES (F-TACLES) approaches are used for this investigation. A Germano like procedure for dynamical calculation of SGS wrinkling is used which ensures the conservation of the total flame surface for both models. It turns out that using the dynamic SGS wrinkling model improves the results. Although the results of both combustion models in terms of statistics, mean and variances show very good agreement, the resolved flame surfaces hide different dynamic behaviour.     

Modeling of flame lift-off length in diesel low-temperature combustion with multi-dimensional CFD based on the flame surface density and extinction concept

Low-Temperature Combustion (LTC) is becoming a promising technology for simultaneously reducing soot and NOx emissions from diesel engines. LTC regimes are evaluated by the flame lift-off length – the distance from the injector orifice to the location of hydroxyl luminescence closest to the injector in the flame jet. Various works have been dedicated to successful simulations of lifted flames of a diesel jet by use of various combustion modeling approaches. In this work, flame surface density and flamelet concepts were used to model the diesel lift-off length under LTC conditions. Numerical studies have been performed with the ECFM3Z model, n-Heptane and diesel fuels to determine the flame lift-off length and its correlation with soot formation under quiescent conditions. The numerical results showed good agreement with experimental data, which were obtained from an optically accessible constant volume chamber and presented at the Engine Combustion Network (ECN) of Sandia National Laboratories. It was shown that at a certain distance downstream from the injector orifice, stoichiometric scalar dissipation rate matched the extinction scalar dissipation rate. This computed extinction scalar dissipation rate correlated well with the flame lift-off length. For the range of conditions investigated, adequate quantitative agreement was obtained with the experimental measurements of lift-off length under various ambient gas O₂ concentrations, ambient gas temperatures, ambient gas densities and fuel injection pressures. The results showed that the computed lift-off length values for most of the conditions lay in a reasonable range within the quasi-steady lift-off length values obtained from experiments. However, at ambient temperatures lower than 1000 K, the lift-off length values were under-predicted by the numerical analysis. This may be due to the use of the droplet evaporation model as it is believed that evaporation has a strong effect on the lift-off length.     

Comparison of laser-induced and planar laser-induced fluorescence measurements of nitric oxide in a high-pressure, swirl-stabilized, spray flame

We report spatially resolved linear laser-induced fluorescence (LIF) and planar laser-induced fluorescence (PLIF) measurements of nitric oxide (NO) in a pre-heated, high-pressure (4.27 atm), lean direct-injection (LDI) spray flame. The feasibility of using PLIF in lieu of LIF is assessed with respect to measuring NO concentrations in high-pressure LDI spray flames. NO is excited via the resonant Q₂(26.5) transition of the γ(0,0) band while a non-resonant wavelength is employed to subtract background interferences. LIF detection is performed in a 2-nm region centered on the γ(0,1) band. PLIF detection is performed in a 68-nm window that captures fluorescence from several vibrational bands. An in situ NO doping scheme for fluorescence calibration is successfully employed to quantify the LIF signals. However, a similar calibration scheme for the reduction of PLIF images to quantitative field measurements is plagued by the laser-excited background. Excitation scans and calibration comparisons have been performed to assess the background contribution for PLIF detection. Quantitative radial NO profiles measured by LIF are presented and analyzed so as to correct the PLIF measurements to within the accuracy bars of the LIF measurements via a single-point scaling of the PLIF image. Received: 23 November 1999 / Revised version: 17 January 2000 / Published online: 27 April 2000     

Experimental and numerical investigations on the interactions of volatile flame and char combustion of a coal particle

The model that takes chemical reactions, heat and mass transfers in the boundary layer of the particle into account simultaneously, is developed for simulating the combustion of a pulverized coal particle. The FTIR in situ temperature-measurements and the comparison between numerical simulations for the pulverized coal and the devolatilized char show that the volatile flame induces the combustion of the primary product of surface oxidation CO. Due to the influence of volatile flame, the char particle can be ignited at temperature lower than its heterogeneous ignition temperature, which elucidates the physical essence of joint hetero-homogeneous ignition mode discovered by Jüntgen.     

Characterizing the fuel-specific combustion chemistry of acetic acid and propanoic acid: Laminar flame propagation and kinetic modeling studies

In order to study the combustion chemistry of carboxyl functionality, the laminar burning velocity of acetic acid/air and propanoic acid/air mixtures was investigated in a high-pressure constant-volume cylindrical combustion vessel at 423 K, 1 atm and equivalence ratios of 0.7–1.4. Experimental results reveal that the flame propagation of propanoic acid flame is much faster than that of acetic acid flame, especially under rich conditions, and the laminar burning velocity of propanoic acid/air mixtures peaks at richer conditions than that of acetic acid. The present theoretical calculations for the isomerization and decomposition of propanoic acid radicals indicate that the primary radical products are HOCO, H and C₂H₅, while those in acetic acid flame are CH₃ and OH based on previous studies. A kinetic model of the two acids was developed mainly based on previous and the present theoretical calculation results. It could reasonably capture the measured laminar burning velocities of acetic acid/air and propanoic acid/air mixtures in this work, as well as the previous experimental data in literature. Based on the present model, CH₃- and ketene-related pathways play an important role in acetic acid flames. Under rich conditions, ketene is mostly converted to CH₃ via CH₂CO+HCH₃+CO, and the chain-termination reaction of CH₃+H(+M)=CH₄(+M) is enhanced, which strongly inhibits the propagation of rich acetic acid flames. In contrast, C₂H₅ and ethylene chemistry play an important role in propanoic acid flames. Rich conditions promote the decomposition of C₂H₅, yielding ethylene and H, which can facilitate the flame propagation. This can explain the shift of the peak laminar burning velocity of propanoic acid/air mixtures towards a slightly richer condition compared with that of acetic acid/air mixtures.