Statistical Analysis and a-priori Modelling of Flame Surface Density Transport in Turbulent Stratified Flames: A Direct Numerical Simulation Study

Statistically planar turbulent premixed and stratified flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios <???> = 0.7 and <???> = 1.0 using three-dimensional simplified chemistry based Direct Numerical Simulations (DNS). The simulation parameters are chosen such that the thin reaction zones regime combustion is realised in all cases and a random bi-modal distribution of equivalence ratio ? is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity for stratified flames. The modelling of the unclosed terms (i.e. the turbulent transport term T ₁, the tangential strain rate term T ₂, the propagation term T ₃, and the curvature term T ₄) of the generalised FSD transport equation has been addressed in the context of RANS simulations. It has been found that the turbulent transport term T ₁ remains small in comparison to the leading order contributions of the tangential strain rate and curvature terms (i.e. T ₂ and T ₄, respectively) in the globally stoichiometric cases, but T ₁ begins to play a more important role in the globally fuel-lean cases. The strain rate term T ₂ remains positive throughout the flame brush and acts as a leading order source term for all the flames considered in this analysis. It is has been found that the magnitude of T ₂ decreases with decreasing root-mean-square velocity fluctuations u ^′ (<???>) for a given value of <???> (u ^′). The contribution of the propagation term T ₃ remains generally positive towards the unburned gas side of the flame brush but assumes generally negative values towards the burned gas side of the flame brush. Moreover, whilst the order of magnitude of the propagation term T ₃ is comparable in all cases, T ₃ remains small in comparison to the leading order contributors (i.e. T ₂ and T ₄) in the globally stoichiometric cases however it plays a more important role in the globally fuel-lean cases. The curvature term T ₄ acts as a leading order sink term in all cases except towards the unburned gas side of the flame brush in low u ^′ globally stoichiometric (i.e. <???> = 1.0) flames. Furthermore, it has been demonstrated that the magnitude of T ₄ decreases with decreasing u ^′ (<???>) for a given value of <???> (u ^′). Appropriate model expressions have been identified for T ₁, T ₂, T ₃ and T ₄ based on an a-priori analysis of the DNS data.     

Presumed Joint-PDF Modelling for Turbulent Stratified Flames

Flow, Turbulence and Combustion - Turbulence–chemistry interaction models such as the conditional moment closure and various flamelet models require a presumed Probability Density Function...     

Statistical Analysis of Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Statistically planar turbulent premixed and partially premixed flames for different initial turbulence intensities are simulated for global equivalence ratios ??>?=?0.7 and ??>?=?1.0 using three-dimensional Direct Numerical Simulations (DNS) with simplified chemistry. For the simulations of partially premixed flames, a random distribution of equivalence ratio following a bimodal distribution of equivalence ratio is introduced in the unburned reactants ahead of the flame. The simulation parameters in all of the cases were chosen such that the combustion situation belongs to the thin reaction zones regime. The DNS data has been used to analyse the behaviour of the dissipation rate transports of both active and passive scalars (i.e. the fuel mass fraction Y _F and the mixture fraction ξ) in the context of Reynolds Averaged Navier–Stokes (RANS) simulations. The behaviours of the unclosed terms of the Favre averaged scalar dissipation rates of fuel mass fraction and mixture fraction (i.e. \(\widetilde {\varepsilon }_Y =\overline {\rho D\nabla Y_F^{\prime \prime } \cdot \nabla Y_F^{\prime \prime } } /\overline{\rho }\) and \(\widetilde {\varepsilon }_\xi =\overline {\rho D\nabla \xi ^{\prime \prime }\cdot \nabla \xi ^{\prime \prime }} /\overline {\rho })\) transport equations have been analysed in detail. In the case of the \(\widetilde {\varepsilon }_Y \) transport, it has been observed that the turbulent transport term of scalar dissipation rate remains small throughout the flame brush whereas the terms due to density variation, scalar–turbulence interaction, reaction rate and molecular dissipation remain the leading order contributors. The term arising due to density variation remains positive throughout the flame brush and the combined contribution of the reaction and molecular dissipation to the \(\widetilde {\varepsilon }_Y \) transport remains negative throughout the flame brush in all cases. However, the behaviour of scalar–turbulence interaction term of the \(\widetilde {\varepsilon }_Y \) transport equation is significantly affected by the relative strengths of turbulent straining and the straining due to chemical heat release. In the case of the \(\widetilde {\varepsilon }_\xi \) transport, the turbulent transport term remains small throughout the flame brush and the density variation term is found to be negligible in all cases, whilst the reaction rate term is exactly zero. The scalar–turbulence interaction term and molecular dissipation term remain the leading order contributors to the \(\widetilde {\varepsilon }_\xi \) transport throughout the flame brush in all cases that have been analysed in the present study. Performances of existing models for the unclosed terms of the transport equations of \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) are assessed with respect to the corresponding quantities obtained from DNS data. Based on this exercise either suitable models have been identified or new models have been proposed for the accurate closure of the unclosed terms of both \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) transport equations in the context of Reynolds Averaged Navier–Stokes (RANS) simulations.     

Statistical Analysis of Cross Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios <????> = 0.7 and <????> = 1.0 using three-dimensional simplified chemistry based Direct Numerical Simulations (DNS). The simulation parameters are chosen such that the combustion situation belongs to the thin reaction zones regime and a random bi-modal distribution of equivalence ratio ?? is introduced in the unburned gas ahead of the flame to account for mixture inhomogeneity. The DNS data has been used to analyse the statistical behaviour of the transport of the cross-scalar dissipation rate based on the fuel mass fraction Y _F and the mixture fraction ?? fluctuations $\,\tilde{\varepsilon}_{Y\xi}={\overline{\rho D\nabla Y_{F}^{\prime\prime}.\nabla \xi^{\prime\prime}} } \big/ {\bar {\rho }}$ (where $\bar{q}$ , $\tilde{q}={\overline{\rho q} } \big/ {\bar {\rho }}$ and $q^{\prime\prime} =q-\tilde {q}$ are Reynolds average, Favre mean and Favre fluctuation of a general quantity q) in the context of Reynolds Averaged Navier?CStokes simulations where ?? is the gas density and D is the gas diffusivity. The statistical behaviours of the unclosed terms in the $\tilde{\varepsilon }_{Y\xi } $ transport equation originating from turbulent transport T ₁, density variation T ₂, scalar?Cturbulence interaction T ₃, chemical reaction rate T ₄ and the molecular dissipation rate D ₂ have been analysed in detail. It has been observed that the contributions of T ₂, T ₃, T ₄ and D ₂ play important roles in the $\tilde{\varepsilon }_{Y\xi } $ transport for the globally stoichiometric cases, but in the globally fuel-lean cases the contributions of T ₂ and T ₄ become relatively weaker in comparison to the contributions of T ₃ and D ₂. The term T ₁ remains small in comparison to the leading order contributions of T ₃ and D ₂ for all cases, but the contribution of T ₁ plays a more important role in the low Damköhler combustion cases. The term T ₂ behaves as a sink term towards the unburned gas side but becomes a source term towards the burned gas side. The scalar?Cturbulence interaction term T ₃ has been found to be generally positive throughout the flame brush, but in globally stoichiometric cases the contribution of T ₃ becomes negative in regions of intense heat release. The combined contribution of (T ₄ ?C D ₂) remains mostly as a sink in all cases studied here. Models are proposed for the unclosed terms of the $\tilde{\varepsilon }_{Y\xi } $ transport equation in the context of Reynolds Averaged Navier?CStokes simulations, which are shown to satisfactorily predict the corresponding quantities extracted from the DNS data for all cases.     

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.

     

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.     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical Study and Lagrangian Modelling of Turbulent Heat Transport

A stochastic Lagrangian model for both fluid velocities and temperature fluctuations is evaluated from Direct Numerical Simulation of heat transport in homogeneous isotropic turbulence submitted to a linear mean temperature gradient. The first stage lies on the study of the Lagrangian fluid turbulence statistics (Lagrangian correlations functions) computed from predictions of DNS. They are crucial for the analysis and the modelling of the fluid turbulent properties along discrete particle trajectories. In the second stage, a velocity-scalar Lagrangian stochastic model is proposed and evaluated from the DNS data. The coefficients of the drift and diffusion terms of the model are determined by only Lagrangian timescales, temperature variance and turbulent flux. The shapes of correlation functions present a good agreement between DNS results and stochastic modelling approach.     

Effects of Lewis Number on Head on Quenching of Turbulent Premixed Flames: A Direct Numerical Simulation Analysis

The head on quenching of statistically planar turbulent premixed flames by an isothermal inert wall has been analysed using three-dimensional Direct Numerical Simulation (DNS) data for different values of global Lewis number Le(0.8, 1.0 and 1.2) and turbulent Reynolds number Re_t. The statistics of head on quenching have been analysed in terms of the wall Peclet number Pe (i.e. distance of the flame from the wall normalised by the Zel’dovich flame thickness) and the normalised wall heat flux Φ. It has been found that the maximum (minimum) value of Φ(Pe) for the turbulent Le=0.8 cases are greater (smaller) than the corresponding laminar value, whereas both Pe and Φ in turbulent cases remain comparable to the corresponding laminar values for Le=1.0 and 1.2. Detailed physical explanations are provided for the observed Le dependences of Pe and Φ. The existing closure of mean reaction rate \(\overline {\dot {\omega }}\) using the scalar dissipation rate (SDR) in the near wall region has been assessed based on a-priori analysis of DNS data and modifications to the existing closures of mean reaction rate and SDR have been suggested to account for the wall effects in such a manner that the modified closures perform well both near to and away from the wall.     

Direct Numerical Simulation of Head-On Quenching of Statistically Planar Turbulent Premixed Methane-Air Flames Using a Detailed Chemical Mechanism

A three-dimensional compressible Direct Numerical Simulation (DNS) analysis has been carried out for head-on quenching of a statistically planar stoichiometric methane-air flame by an isothermal inert wall. A multi-step chemical mechanism for methane-air combustion is used for the purpose of detailed chemistry DNS. For head-on quenching of stoichiometric methane-air flames, the mass fractions of major reactant species such as methane and oxygen tend to vanish at the wall during flame quenching. The absence of \(\text {OH}\) at the wall gives rise to accumulation of carbon monoxide during flame quenching because \(\text {CO}\) cannot be oxidised anymore. Furthermore, it has been found that low-temperature reactions give rise to accumulation of \(\text {HO}_{2}\) and \(\mathrm {H}_{2}\mathrm {O}_{2}\) at the wall during flame quenching. Moreover, these low temperature reactions are responsible for non-zero heat release rate at the wall during flame-wall interaction. In order to perform an in-depth comparison between simple and detailed chemistry DNS results, a corresponding simulation has been carried out for the same turbulence parameters for a representative single-step Arrhenius type irreversible chemical mechanism. In the corresponding simple chemistry simulation, heat release rate vanishes once the flame reaches a threshold distance from the wall. The distributions of reaction progress variable c and non-dimensional temperature T are found to be identical to each other away from the wall for the simple chemistry simulation but this equality does not hold during head-on quenching. The inequality between c (defined based on \(\text {CH}_{4}\) mass fraction) and T holds both away from and close to the wall for the detailed chemistry simulation but it becomes particularly prominent in the near-wall region. The temporal evolutions of wall heat flux and wall Peclet number (i.e. normalised wall-normal distance of \(T = 0.9\) isosurface) for both simple and detailed chemistry laminar and turbulent cases have been found to be qualitatively similar. However, small differences have been observed in the numerical values of the maximum normalised wall heat flux magnitude \(\left ({\Phi }_{\max } \right )_{\mathrm {L}}\) and the minimum Peclet number \((Pe_{\min })_{\mathrm {L}}\) obtained from simple and detailed chemistry based laminar head-on quenching calculations. Detailed explanations have been provided for the observed differences in behaviours of \(\left ({\Phi }_{\max }\right )_{\mathrm {L}}\) and \((Pe_{\min })_{\mathrm {L}}\). The usual Flame Surface Density (FSD) and scalar dissipation rate (SDR) based reaction rate closures do not adequately predict the mean reaction rate of reaction progress variable in the near-wall region for both simple and detailed chemistry simulations. It has been found that recently proposed FSD and SDR based reaction rate closures based on a-priori DNS analysis of simple chemistry data perform satisfactorily also for the detailed chemistry case both away from and close to the wall without any adjustment to the model parameters.     

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.     

Direct Numerical Simulation of Subsonic Round Turbulent Jet

Direct numerical simulation(DNS) of spatially developing round turbulent jet flow with Reynolds number 4,700 was carried out. Over 20 million grid points were used in this simulation. Fully compressible three-dimensional Navier–Stokes equations were solved. High order explicit spatial difference schemes and Runge–Kutta time integration scheme were used to calculate derivatives and time marching, respectively. Non-reflecting boundary conditions and exit zone techniques were adopted. Some refined computational grids were used in order to capture the smallest turbulent structures near the centerline of the jet. Low level disturbance were imposed on the jet inflow velocity to trigger the developing of turbulence. Turbulent statistics such as mean velocity, Reynolds stresses, third order velocity moments were obtained and compared with experimental data. One-dimensional velocity autospectra was also calculated. The inertial region where the spectra decays according to the k^− 5/3 was observed. The quantitative profiles of mean velocity and all of the third order velocity moments which were difficult to measure via experimental techniques were presented here in detail. The jet flow was proven to be close to fully self-similar around 19 jet diameters downstream of jet exit. The statistic data and revealed flow feature obtained in this paper can provide valuable reference for round turbulent jet research.     

Direct Numerical Simulation of Turbulent Compression Ramp Flow

A numerical procedure for the direct numerical simulation of compressible turbulent flow and shock–turbulence interaction is detailed and analyzed. An upwind-biased finite-difference scheme with a compact centered stencil is used to discretize the convective part of the Navier–Stokes equations. The scheme has a uniformly high approximation order and allows for a spectral-like wave resolution while dissipating nonresolved wave numbers. When hybridized with an essentially nonoscillatory scheme near discontinuities, the scheme becomes shock–capturing and its resolution properties are preserved. Diffusive parts are discretized with symmetric compact finite differences and an explicit Runge–Kutta scheme is used for time-advancement. The peculiarities of efficient upwinding and coupling procedures are described and validation results are given. Using direct numerical simulation data, some aspects of turbulent supersonic compression ramp flow are studied to demonstrate the effectiveness of the simulation procedure. Received 13 November 1997 and accepted 14 May 1998     

Large Eddy Simulation of Stratified and Sheared Flames of a Premixed Turbulent Stratified Flame Burner Using a Flamelet Model with Heat Loss

This paper presents large eddy simulations (LES) of the Darmstadt turbulent stratified flame burner (TSF) at different operating conditions including detailed heat loss modeling. The target cases are a non-reacting and two reacting cases. Both reacting cases are characterized by stratification, while one flame additionally features shear. In the regime diagram for premixed combustion, the studied flames are found at the border separating the thin reaction zones regime and the broken reaction zones regime. A coupled level set/progress variable model is utilized to describe the combustion process. To account for heat loss, an enthalpy defect approach is adopted and reformulated to include differential diffusion effects. A novel power-law rescaling methodology is proposed to integrate the enthalpy defect approach into the level set/progress variable model which is extensively validated in two validation scenarios. It is demonstrated that the LES with the newly developed model captures the influence of heat loss well and that the incorporation of heat loss effects improves the predictions of the TSF-burner over adiabatic simulations, while reproducing the experimentally observed flame lift-off from the pilot nozzle.     

Direct Numerical Simulations of Premixed Turbulent Flames with Reduced Chemistry: Validation and Flamelet Analysis

Direct numerical simulation is a very powerful tool to evaluate the validity of new models and theories for turbulent combustion. In this paper, direct numerical simulations of spherically expanding premixed turbulent flames in the thin reaction zone regime and in the broken reaction zone regime are performed. The flamelet-generated manifold method is used in order to deal with detailed reaction kinetics. The computational results are analyzed by using an extended flame stretch theory. It is investigated whether this theory is able to describe the influence of flame stretch and curvature on the local burning velocity of the flame. It is found that if the full profiles of flame stretch and curvature through the flame front are included in the theory, the local mass burning rate is well predicted. The influence of several approximations, which are used in other existing theories, is studied. When flame stretch is assumed constant through the flame front or when curvature of the flame front is neglected, the theory fails to predict the local mass burning rate. The influence of using a reduced chemistry model is investigated by comparing flamelet simulations with reduced and detailed chemistry.     

Large Eddy Simulations of the Darmstadt Turbulent Stratified Flames with REDIM Reduced Kinetics

Turbulent stratified combustion is often found in practical combustion devices, however, for large eddy simulations (LES) of it is still a challenge. In the present work, LES of the Darmstadt turbulent stratified flame (TSF) cases are conducted. In total, one isothermal flow case A-i2 and four reacting cases with different combinations of stratification and shear, i.e., A-r, C-r, E-r, G-r cases, are simulated. The employed sub-grid scale (SGS) combustion model is the REDIM-PFDF model, in which the chemical kinetics is reduced into a two-dimensional chemistry look-up table by the reaction-diffusion manifolds (REDIM) method, which performs a model reduction based on the coupling of the chemical kinetics with molecular transport. The fluctuation of scalars within the LES filter volume is modeled by the presumed filtered density function (PFDF). The overall good agreement of the statistics of velocity, temperature and species with the experimental data demonstrates the capability of the REDIM-PFDF model for TSF. Additionally, the probability distributions of the alignment angle, α, between the reaction layer and mixing layer, are analyzed in detail. It is shown that the probability distributions of the alignment angel vary with the axial distance from the jet nozzle. It also reveals that, with a stronger turbulence, the stratification effect can be weakened and the probability difference for finding ‘back-supported’ and ‘front-supported’ flame modes tends to decrease.     

各向同性湍流内颗粒碰撞率的直接模拟研究

对 Re_{\lambda } 约为51均匀各向同性湍流内 St_{k}(=\tau_{p}/\tau_{k}) 为 0 ~10.0 的 有限惯性颗粒的碰撞行为进行了直接数值模拟，以研究湍流对有限惯性 颗粒碰撞的影响. 结果表明，具有一定惯性颗粒的湍流碰撞率完全不同于零惯性的轻颗粒 (St_{k}=0) 和可忽略湍流作用的重颗粒 (St{k} \to \infty) , 其变化趋势极其复杂： 在Stk为 0~1.0 之间，颗粒的碰撞率随 St 的增加而近乎线性地剧烈增长，在 Stk≈1.0 3.0(对应的StE=τp/Te≈0.5)附近，颗粒碰撞率出现两个峰值，在Stk>3.0以后，颗粒的碰撞率随惯性增 大而逐渐趋向于重颗粒极限；在峰值处，有限惯性颗粒的平均碰撞率的峰值较轻颗粒增强了 30倍左右. 为进一步分析湍流作用下颗粒碰撞率的影响因素，分别使用可能发生碰撞 的颗粒对的径向分布函数和径向相对速度来量化颗粒的局部富集效应和湍流掺混效应，表明 St_{k} \approx 1.0 时局部富集效应最为强烈，使得颗粒的碰撞率出现第1个峰值； 湍流掺混效应则随着颗粒Stk的增大而渐近增大；局部富集和湍流掺混联合作用的结果， 使得颗粒碰撞率在 St_{k} \approx 3.0 附近出现另一个峰值.     

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.     

Effects of Body Forces on Turbulent Kinetic Energy Transport in Premixed Flames

Flow, Turbulence and Combustion - The effects of buoyancy on turbulent premixed flames are expected to be strong due to the large changes in density between the unburned and fully burned gases. The...