Comparative analysis of methods for heat losses in turbulent premixed flames using physically-derived reduced-order manifolds

Heat losses have the potential to substantially modify turbulent combustion processes, especially the formation of pollutants such as nitrogen oxides. The chemistry governing these species is strongly temperature sensitive, making heat losses critical for an accurate prediction. To account for the effects of heat loss in large eddy simulation (LES) using a precomputed reduced-order manifold approach, thermochemical states must be precomputed not only for adiabatic conditions but also over a range of reduced enthalpy states. However, there are a number of methods for producing reduced enthalpy states, which invoke different implicit assumptions. In this work, a set of a priori and a posteriori LES studies have been performed for turbulent premixed flames considering heat losses within a precomputed reduced-order manifold approach to determine the sensitivity to the method by which reduced enthalpy states are generated. Two general approaches are explored for generating these reduced enthalpy states and are compared in detail to assess any effects on turbulent flame structure and emissions. In the first approach, the enthalpy is reduced at the boundary of the one-dimensional (1D) premixed flame solution, resulting in a single enthalpy deficit for a single premixed flame solution. In the second approach, a variable heat loss source term is introduced into the 1D flame solutions by mimicking a real heat loss to reduce the post-flame enthalpy. The two approaches are compared in methane–air piloted turbulent premixed planar jet flames with different diluents that maintain a constant adiabatic flame temperature but experience different radiation heat losses. Both a priori and a posteriori results, as well as a chemical pathway analysis, indicate that the manner by which the heat loss is accounted for in the manifold is of secondary importance compared to other model uncertainties such as the chemical mechanism, except in situations where heat loss is unphysically fast compared to the flame time scale. A new theoretical framework to explain this insensitivity is also proposed, and its validity is briefly assessed.     

On the two-reactant one-step activation-energy asymptotics for steady,adiabatic, planar flames with Lewis numbers of unity

Aspects of predictions of activation-energy asymptotics concerning the dependence of the burning velocity on the equivalence ratio are examined here through both asymptotic analyses and numerical computation. In typical hydrocarbon–air flames, the burning velocity achieves its maximum value for fuel-rich mixture, the cause being generally attributed to the effects of detailed chemical kinetics and unequal diffusivities of the reactants. The present results demonstrate the possibility of this attribute of the burning velocity occurring even when these two effects are absent. This is accomplished by parametrically studying the burning-velocity formula valid for all equivalence ratios under the conditions specified in the title of this article, with special attention paid to implications for hydrocarbon–air flames.     

The response of turbulent premixed flames to normal acoustic excitation

To model the thermo-acoustic excitation of flames in practical combustion systems, it is necessary to know how a turbulent flame front responds to an incident acoustic wave. This will depend partly on the way in which the burning velocity responds to the wave. In this investigation, the response of CH₄/air and CH₄/H₂/air mixtures has been observed in a novel flame stabilisation configuration, in which the premixture of fuel and air is made to decelerate under controlled conditions in a wide-angle diffuser. Control is provided by an annular wall-jet of air and by turbulence generators at the inlet. Ignition from the outlet of the diffuser allows an approximately flat flame to propagate downwards and stabilise at a height that depends on the turbulent burning velocity. When the flow is excited acoustically, the ensemble-averaged height oscillates. The fluctuations in flow velocity and flame height are monitored by phase-locked particle image velocimetry and OH-planar laser induced fluorescence, respectively. The flame stabilised against a lower incident velocity as the acoustic amplitude increased. In addition, at the lowest frequency of 52 Hz, the fluctuations in turbulent burning velocity (as represented by the displacement speed) were out-of-phase with the acoustic velocity. Thus, the rate of displacement of the flame front relative to the flow slowed as the flow accelerated, and so the flame movement was bigger than it would have been if the burning velocity had not responded to the acoustic fluctuation. With an increase in frequency to 119 Hz, the relative flame movement became even larger, although the phase-difference was reduced, so the effect on burning velocity was less dramatic. The addition of hydrogen to the methane, so as to maintain the laminar burning velocity at a lower equivalence ratio, suppressed the response at low amplitude, but at a higher amplitude, the effect was reversed.     

Scalar length scales and spatial averaging effects in turbulent piloted methane/air jet flames

The effects of spatial averaging in measurements of scalar variance and scalar dissipation in three piloted methane/air jet flames (Sandia flames C, D, and E) are investigated. Line imaging of Raman scattering, Rayleigh scattering, and laser-induced CO fluorescence is applied to obtain simultaneous single-shot measurements of temperature, the mass fractions of all major species, and mixture fraction, ξ, along 7-mm segments. Spatial filters are applied to ensembles of instantaneous profiles to quantify effects of spatial averaging on the Favre mean and variance of mixture fraction and scalar dissipation at several locations in the three flames. The radial contribution to scalar dissipation, χ_r = 2D_ξ (∂ξ/∂r)², is calculated from the filtered instantaneous profiles. The variance of mixture fraction tends to decrease linearly with increasing filter width, while the mean and variance of scalar dissipation are observed to follow an exponential dependence. In each case, the observed functional dependence is used to extrapolate to zero filter width, yielding estimates of the “fully resolved” profiles of measured quantities. Length scales for resolution of scalar variance and scalar dissipation are also extracted from the spatial filtering analysis and compared with length scales obtained from spatial autocorrelations. These results provide new insights on the small scale structure of turbulent jet flames and on the spatial resolution requirements for measurements of scalar variance and scalar dissipation.     

On Flame-Edge Propagation

Contemporary interest exists for understanding how reaction zones stabilize and counter-propagate against incoming reactants. Images of flame position, morphology and dynamics are presented primarily from CH planar laser-induced fluorescence (CH-PLIF) measurements. Observations of the leading-edge flame behavior with respect to upstream propagation and recession downstream are made with sequential CH-PLIF imaging, and data have been revisited in light of the recent research of McCraw et al. (Flow Turbul Combust 70(1):83–97, 2007). It is found that in cases where a distinct branch of the outer (fuel-lean) edge of the reaction zone is present, the edge of the flame is either witnessed to propagate upstream or locally disappear. In cases where no distinct branch other than the main branch is observed, the flame is witnessed to either remain stationary or drop back downstream. These observations support the notion that structures in the low speed, outer edge of the reaction zone are involved in the upstream phase of the flame propagation.     

Effects of equivalence ratio variation on lean,stratified methane–air laminar counterflow flames

The effects of equivalence ratio variations on flame structure and propagation have been studied computationally. Equivalence ratio stratification is a key technology for advanced low emission combustors. Laminar counterflow simulations of lean methane–air combustion have been presented which show the effect of strain variations on flames stabilized in an equivalence ratio gradient, and the response of flames propagating into a mixture with a time-varying equivalence ratio. ‘Back supported’ lean flames, whose products are closer to stoichiometry than their reactants, display increased propagation velocities and reduced thickness compared with flames where the reactants are richer than the products. The radical concentrations in the vicinity of the flame are modified by the effect of an equivalence ratio gradient on the temperature profile and thermal dissociation. Analysis of steady flames stabilized in an equivalence ratio gradient demonstrates that the radical flux through the flame, and the modified radical concentrations in the reaction zone, contribute to the modified propagation speed and thickness of stratified flames. The modified concentrations of radical species in stratified flames mean that, in general, the reaction rate is not accurately parametrized by progress variable and equivalence ratio alone. A definition of stratified flame propagation based upon the displacement speed of a mixture fraction dependent progress variable was seen to be suitable for stratified combustion. The response times of the reaction, diffusion, and cross-dissipation components which contribute to this displacement speed have been used to explain flame response to stratification and unsteady fluid dynamic strain.     

Ignition and extinction of low molecular weight esters in nonpremixed flows

An experimental and kinetic modeling study is carried out to characterize combustion of low molecular weight esters in nonpremixed, nonuniform flows. An improved understanding of the combustion characteristics of low molecular weight esters will provide insights on combustion of high molecular weight esters and biodiesel. The fuels tested are methyl butanoate, methyl crotonate, ethyl propionate, biodiesel, and diesel. Two types of configuration – the condensed fuel configuration and the prevaporized fuel configuration – are employed. The condensed fuel configuration is particularly useful for studies on those liquid fuels that have high boiling points, for example biodiesel and diesel, where prevaporization, without thermal breakdown of the fuel, is difficult to achieve. In the condensed fuel configuration, an oxidizer, made up of a mixture of oxygen and nitrogen, flows over the vaporizing surface of a pool of liquid fuel. A stagnation-point boundary layer flow is established over the surface of the liquid pool. The flame is stabilized in the boundary layer. In the prevaporized fuel configuration, the flame is established in the mixing layer formed between two streams. One stream is a mixture of oxygen and nitrogen and the other is a mixture of prevaporized fuel and nitrogen. Critical conditions of extinction and ignition are measured. The results show that the critical conditions of extinction of diesel and biodiesel are nearly the same. Experimental data show that in general flames burning the esters are more difficult to extinguish in comparison to those for biodiesel. At the same value of a characteristic flow time, the ignition temperature for biodiesel is lower than that for diesel. The ignition temperatures for biodiesel are lower than those for the methyl esters tested here. Critical conditions of extinction and ignition for methyl butanoate were calculated using a detailed chemical kinetic mechanism. The results agreed well with the experimental data. The asymptotic structure of a methyl butanoate flame is found to be similar to that for many hydrocarbon flames. This will facilitate analytical modeling, of structures of ester flames, using rate-ratio asymptotic techniques, developed previously for hydrocarbon flames.     

Dilution effects analysis of opposed-jet H2/CO syngas diffusion flames

This paper reported the analysis of dilution effects on the opposed-jet H₂/CO syngas diffusion flames. A computational model, OPPDIF coupled with narrowband radiation calculation, was used to study one-dimensional counterflow syngas diffusion flames with fuel side dilution from CO₂, H₂O and N₂. To distinguish the contributing effects from inert, thermal/diffusion, chemical, and radiation effects, five artificial and chemically inert species XH₂, XCO, XCO₂, XH₂O and XN₂ with the same physical properties as their counterparts were assumed. By comparing the realistic and hypothetical flames, the individual dilution effects on the syngas flames were revealed. Results show, for equal-molar syngas (H₂/CO = 1) at strain rate of 10 s^?1, the maximum flame temperature decreases the most by CO₂ dilution, followed by H₂O and N₂. The inert effect, which reduces the chemical reaction rates by behaving as the inert part of mixtures, drops flame temperature the most. The thermal/diffusion effect of N₂ and the chemical effect of H₂O actually contribute the increase of flame temperature. However, the chemical effect of CO₂ and the radiation effect always decreases flame temperature. For flame extinction by adding diluents, CO₂ dilution favours flame extinction from all contributing effects, while thermal/diffusion effects of H₂O and N₂ extend the flammability. Therefore, extinction dilution percentage is the least for CO₂. The dilution effects on chemical kinetics are also examined. Due to the inert effect, the reaction rate of R84 (OH+H₂ = H+H₂O) is decreasing greatly with increasing dilution percentage while R99 (CO+OH→CO₂+H) is less affected. When the diluents participate chemically, reaction R99 is promoted and R84 is inhibited with H₂O addition, but the trend reverses with CO₂ dilution. Besides, the main chain-branching reaction of R38 (H+O₂→O+OH) is enhanced by the chemical effect of H₂O dilution, but suppressed by CO₂ dilution. Relatively, the influences of thermal/diffusion and radiation effects on the reaction kinetics are then small.     

On nonlinear model equations for the response of premixed flames to acoustic like accelerations

The response of a dynamical flame model to imposed acoustic accelerations is studied analytically and numerically. Through linear stability analyses, two analytical approximations for the primary and the parametric stability boundaries are found. The approximation for the primary instability boundary is accurate for any periodic accelerations, in the limit of large acoustic frequencies. The critical acoustic amplitude u _a for Landau–Darrieus instability suppression is identified and found to depend only on the density contrast and the shape of the periodic acoustic stimuli. The proposed model evolution equation is next integrated numerically with various imposed acoustic accelerations; the primary and parametric flame responses are identified. It is shown analytically and numerically that in the presence of a fully developed, yet weakened by acoustics, Landau–Darrieus (or primary) instability the wrinkle amplitude and the mean flame speed oscillate at the same frequency as the acoustic stimuli; the threshold for suppression of primary instability by acoustic forcing is determined exactly. Increasing the acoustic amplitude allows the flame to respond parametrically to the acoustics. This response is characterised by troughs and crests interchanging their roles while the mean flame speed again oscillates with the same frequency as the acoustic stimuli and at twice that of wrinkle amplitude oscillations.