LES/CMC Simulations of Swirl-Stabilised Ethanol Spray Flames Approaching Blow-Off

Large Eddy Simulations (LES) with the Conditional Moment Closure (CMC) combustion model of swirling ethanol spray flames have been performed in conditions close to blow-off for which a wide database of experimental measurements is available for both flame and spray characterization. The solution of CMC equations exploits a three-dimensional unstructured code with a first order closure for chemical source terms. It is shown that LES/CMC is able to properly capture the flame structure at different conditions and agrees reasonably well with the measurements both in terms of mean flame shape and dynamic behaviour of the flame evaluated in terms of local extinctions and statistics of the lift-off height. Experimental measurements of the overall (liquid plus gaseous) mixture fraction, performed using the Laser-Induced Breakdown Spectroscopy technique, are also included allowing further assessment and validation of the numerical method. The sensitivity of the simulation results to the various boundary conditions is discussed.     

One-Dimensional Modeling of Turbulent Premixed Jet Flames - Comparison to DNS

The one-dimensional turbulence (ODT) model, formulated in an Eulerian reference frame, is applied to a temporally-evolving premixed turbulent hydrogen plane-jet flame and results are compared with direct numerical simulation (DNS) data. This is the first published study to perform direct comparisons of ODT to DNS for premixed flames. The ODT model solves the full set of conservation equations for mass, momentum, energy, and species on a one-dimensional domain corresponding to the transverse jet direction. The effects of turbulent mixing are modeled via a stochastic process, while the full range of diffusive-reactive length and time scales are resolved directly on the one-dimensional domain. A detailed chemical mechanism for hydrogen combustion consisting of 9 species and 21 reactions and a mixture-averaged transport model are used (consistent with the DNS). Cases with two different Damköhler numbers are considered and comparisons between the ODT and DNS data are shown with respect to flow dynamics and thermochemistry. The ODT compared favorably with the DNS in terms of the overall entrainment as judged by the streamwise velocity profile and in terms of local flamelet structure as judged by progress-variable conditional reaction and scalar dissipation rates. While the ODT agreed qualitatively with the overall flame evolution, the net fuel consumption rate was somewhat over-predicted for a brief early period and under-predicted later on, leading to an overly long flame burnout time. It was demonstrated that adjusting a parameter controlling the selection of large eddies improved the prediction of the peak fuel consumption rate and overall reaction progress but worsened the prediction of jet entrainment. An analysis of the 1D nature of ODT is presented that suggests the FSD in ODT needs to be much higher than the FSD in the DNS in order to achieve the same overall burning rate, suggesting that the FSD is under-predicted by a significant fraction. While the success of the ODT in reproducing many of the salient features of nonpremixed flames has been demonstrated, the current study suggests that improvements are needed when applied to premixed flames. It is also important to note that the DNS required approximately 40×10⁶ CPU hours while the ODT required approximately 10³ CPU hours.     

Joint PDF Closure of Turbulent Premixed Flames

In this paper, a novel model for turbulent premixed combustion in the corrugated flamelet regime is presented, which is based on transporting a joint probability density function (PDF) of velocity, turbulence frequency and a scalar vector. Due to the high dimensionality of the corresponding sample space, the PDF equation is solved with a Monte-Carlo method, where individual fluid elements are represented by computational particles. Unlike in most other PDF methods, the source term not only describes reaction rates, but accounts for “ignition” of reactive unburnt fluid elements due to propagating embedded quasi laminar flames within a turbulent flame brush. Unperturbed embedded flame structures and a constant laminar flame speed (as expected in the corrugated flamelet regime) are assumed. The probability for an individual particle to “ignite” during a time step is calculated based on an estimate of the mean flame surface density (FSD), latter gets transported by the PDF method. Whereas this model concept has recently been published [21], here, a new model to account for local production and dissipation of the FSD is proposed. The following particle properties are introduced: a flag indicating whether a particle represents the unburnt mixture; a flame residence time, which allows to resolve the embedded quasi laminar flame structure; and a flag indicating whether the flame residence time lies within a specified range. Latter is used to transport the FSD, but to account for flame stretching, curvature effects, collapse and cusp formation, a mixing model for the residence time is employed. The same mixing model also accounts for molecular mixing of the products with a co-flow. To validate the proposed PDF model, simulation results of three piloted methane-air Bunsen flames are compared with experimental data and very good agreement is observed.     

Turbulence-Driven Blowout Instabilities of Premixed Bluff-Body Flames

Bluff-body flame instabilities are experimentally investigated under varying turbulence conditions during lean blowout. For all turbulence conditions, the blowout process is induced through a temporal reduction of the fuel flow rate to capture the flame-flow dynamics approaching blowout. Simultaneous high-speed particle image velocimetry (PIV), stereoscopic PIV, and C₂*/CH* chemiluminescence imaging are employed, along with an independent CH* imaging system, to capture flame-flow instabilities. Proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) techniques are used to identify prominent flame oscillations and evaluate recurring spatiotemporal modes during blowout. The results reveal that the dominant flame oscillations and wrinkling characteristics are directly dependent on the turbulence conditions in the combustor. Specifically, the flame-flow oscillations are strongly coupled with the integral length scales, which were able to collapse the oscillation frequencies to a unified value. The turbulence-driven flame-flow oscillations are shown to largely impact the magnitude, temporal evolution, and oscillatory behavior of the flame strain rate. As the turbulence intensity is increased, the oscillation of the flame strain rate increases in frequency, making it more likely for localized extinctions to occur. Additionally, the magnitude of the flame strain rate increases at high turbulence intensities and accelerates the lean blowout process.

     

A-Priori Validation of Scalar Dissipation Rate Models for Turbulent Non-Premixed Flames

The modelling of scalar dissipation rate in conditional methods for large-eddy simulations is investigated based on a priori direct numerical simulation analysis using a dataset representing an igniting non-premixed planar jet flame. The main objective is to provide a comprehensive assessment of models typically used for large-eddy simulations of non-premixed turbulent flames with the Conditional Moment Closure combustion model. The linear relaxation model gives a good estimate of the Favre-filtered scalar dissipation rate throughout the ignition with a value of the related constant close to the one deduced from theoretical arguments. Such value of the constant is one order of magnitude higher than typical values used in Reynolds-averaged approaches. The amplitude mapping closure model provides a satisfactory estimate of the conditionally filtered scalar dissipation rate even in flows characterised by shear driven turbulence and strong density variation.

     

Simulations of Turbulent Non-Premixed Counterflow Flames with First-Order Conditional Moment Closure

Simulations of turbulent CH₄-air counterflow flames are presented, obtained in terms of zero and two-dimensional first-order Conditional Moment Closure (CMC) to study the flame structure and extinction limits. The CMC equation with detailed chemistry is solved without the need for operator splitting, while the accompanying flow field is determined using a commercial CFD software employing a Reynolds stress turbulence model and additional transport equations for the turbulent scalar flux and for the mean scalar dissipation rate. Two detailed chemical mechanisms and different conditional scalar dissipation rate models have been examined and small differences were found.The first-order CMC captures the overall structure of the counterflow flame accurately for the unconditional averages. The calculated conditional averages behave as if the scalar dissipation rate were under-predicted, although a comparison with measurement of the conditional scalar dissipation rate is reasonable. The calculated extinction velocity is found to be much higher than the experimental value, but the trend of increasing extinction velocity with air dilution of the fuel stream is captured well. The discrepancies with the data are mostly attributed to the neglect of conditional fluctuations.     

On the Flame-generated Vorticity Dynamics of Bluff-body-stabilized Premixed Flames

This investigation considers the dynamics of flame-generated vorticity for a premixed, submerged bluff-body stabilized flame. Experimentation characterizes the far-field region in particular with a level of detail not previously afforded to this type of flow. Simultaneous particle imaging velocimetry (PIV), Mie scattering and CH ^? chemiluminescence are used to obtain velocity fields and flame location. Mean static pressure measurements at the combustion chamber wall capture the pressure field. Analysis of the flame fronts in relation to the mean velocity and vorticity fields provides useful insight into the interaction of the flame and the flow. The unique nature of the velocity and vorticity fields and their effect on downstream flame structures are explained by the baroclinic torque generation of vorticity. The coupling that exists among pressure, heat release, and baroclinic generation is acknowledged and will influence strategies for control of the baroclinic mechanism.     

Topology and Brush Thickness of Turbulent Premixed V-shaped Flames

Topology and brush thickness of turbulent premixed V-shaped flames were investigated using Mie scattering and Particle Image Velocimetry techniques. Mean bulk flow velocities of 4.0, 6.2, and 8.3 m/s along with two fuel-air equivalence ratios of 0.6 and 0.7 were tested in the experiments. Using a novel experimental turbulence generating apparatus, three turbulence intensities of approximately 2 %, 6 %, and 17 % were tested in the experiments. The results show that topology of the flame front is significantly altered by changing the turbulence intensity. Specifically, at relatively small turbulence intensities, the flame fronts feature wrinkles which are symmetric with respect to the vertical axis. At moderate values of turbulence intensities, the flame fronts form cusps. The formation of cusps is more pronounced at large mean bulk flow velocities. The results associated with relatively large turbulence intensity show that flame surfaces feature: mushroom-shaped structures, freely propagating sub-flames, pocket formation, localized extinction, and horn-shaped structures. Analysis of the results show that the flame brush thickness follows a linear correlation with the root-mean-square of the flame front position. The correlation is in agreement with the results of past experimental investigations associated with moderately turbulent premixed V-shaped flames, and holds for the range of turbulence conditions tested. This suggests that the underlying mechanism associated with the dynamics of moderately turbulent premixed V-shaped flames proposed in past studies can potentially be valid for the the wide range of turbulence conditions examined in the present investigation.     

Turbulence Statistics and Scalar Transport in Highly-Sheared Premixed Flames

This paper describes recent progress in the analysis of the nature of turbulent premixed flames stabilised behind an axisymmetric baffle which are of fundamental interest in the development of new and cleaner combustion systems. The work includes the use of laser-based diagnostics for velocity and temperature measurements, which are extended to the analysis of turbulence statistics, including the energy spectrum and typical length scales in a reacting shear layer. The results provided experimental evidence of the extension of the flamelet regime beyond the Klimov--Williams criterion. Arguments based on the shape of the weighted-joint-probability distributions of axial velocity and temperature fluctuations show that the counter-gradient nature of heat flux is associated with the preferential deceleration of products of combustion in relation to the cold reactants.     

A Comparison of Strategies for Direct Numerical Simulation of Turbulence Chemistry Interaction in Generic Planar Turbulent Premixed Flames

Three different methods to introduce turbulence in the computational domain of Direct Numerical Simulations (DNS) of statistically planar turbulent premixed flame configurations have been reviewed and their advantages and disadvantages in terms of run time, natural flame development, control of turbulence parameters and convergence of statistics extracted from the simulations have been discussed in detail. It has been found that there is no method, which is clearly superior to the other two alternative methods. An analysis has been performed to explain why Lundgren’s physical space linear forcing results in an integral length scale which is, independent of the Reynolds number, a constant fraction of the domain size. Furthermore, an evolution equation for the integral length scale has been derived, and a scaling analysis of its terms has been performed to explain the evolution of the integral length scale in the context of Lundgren’s physical space linear forcing. Finally, a modification to Lundgren’s forcing approach has been suggested which ensures that the integral length scale settles to a predetermined value so that DNS of statistically planar turbulent premixed flames with physical space forcing can be conducted for prescribed values of Damköhler and Karlovitz numbers.     

Subgrid Modeling of Turbulent Premixed Flames in the Flamelet Regime

Large eddy simulation (LES) models for flamelet combustion are analyzed by simulating premixed flames in turbulent stagnation zones. ALES approach based on subgrid implementation of the linear eddy model(LEM) is compared with a more conventional approach based on the estimation of the turbulent burning rate. The effects of subgrid turbulence are modeled within the subgrid domain in the LEM-LES approach and the advection (transport between LES cells) of scalars is modeled using a volume-of-fluid (VOF) Lagrangian front tracking scheme. The ability of the VOF scheme to track the flame as a thin front on the LES grid is demonstrated. The combined LEM-LES methodology is shown to be well suited for modeling premixed flamelet combustion. The geometric characteristics of the flame surfaces, their effects on resolved fluid motion and flame-turbulence interactions are well predicted by the LEM-LES approach. It is established here that local laminar propagation of the flamelets needs to be resolved in addition to the accurate estimation of the turbulent reaction rate. Some key differences between LEM-LES and the conventional approach(es) are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Experimental Characterization of Gelled Jet A1 Spray Flames

Gelled propellants provide energetic performance similar to conventional liquid propellants and safety during storage and handling like a solid propellant. Experiments on unconfined gelled Jet A1 spray flames and the comparison with ungelled spray flames are reported for the first time in this paper in terms of the global features, burning regimes, stability limits, visible flame height, emission spectra, natural luminosity, and CH ^? chemiluminescence. Propellants were atomized by an internally impinging two-fluid atomizer, developed specifically for efficient atomization of non-Newtonian gels. Swirling and non-swirling spray flames were successfully stabilized on a burner incorporating bluff body and annular jet of combustion air over a wide range of operating parameters. Structural features of the atomizer impart high momentum to the (central) spray jet, such that the recirculation zone could be penetrated under all conditions. Long-exposure smoke and high-speed visualizations were employed to study cold flow structures and droplet-vortex interactions. Short-exposure direct and backlit imaging were used to observe global features of spray flames. Stability limits and visible flame heights were mapped for different thermal inputs, swirl numbers, and flow rates of atomizing and combustion air jets. Non-swirling stable anchored, partially blown off, and neck-blown off flames were observed. Lifted, and a transition regime, in which the flame could burn in stable and lifted mode repetitively, were observed for the swirling flames. Interactions between central and annular jets are important in these regimes, determining flame shape, symmetry, and flame height. Jet-like propagation zone determines the flame height through its dependence on momentum of spray jet. The length of this zone is affected by variations in thermal input, gas-liquid ratio, and air-fuel ratio. The gelled Jet A1 flames are remarkably shorter despite having a larger average droplet size than ungelled Jet A1. This experimental observation directly supports theoretical predictions reported in literature. These flames are more luminous than ungelled Jet A1, especially at the base and the neck regions. While, majority of the heat is released in the jet-like propagation zone for both the flames, significant heat is released in the neck zone of ungelled Jet A1 spray flame in comparison to ungelled Jet A1 spray flame due to intense turbulence and smaller droplet size.     

Transient Behavior of Turbulent Scalar Transport in Premixed Flames

Transition from gradient to countergradient scalar transport in a statistically planar, one-dimensional, developing, premixed turbulent flame is studied both theoretically and numerically. A simple criterion of the transition referred to is derived from the balance equation for the combustion progress variable, with the criterion highlighting an important role played by flame development. A balance equation for the difference in velocities $\bar{u}_b$ and $\bar{u}_u$ conditioned on burned and unburned mixture, respectively, is numerically integrated. Both analytical and computed results show that; (1) The flux $\overline{\rho u'' c''}$ is gradient during an early stage of flame development followed by transition to countergradient scalar transport at certain instant t _tr . (2) The transition time is increased when turbulence length scale L is increased or when the laminar flame speed S _L and/or the density ratio are decreased. (3) The transition time normalized using the turbulence time scale is increased by u??. Moreover, the numerical simulations have shown that the transition time is increased by u?? if a ratio of u??/S _L is not large. This dependence of t _tr on u?? is substantially affected by (i) the mean pressure gradient induced within the flame due to heat release and (ii) by the damping effect of combustion on the growth rate of mean flame brush thickness. The reasonable qualitative agreement between the computed trends and available experimental and DNS data, as well as the agreement between the computed trends and the present theoretical results, lends further support to the conditioned balance equation used in the present work.     

Resolution Requirements in Stochastic Field Simulation of Turbulent Premixed Flames

The spatial resolution requirements of the Stochastic Fields probability density function approach are investigated in the context of turbulent premixed combustion simulation. The Stochastic Fields approach is an attractive way to implement a transported Probability Density Function modelling framework into Large Eddy Simulations of turbulent combustion. In premixed combustion LES, the numerical grid should resolve flame-like structures that arise from solution of the Stochastic Fields equation. Through analysis of Stochastic Fields simulations of a freely-propagating planar turbulent premixed flame, it is shown that the flame-like structures in the Stochastic Fields simulations can be orders of magnitude narrower than the LES filter length scale. The under-resolution is worst for low Karlovitz number combustion, where the thickness of the Stochastic Fields flame structures is on the order of the laminar flame thickness. The effect of resolution on LES predictions is then assessed by performing LES of a laboratory Bunsen flame and comparing the effect of refining the grid spacing and filter length scale independently. The usual practice of setting the LES filter length scale equal to grid spacing leads to severe under-resolution and numerical thickening of the flame, and to substantial error in the turbulent flame speed. The numerical resolution required for accurate solution of the Stochastic Fields equations is prohibitive for many practical applications involving high-pressure premixed combustion. This motivates development of a Thickened Stochastic Fields approach (Picciani et al. Flow Turbul. Combust. X, YYY (2018) in order to ensure the numerical accuracy of Stochastic Fields simulations.     

Influence of Differential Diffusion on Super-Equilibrium Temperature in Turbulent Non-Premixed Hydrogen/Air Flames

In this paper we wish to investigate the occurrence of super-equilibrium temperature values, observed in many experimental configurations. We would like to understand the origin of this phenomenon. Previous authors have already shown that differential diffusion can lead to considerable changes in the temperature field and we would like to build on top of this observation. We investigate numerically super-equilibrium combustion by considering both laminar counter-flow and turbulent diluted hydrogen/air diffusion flames. These turbulent flames are computed using direct numerical simulations (DNS). A detailed reaction mechanism is employed and the transport properties are modeled using multicomponent diffusion velocities, including the Soret effect. Analyzing these results we introduce three complementary parameters (dilution-free mixture fraction, dilution excess and local enthalpy) to describe the local combustion conditions. Introducing a measure of dilution separately from the mixture fraction is necessary for a proper analysis. Using this set of parameters it becomes possible to explain super-equilibrium temperature levels as a consequence of differential diffusion.     

A Priori Evaluation of the Laminar Flamelet Decomposition Model for Turbulent Premixed Flames using DNS Data

Laminar flamelet decomposition (LFD) is a dynamic approach for modelling sub-filter scale turbulence-chemistry interactions in Large-Eddy Simulations using a stretched flamelet library. In this work, the performance of the LFD model – that was previously used only in non-premixed combustion—is investigated a priori for premixed combustion using positively-strained flamelets in the reactant-to-product configuration. For this purpose, a DNS database of methane-air premixed flames is utilized. The flames are propagating in a rectangular box under homogeneous isotropic turbulence conditions over a wide range of Karlovitz numbers. The results show that the LFD model can correctly account for the sub-filter scale turbulence-chemistry interactions to predict the filtered reaction rates and the filtered scalar field, provided that turbulent and laminar mixing are well predicted. The deviations from the DNS results are attributed to the shortcomings of the strained flamelet library and the non-flamelet effects. Finally, the LFD results are compared with a different sub-filter scale model using the same strained flamlelet library, and the relative performances of the two models are discussed.—

     

Effects of Lewis Number on Scalar Variance Transport in Premixed Flames

The influences of differential diffusion rates of heat and mass on the transport of the variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been studied using three-dimensional simplified chemistry based Direct Numerical Simulation (DNS) data of statistically planar turbulent premixed flames with global Lewis number ranging from Le?= 0.34 to 1.2. The Lewis number effects on the statistical behaviours of the various terms of the transport equations of variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been analysed in the context of Reynolds Averaged Navier Stokes (RANS) simulations. It has been found that the turbulent fluxes of the progress variable and temperature variances exhibit counter-gradient transport for the flames with Lewis number significantly smaller than unity whereas the extent of this counter-gradient transport is found to decrease with increasing Lewis number. The Lewis number is also shown to have significant influences on the magnitudes of the chemical reaction and scalar dissipation rate contributions to the scalar variance transport. The modelling of the unclosed terms in the scalar variance equations for the non-unity Lewis number flames have been discussed in detail. The performances of the existing models for the unclosed terms are assessed based on a-priori analysis of DNS data. Based on the present analysis, new models for the unclosed terms of the active scalar variance transport equations are proposed, whenever necessary, which are shown to satisfactorily capture the behaviours of unclosed terms for all the flames considered in this study.     

Large Eddy Simulation of Premixed Turbulent Flames Using the Probability Density Function Approach

In the present study a Large Eddy Simulation and Filtered Density Function model is applied to three premixed piloted turbulent methane flames at different Reynolds Numbers using the Eulerian stochastic fields approach. The model is able to reproduce the flame structure and flow characteristics with a low number of fields (between 4 and 16 fields). The results show a good agreement with experimental data with the same closures employed in non-premixed combustion without any adjustment for combustion regime. The effect of heat release on the flow field is captured correctly. A wide range of sensitivity studies is carried out, including the number of fields, the chemical mechanism, differential diffusion effects and micro-mixing closures. The present work shows that premixed combustion (at least in the conditions under study) can be modelled using LES-PDF methods.. Finally, the ability of the model to predict flame quenching is studied. The model can accurate capture the conditions at which combustion is not sustainable and large pockets of extinction appear.     

Leading-Edge Reaction Zones in Lifted-Jet Gas and Spray Flames

An investigation of the leading edge characteristics in lifted turbulent methane-air (gaseous) and ethanol-air (spray) diffusion flames is presented. Both combustion systems consist of a central nonpremixed fuel jet surrounded by low-speed air co-flow. Non-intrusive laser-based diagnostic techniques have been applied to each system to provide information regarding the behavior of the combustion structures and turbulent flow field in the regions of flame stabilization. Simultaneous sequential CH-PLIF/particle image velocimetry and CH-PLIF/Rayleigh scattering measurements are presented for the lifted gaseous flame. The CH-PLIF data for the lifted gas flame reveals the role that ``leading-edge' combustion plays as the stabilization mechanism in gaseous diffusion flames. This phenomenon, characterized by a fuel-lean premixed flame branch protruding radially outward at the flame base, permits partially premixed flame propagation against the incoming flow field. In contrast, the leading edge of the ethanol spray flame, examined using single-shot OH-PLIF imaging and smoke-based flow visualization, does not exhibit the same variety of leading-edge combustion structure, but instead develops a dual reaction zone structure as the liftoff height increases. This dual structure is a result of the partial evaporation (hence partial premixing) of the polydisperse spray and the enhanced rate of air entrainment with increased liftoff height (due to co-flow). The flame stabilizes in a region of the spray, near the edge, occupied by small fuel droplets and characterized by intense mixing due to the presence of turbulent structures. This revised version was published online in July 2006 with corrections to the Cover Date.