Lagrangian conditional moment closure model with flame group interaction for lifted turbulent spray jet flames

A new Lagrangian conditional moment closure (CMC) model is developed for multiple Lagrangian groups of sequentially evaporating fuel in turbulent spray combustion. Flame group interaction is taken into account as premixed combustion by the eddy breakup (EBU) model in terms of the probability of finding flame groups in the burned and the unburned state. Evaporation source terms are included in the two phase conditional transport equations, although they turn out to have negligible influence on the mean temperature field during combustion. The Lagrangian CMC model is implemented in OpenFOAM [1 H.G. Weller, G. Tabor, H. Jasak, and C. Fureby, A tensorial approach to computational continuum mechanics using object-oriented techniques, Comput. Phys. 12 (1998), pp. 620–631.[Crossref] [Google Scholar]] and validated for test cases in the Engine Combustion Network (ECN) [2 Engine Combustion Network. (2011). Available at http://www.sandia.gov/ecn. [Google Scholar],3 L.M. Pickett, C.L. Genzale, G. Bruneaux, L. Malbec, L. Hermant, C. Christian, and J. Schramm, Comparison of diesel spray combustion in different high-temperature, high-pressure facilities, SAE technical paper 2010-01-2106 (2010). [Google Scholar]]. Similar ignition delays and lift-off lengths are predicted by the incompletely stirred reactor (ISR) and the Eulerian CMC models due to relatively uniform conditional flame structure in the domain. The improved Lagrangian CMC model shows no abrupt reaction or oscillatory behaviour with an appropriate model constant K and gives results in better agreement with measurements lying between the predictions by ISR and Lagrangian CMC without flame group interaction.     

Numerical study of transient evolution of lifted jet flames: partially premixed flame propagation and influence of physical dimensions

Three-dimensional (3D) unsteady Reynolds-averaged Navier–Stokes simulations of a spark-ignited turbulent methane/air jet flame evolving from ignition to stabilisation are conducted for different jet velocities. A partially premixed combustion model is used involving a correlated joint probability density function and both premixed and non-premixed combustion mode contributions. The 3D simulation results for the temporal evolution of the flame's leading edge are compared with previous two-dimensional (2D) results and experimental data. The comparison shows that the final stabilised flame lift-off height is well predicted by both 2D and 3D computations. However, the transient evolution of the flame's leading edge computed from 3D simulation agrees reasonably well with experiment, whereas evident discrepancies were found in the previous 2D study. This difference suggests that the third physical dimension plays an important role during the flame transient evolution process. The flame brush's leading edge displacement speed resulting from reaction, normal and tangential diffusion processes are studied at different typical stages after ignition in order to understand the effect of the third physical dimension further. Substantial differences are found for the reaction and normal diffusion components between 2D and 3D simulations especially in the initial propagation stage. The evolution of reaction progress variable scalar gradients and its interaction with the flow and mixing field in the 3D physical space have an important effect on the flame's leading edge propagation.     

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.     

Simulations of sooting turbulent jet flames using a hybrid flamelet/stochastic Eulerian field method

The stochastic Eulerian field method is applied to simulate 12 turbulent C₁?C₃ hydrocarbon jet diffusion flames covering a wide range of Reynolds numbers and fuel sooting propensities. The joint scalar probability density function (PDF) is a function of the mixture fraction, enthalpy defect, scalar dissipation rate and representative soot properties. Soot production is modelled by a semi-empirical acetylene/benzene-based soot model. Spectral gas and soot radiation is modelled using a wide-band correlated-k model. Emission turbulent radiation interactions (TRIs) are taken into account by means of the PDF method, whereas absorption TRIs are modelled using the optically thin fluctuation approximation. Model predictions are found to be in reasonable agreement with experimental data in terms of flame structure, soot quantities and radiative loss. Mean soot volume fractions are predicted within a factor of two of the experiments whereas radiant fractions and peaks of wall radiative fluxes are within 20%. The study also aims to assess approximate radiative models, namely the optically thin approximation (OTA) and grey medium approximation. These approximations affect significantly the radiative loss and should be avoided if accurate predictions of the radiative flux are desired. At atmospheric pressure, the relative errors that they produced on the peaks of temperature and soot volume fraction are within both experimental and model uncertainties. However, these discrepancies are found to increase with pressure, suggesting that spectral models describing properly the self-absorption should be considered at over-atmospheric pressure.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.     

Large Eddy Simulation of flame stabilisation dynamics and vortex control in a lifted H2/N2 jet flame

Flame stabilisation in (highly) preheated mixture is common in several industrial applications. When the reactants are injected separately in the device (usually at high-speed), the flame is lifted so that the fuel and oxidant first mix to give an ignitable mixture. If the temperature of the mixture is adequate, it auto-ignites stabilizing the flame. Here we focus on an academic lifted jet flame and Large Eddy Simulation (LES) is used to capture the flame and auto-ignition dynamics. Comparisons with experimental data show that LES simulates accurately high OH fluctuation levels at the stabilisation location. The vortex dynamics linked to these fluctuations is analyzed and it is found that small scale coherent structures play a vital role in the auto-ignition process. These structures are axial vorticity tubes (braids) and are located relatively far (in the radial direction) from the shear-layer. As a consequence, the lift-off height varies dramatically in time leading to OH fluctuations of the order of the mean OH concentration. This scenario is monitored in the compositional space highlighting the simultaneous evolution of OH, HO ₂ and temperature. Further, different strategies for open-loop control of the flame lift-off height are tested. In order to anchor the flame at different positions downstream of the nozzle, the vortex dynamics in the shear-layer was modified. Promoting successively vortex ring and braids, the auto-ignition region was moved significantly. In particular, modified nozzle geometries impacted the formation of braids and ensured a good premixing very close to the nozzle. As a consequence, it was possible to reduce significantly the lift-off height and stabilise the flame few diameters downstream of the nozzle.     

A priori analysis of sub-grid variance of a reactive scalar using DNS data of high Ka flames

Direct numerical simulations (DNS) of low and high Karlovitz number (Ka) flames are analysed to investigate the behaviour of the reactive scalar sub-grid scale (SGS) variance in premixed combustion under a wide range of combustion conditions (regimes). An order of magnitude analysis is performed to assess the importance of various terms in the variance evolution equation and the analysis is validated using the DNS results. This analysis sheds light on the relative behaviour among turbulent transport and production, scalar dissipation and chemical processes involved in the evolution of the SGS variance at different Ka. The common expectation is that the variance equation shifts from a reaction-dissipation balance at low Ka to a production–dissipation balance at high Ka with diminishing reaction contribution. However, in large eddy simulation (LES), a high Ka alone does not make the reaction term negligible, as the relative importance of the reaction term has a concurrent increase with filter size. The filter size can be relatively large compared with the Kolmogorov length scale in practical LES of high Ka flames, and as a consequence a reaction–production–dissipation balance may prevail in the variance equation even in a high Ka configuration, and this possibility is quantified using the DNS analysis in this work. This has implications from modelling perspectives, and therefore two commonly used closures in LES for the SGS scalar dissipation rate are investigated a priori to estimate the importance of the above balance in LES modelling. The results are explained to highlight the interplay among turbulence, chemistry and dissipation processes as a function of Ka.     

RANS predictions of turbulent diffusion flames: comparison of a reactor and a flamelet combustion model to the well stirred approach

The flame stabilisation process in turbulent non-premixed flames is not fully understood and several models have been developed to describe the turbulence–chemistry interaction. This work compares the performance of the multiple Representative Interactive Flamelet (mRIF) model, the Volume Reactor Fraction Model (VRFM), and the Well-Stirred reactor (WS) model in describing such flames. The predicted ignition delay and flame lift-off length of n-heptane sprays are compared to experimental results published within the Engine Combustion Network (ECN). All of the models predict the trend of ignition delay reasonably well. At a low gas pressure (42 bar) the ignition delay is overpredicted compared to the experimental data, but the difference between the models is not significant. However, the predicted lift-off lengths differ. At high pressure (87 bar) the difference between the models is small. All models slightly underpredict the lift-off length compared to the experimental data. At low gas pressure (42 bar) the mRIF model gives the best results. The VRFM and WS models predict excessively short lift-off lengths, but the VRFM model gives better results than the WS model. The flame structures of the models are also compared. The WS model and the VRFM model yield a well defined flame stabilisation point whereas the mRIF model does not. The flame of the mRIF model is more diffuse and the model is not able to predict flame propagation. All models were able to predict the experimental trends in lift-off and ignition delay, but certain differences between them are demonstrated.     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

A DNS assessment of linear relations between filtered reaction rate,flame surface density,and scalar dissipation rate in a weakly turbulent premixed flame

Linear relations between (i) filtered reaction rate and filtered flame surface density (FSD) and (ii) filtered reaction rate and filtered scalar dissipation rate (SDR), which are widely used in Large Eddy Simulation (LES) research into premixed turbulent combustion, are examined by processing DNS data obtained from a statistically 1D planar flame under weakly turbulent conditions that are most favourable for the two approaches (flamelet combustion regime, single-step chemistry, equidiffusive mixture, adiabatic burner, and low Mach number). The analysis well supports the former approach provided that the filtered reaction rate is combined with filtered molecular transport term. In such a case, both the RANS and LES FSD approaches are based on local relations valid within weakly perturbed flamelets. Accordingly, simply recasting RANS expressions to a filtered form works well. On the contrary, while the FSD and SDR approaches appear to be basically similar at first glance, the analysis does not support the latter one, but shows that a ratio of the filtered reaction rate to the filtered SDR is strongly scattered within the studied flame brush, with its conditionally mean value varying significantly with Favre-filtered combustion progress variable. As argued in the paper, these limitations of the LES SDR approach stem from the fact that it is based on a relation valid after integration over weakly perturbed flamelets, but this relation does not hold locally within such flamelets. Consequently, when a sufficiently small filter is applied to instantaneous fields, the filter may contain only a part of the local flamelet, whereas the linear relation holds solely for the entire flamelet and may not hold within the filtered flamelet volume. Thus, the present study implies that straightforwardly recasting well-established RANS equations to a filtered form is a flawed approach if the equations are based on integral features of local burning.