Conditional statistics of inert droplet effects on turbulent combustion in reacting mixing layers

Direct numerical simulation (DNS) of turbulent reacting mixing layers laden with evaporating inert droplets is used to assess the droplet effects in the context of the conditional moment closure (CMC) for multiphase turbulent combustion. The temporally developing mixing layer has an initial Reynolds number of 1000 based on the vorticity thickness with more than 16 million Lagrangian droplets traced. An equivalent mixture fraction incorporating the inert vapour mass fractions clearly demonstrates the effects of vapour dilution on the mixture. Instantaneous fields and conditional statistics, such as the singly conditioned scalar dissipation rate, the gas temperature 〈 T _g|η 〉, conditional variance of the gas temperature 〈 T _g ^”2|η 〉 and conditional covariance between the fuel mass fraction and gas temperature 〈 Y _f ^” T _g ^”|η 〉 show considerable droplet effects. Comparison between the droplet-free and droplet-laden reacting mixing layer cases suggests significant extinction in the latter case. The resulting large conditional fluctuations around the conditional means contradict the basic assumption behind the first-order singly conditioned CMC. More sophisticated CMC approaches, such as doubly conditioned or second-order CMCs are, in principle, better able to incorporate the effects of evaporating droplets, but significant modelling challenges exist. The scalar dissipation rate doubly conditioned on the mixture fraction and a normalized gas temperature 〈 χ | η, ζ 〉 exemplifies the modelling complexity in the CMC of multiphase combustion.     

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.     

A LES-CMC formulation for premixed flames including differential diffusion

A finite volume large eddy simulation–conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane–air flame with Le_eff = 0.99 and a lean hydrogen–air flame with Le_eff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane–air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.     

The effect of timescale variation in multiple mapping conditioning mixing of PDF calculations for Sandia Flame series D–F

A stochastic implementation of the multiple mapping conditioning (MMC) model has been used for the modelling of turbulence–chemistry interactions in a series of turbulent jet diffusion flames with varying degrees of local extinction (Sandia Flames D–F). The mapping function approximates the cumulative probability distribution of mixture fraction and the corresponding variance can be controlled by a standard implementation of the scalar mixing timescale. The conditional fluctuations are controlled by a minor dissipation timescale, τ_min. The results show a clear dependence of the conditional fluctuations on the choice of the minor timescale, and the appropriate value for turbulent jet flames is similar to values determined in related direct numerical simulation (DNS) studies of homogeneous turbulent reacting flows. The predictions of means and variances of temperature and species mass fractions are very good for all flames, indicating an appropriate modelling of the conditional variances. Further sensitivity studies with respect to particle number density demonstrate a relative insensitivity of the results to the particle number in the numerical solution procedure. Good results can be obtained with as few as 10 particles per cell, allowing for a computationally inexpensive implementation of a Monte Carlo/probability density function (PDF) method.     

Measurements of flame orientation and scalar dissipation in turbulent partially premixed methane flames

Three turbulent flames were studied using a new experimental facility developed at Sandia National Laboratories. Line imaging of Raman and Rayleigh scattering and CO laser-induced fluorescence (LIF) yielded information on all major species, temperature, mixture fraction, and a 1D surrogate measure of scalar dissipation. Simultaneously, crossed planar OH LIF imaging provided information on the instantaneous flame orientation, allowing estimation of the full 3D (flame-normal) scalar dissipation rate. The three flames studied were methane–air piloted jet flames (Sandia flames C, D, and E), which cover a range in Reynolds number from 13,400 to 33,600. The statistics of the instantaneous flame orientation are examined in the different flames, with the purpose of studying the prevailing kinematics of isoscalar contours. The 1D and 3D results for scalar dissipation rate are examined in detail, both in the form of conditional averages and in the form of probability density functions. The effect of overall strain and Reynolds number on flame suppression and eventual extinction is also investigated, by examining the doubly conditional statistics of temperature in the form of S-shaped curves. This latter analysis reveals that double conditioning of temperature on both mixture fraction and scalar dissipation does not collapse the data from these flames onto the same curve at low scalar dissipation rates, as might be expected from simple flamelet concepts.     

Simulations of turbulent lifted jet flames with two-dimensional conditional moment closure

Simulations of H₂ air lifted jet flames are presented, obtained in terms of two-dimensional, first-order conditional moment closure (CMC). The unsteady CMC equation with detailed chemistry is solved without the need for operator splitting, while the accompanying flow field is determined using commercial CFD software employing a k − ε turbulence model. Computed lift-off heights and Favre-averaged species mole fractions are found to be very close to values obtained experimentally for a wide range of jet velocities and fuel–air mixtures. Simulations for which the initial condition is an attached flame and the jet velocity gradually increased do not result in lift-off, a result fully consistent with experimental observation and capturing the hysteresis behaviour of lifted flames. The stabilisation mechanism is explored by quantifying the balance of terms comprising the CMC in the lift-off region. In line with experimental data, it is found that the scalar dissipation rate at the stabilisation height is well below the extinction value, and that axial transport and molecular diffusion play a major role. The radial components of spatial convection and diffusion are always small, fully justifying the alternative approach of employing a cross-stream averaged CMC.     

Second-order conditional moment closure modeling of a turbulent CH4/H2/N2 jet diffusion flame

The second-order CMC model for a detailed chemical mechanism is used to model a turbulent CH₄/H₂/N₂ jet diffusion flame. Second-order corrections are made to the three rate limiting steps of methane–air combustion, while first-order closure is employed for all the other steps. Elementary reaction steps have a wide range of timescales with only a few of them slow enough to interact with turbulent mixing. Those steps with relatively large timescales require higher-order correction to represent the effect of fluctuating scalar dissipation rates. Results show improved prediction of conditional mean temperature and mass fractions of OH and NO. Major species are not much influenced by second-order corrections except near the nozzle exit. A parametric study is performed to evaluate the effects of the variance parameter in log-normal scalar dissipation PDF and the constants for the dissipation term in conditional variance and covariance equations.     

Analysis of auto-ignition of heated hydrogen–air mixtures with different detailed reaction mechanisms

Auto-ignition processes of hydrogen, diluted with nitrogen, in heated air are numerically investigated by means of an unsteady laminar flamelet approach in mixture fraction space. The focus is on the auto-ignition delay time and the most reactive mixture fraction as obtained with five chemical mechanisms. Two strongly different levels of dilution, corresponding to experiments in the open literature, are considered. This concerns low-temperature chemistry at atmospheric pressure. The temperature of the air stream is much higher than the temperature of the fuel stream in the cases under study. We extensively investigate the effect of the co-flow temperature, the conditional scalar dissipation rate and the resolution in mixture fraction space for one case. With respect to the conditional scalar dissipation rate, we discuss the Amplitude Mapping Closure (AMC) model with imposed maximum scalar dissipation rate at mixture fraction equal to 0.5, as well as a constant conditional scalar dissipation rate value over the entire mixture fraction value range. We also illustrate that an auto-ignition criterion, based on a temperature rise, leads to similar results as an auto-ignition criterion, based on OH mass fraction, provided that the hydrogen is not too strongly diluted.     

Direct numerical simulation of a statistically stationary,turbulent reacting flow

An inhomogeneous, non-premixed, stationary, turbulent, reacting model flow that is accessible to direct numerical simulation (DNS) is described for investigating the effects of mixing on reaction and for testing mixing models. The mixture-fraction-progress-variable approach of Bilger is used, with a model, finite-rate, reversible, single-step thermochemistry, yielding non-trivial stationary solutions corresponding to stable reaction and also allowing local extinction to occur. There is a uniform mean gradient in the mixture fraction, which gives rise to stationarity as well as a flame brush. A range of reaction zone thicknesses and Damkohler numbers are examined, yielding a broad spectrum of behaviour, including thick and thin flames, local extinction and near equilibrium. Based on direct numerical simulations, results from the conditional moment closure (CMC) and the quasi-equilibrium distributed reaction (QEDR) model are evaluated. Large intermittency in the scalar dissipation leads to local extinction in the DNS. In regions of the flow where local extinction is not present, CMC and QEDR based on the local scalar dissipation give good agreement with the DNS.M This article features multimedia enhancements available from the supplemental page in the online journal.     

Direct numerical simulations of autoignition in turbulent two-phase flows

Three-dimensional direct numerical simulations (DNS) were carried out to investigate the impact of evaporation of droplets on the autoignition process under decaying turbulence. The droplets were taken as point sources and were tracked in a Lagrangian manner. Three cases with the same initial equivalence ratio but different initial droplet size were simulated and the focus was to examine the influence of the droplet evaporation process on the location of autoignition. It was found that an increase in the initial droplet size results in an increase in the autoignition time, that highest reaction rates always occur at a specific mixture fraction ξ_MR, as in purely gaseous flows, and that changes in the initial droplet size did not affect the value of ξ_MR. The conditional correlation coefficient between scalar dissipation rate and reaction rates was only mildly negative, contrary to the strongly negative values for purely gaseous autoigniting flows, possibly due to the continuous generation of mixture fraction by the droplet evaporation process that randomizes both the mixture fraction and the scalar dissipation fields.     

Direct numerical simulations and analysis of three-dimensional n-heptane spray flames in a model swirl combustor

Three-dimensional n-heptane spray flames in a swirl combustor are investigated by means of direct numerical simulation (DNS) to provide insight into realistic spray evaporation and combustion as well as relevant modeling issues. The variable-density, low-Mach number Navier–Stokes equations are solved using a fully conservative and kinetic energy conserving finite difference scheme in cylindrical coordinates. Dispersed droplets are tracked in a Lagrangian framework. Droplet evaporation is described by an equilibrium model. Gas combustion is represented using an adaptive one-step irreversible reaction. Two different cases are studied: a lean case that resembles a lean direct injection combustion, and a rich case that represents the primary combustion region of a rich-burn/quick-quench/lean-burn combustor. The results suggest that premixed combustion contribute more than 70% to the total heat release rate, although diffusion flame have volumetrically a higher contribution. The conditional mean scalar dissipation rate is shown to be strongly influenced, especially in the rich case. The conditional mean evaporation rate increases almost linearly with mixture fraction in the lean case, but shows a more complex behavior in the rich case. The probability density functions (PDF) of mixture fraction in spray combustion are shown to be quite complex. To model this behavior, the formulation of the PDF in a transformed mixture fraction space is proposed and demonstrated to predict the DNS data reasonably well.     

Modeling evaporation effects in conditional moment closure for spray autoignition

Simulations of an n-heptane spray autoigniting under conditions relevant to a diesel engine are performed using two-dimensional, first-order conditional moment closure (CMC) with full treatment of spray terms in the mixture fraction variance and CMC equations. The conditional evaporation term in the CMC equations is closed assuming interphase exchange to occur at the droplet saturation mixture fraction values only. Modeling of the unclosed terms in the mixture fraction variance equation is done accordingly. Comparison with experimental data for a range of ambient oxygen concentrations shows that the ignition delay is overpredicted. The trend of increasing ignition delay with decreasing oxygen concentration, however, is correctly captured. Good agreement is found between the computed and measured flame lift-off height for all conditions investigated. Analysis of source terms in the CMC temperature equation reveals that a convective–reactive balance sets in at the flame base, with spatial diffusion terms being important, but not as important as in lifted jet flames in cold air. Inclusion of droplet terms in the governing equations is found to affect the mixture fraction variance field in the region where evaporation is the strongest, and to slightly increase the ignition delay time due to the cooling associated with the evaporation. Both flame propagation and stabilization mechanisms, however, remain unaffected.     

Modelling of lifted turbulent diffusion flames in a channel mixing layer by the flame hole dynamics

The partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics as an effort to develop a prediction model for the turbulent flame lift off. The essence of the flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of the quenching holes initially created by the local quenching events. The numerical simulation for the flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for a constant-density fuel–air channel mixing layer to obtain the background turbulent flow and mixing fields, from which a time series of two-dimensional scalar-dissipation-rate array is extracted. Subsequently, a Lagrangian simulation of the flame hole random walk mapping, projected to the scalar dissipation rate array, yields a temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. In particular, the probability of encountering the reacting state, while conditioned with the instantaneous scalar dissipation rate, is examined to reveal that the conditional probability has a sharp transition across the crossover scalar dissipation rate, at which the flame edge changes its direction of propagation. This statistical characteristic implies that the flame edge propagation instead of the local quenching event is the main mechanism controlling the partial quenching events in turbulent flames. In addition, the conditional probability can be approximated by a heavyside function across the crossover scalar dissipation rate.     

Air/fuel mixing in jet flames

The paper examines eight diverse regimes in which fuels can mix and react with air. These comprise: (i) Lifted subsonic; and (ii) supersonic jet flames, with (iii) and without (iv) cross flows; (v) Rim-attached flames; (vi) Early Downwash flames; (vii) Downwash-attached jet flames; and (viii) Fire Whirls.Correlations of characteristics within these regimes are principally in terms of a dimensionless Flow Number, U*, Cross Flow Reynolds number, Re_c, and, for Fire Whirls, a dimensionless Critical Velocity, CV. Boundaries of seven of the eight regimes are identified, through plots of U*, against Re_c, and of the eighth through a plot of CV against U*. The circumstances of transitions between regimes are identified. The study involves a variety of CH₄ cross flow flame measurements, in a wind tunnel. Cross flows can initially create a small lee-side flame downwash, due to the depression in pressure. With increasing fuel flow this might extend 1.3 m downwards from the horizontal tip of the vertical burner. Jet flames can attach to the downwash, which can become significant above Re_c ≈ 2000. More extensive downwash might further delay blow-off. Regime boundaries are constructed on the U*/Re_c diagram covering lifted flames, early downwash, and downwash-attached flames. The most powerful flames tend to be lifted, choked, flames, with cross flow, and fire whirls. Combustion becomes less efficient at high Re_c and low U*, although CH₄ was efficiently reacted.Experimental values of the ratio of fuel to air velocity, u/u_c, of CH₄ flames ranged between about 10 and 30 for lifted flames, and between 0.3 and 3.6, at blow-off, for rim-attached flames. The latter comprise an important category, often intermediate between lifted flames and downwash-attached flames.     

Simulations of laminar non-premixed flames of kerosene with hot combustion products as oxidiser

The effects of hot combustion product dilution in a pressurised kerosene-burning system at gas turbine conditions were investigated with laminar counterflow flame simulations. Hot combustion products from a lean (φ = 0.6) premixed flame were used as an oxidiser with kerosene surrogate as fuel in a non-premixed counterflow flame at 5, 7, 9 and 11 bar. Kerosene-hot product flames, referred to as ‘MILD’, exhibit a flame structure similar to that of kerosene–air flames, referred to as ‘conventional’, at low strain rates. The Heat Release Rate (HRR) of both conventional and MILD flames reflects the pyrolysis of the primary and intermediate fuels on the rich side of the reaction zone. Positive HRR and OH regions in mixture fraction space are of similar width to conventional kerosene flames, suggesting that MILD flames are thin fronts. MILD flames do not exhibit typical extinction behaviour, but gradually transition to a mixing solution at very high rates of strain (above A = 160, 000 s^?1 for all pressures). This is in agreement with literature that suggests heavily preheated and diluted flames have a monotonic S-shaped curve. Despite these differences in comparison with kerosene–air flames, MILD flames follow typical trends as a function of both strain and pressure. Further still, the peak locations of the overlap of OH and CH₂O mass fractions in comparison with the peak HRR indicate that the pixel-by-pixel product of OH- and CH₂O-PLIF signals is a valid experimental marker for non-premixed kerosene MILD and conventional flames.     

Modeling of filtered heat release for large eddy simulation of compressible infinitely fast reacting flows

A priori and a posteriori studies for large eddy simulation of the compressible turbulent infinitely fast reacting shear layer are presented. The filtered heat release appearing in the energy equation is unclosed and the accuracy of different models for the filtered scalar dissipation rate and the conditional filtered scalar dissipation rate of the mixture fraction in closing this term is analyzed. The effect of different closures of the subgrid transport of momentum, energy and scalars on the modeling of the filtered heat release via the resolved fields is also considered. Three explicit models of these subgrid fluxes are explored, each with an increasing level of reconstruction and all of them regularized by a Smagorinsky-type term. It is observed that a major part of the error in the prediction of the conditional filtered scalar dissipation comes from the unsatisfactory modeling of the filtered dissipation itself. The error can be substantial in the turbulent fluctuation (rms) of the dissipation fields. It is encouraging that all models give good predictions of the mean and rms density in a posteriori LES of this flow with realistic heat release corresponding to large density change. Although a posteriori results show a small sensitivity to subgrid modeling errors in the current problem, extinction–reignition phenomena involving finite-rate chemistry would demand more accurate modeling of the dissipation rates. A posteriori results also show that the resolved fields obtained with the approximate reconstruction using moments (ARM) agree better with the filtered direct numerical simulation since the level of reconstruction in the modeled subfilter fluxes is increased.     

Scalar dissipation rates in isothermal and reactive turbulent opposed-jets: 1-D-Raman/Rayleigh experiments supported by LES

Isothermal and reactive turbulent opposed flows are presented, which are appropriate to test the applicability and performance of models for turbulence, mixing, chemical reaction, and turbulence-chemistry interaction. Transient flow and scalar fields are measured using laser Doppler velocimetry and one-dimensionally resolved Raman/Rayleigh spectroscopy. Aside of statistical moments of temperature, mean species, and velocity components, scalar dissipation rate across the mixing and reaction layer is determined on a single-shot base. Using large eddy simulation in connection with a steady flamelet model, it is shown how numerical data can serve to estimate the influence of experimental noise upon a measured quantity, such as scalar dissipation. As a key result, it is shown that an increase in scalar rate of dissipation by chemical reactions is caused by a significant increase in the mixture fraction diffusivity, which outweighs the decrease in mixture fraction gradients. In mixture fraction space, local maxima of scalar dissipation rate are found on the rich side, which cannot be correctly reproduced by the steady flamelet model assuming equal species diffusivity. Furthermore, the impact of experimental noise on conditional probability density functions of scalar dissipation rate is shown (exemplary) to lead to erroneous conclusions from experimental data.     

Simulation of spray development and turbulent combustion processes in low and high speed diesel engines by the CMC-ISR model

Simulation is performed to analyse the characteristics of turbulent spray combustion in conventional low and high speed diesel engine conditions. Turbulence–chemistry interaction is resolved by the Conditional Moment Closure (CMC) model in the spatially integrated form of an Incompletely Stirred Reactor (ISR). After validation against measured pressure traces, characteristic length and time scales and dimensionless numbers are estimated at the locations of sequentially injected fuel groups. Conditional flame structures are calculated for sequentially evaporated fuel groups to consider different available periods for ignition chemistry. Injection overlaps the combustion period in the high rpm engine, while most combustion occurs after injection and evaporation are complete in the low rpm engine. Ignition occurs in rich premixture with the initial peak temperature at the equivalence ratio around 2–4 as observed in Dec [2]. It corresponds to the most reactive mixture fraction of the minimum ignition delay for the given mixture states. Combustion proceeds to lean and rich sides in the mixture fraction space as a diffusion process by turbulence. The mean scalar dissipation rates (SDRs) are lower than the extinction limit to show stability of diffusion flames throughout the combustion period.