LES-CMC of a Partially Premixed,Turbulent Dimethyl Ether Jet Diffusion Flame

Large-eddy simulations have been coupled with a conservative formulation of the conditional moment closure (CMC) approach for the computation of a turbulent, partially-premixed dimethyl-ether jet flame. Two different numerical setups and 3 different detailed chemical mechanisms were investigated. The results are compared with measurements of velocity, temperature, and major and intermediate species. The general agreement between simulations and experiments is very good, and differences between the different mechanisms are limited to the predicted concentrations of intermediates only. Larger differences can be observed if the CMC grid size is reduced. This is due to reduced averaging effects on the conditionally averaged dissipation rates that allow to better capture high dissipation events that lead to larger deviations from a fully burning solution. A high CMC resolution provides excellent agreement with experiments throughout the flame and the results demonstrate CMC’s capability to accurately predict turbulence-chemistry interactions in partially-premixed flames involving complex chemistry.     

Numerical Study of Ignition Process in Turbulent Shear-less Methane-air Mixing Layer

In this work, ignition process in a turbulent shear-less methane-air mixing layer is numerically investigated. A compressible large eddy simulation method with Smagorinsky sub-grid scale model is used to solve the flow field. Also, a thickened flame combustion model and DRM-19 reduced mechanism are used to compute species distribution and the heat release. Non-reacting mean and RMS axial velocity profiles and mean mixture fraction are validated against experimental data. Instantaneous mixture fraction contours show that the large bursts penetrate from the fuel stream into that of the oxidizer and vice versa and a random behaviour in the cross-stream direction. Flame kernel initiation, growth and propagation are analysed and compared with the experimental data. The ignition results show that the flame is not stable and blow-off occurs, but a more detailed investigation shows that local and short time flame stabilization exist during blow-off. During these local stabilization, heat release increased at the upstream edge of the flame. Most_upstream flame edge scalar analysis shows that the methane mass fraction has a dominant role in the local flame stabilization. OH, HO₂, CH₂O and heat release contours demonstration reveal that HO₂ and CH₂O mass fraction as well as the heat release reach a maximum on the border of the flame, but the maximum OH concentration is located in the middle of flame kernel.     

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.     

半封闭狭缝湍流冲击射流的数值模拟

将Yakhot和Orszag提出的RNGk-ε模型推广应用于半封闭狭缝冲击射流场的数值模拟，以评价该模型对这种复杂湍流的预测能力。将计算得到的流场平均速度分布、湍流强度分布和流函数分布与标准k—ε模型的预测结果以及相应的实验数据进行了比较，结果表明：RNGk—ε模型的预测结果总体上要好于标准k—ε模型，但与实验值相比，所有预测结果都还存在不同程度的误差，尤其是近壁区和滞止点较远下游处的湍流强度分布。说明RNG模型虽然已在某些湍流的预测中取得了一定的成功，但要定量准确地预测冲击射流场，还必须针对其流动特征对模型加以改进。     

Large-eddy Simulation of Triangular-stabilized Lean Premixed Turbulent Flames: Quality and Error Assessment

In this numerical study, an algebraic flame surface wrinkling (AFSW) reaction submodel based on the progress variable approach is implemented in the large-eddy simulation (LES) context and validated against the triangular stabilized bluff body flame configuration measurements i.e. in VOLVO test rig. The quantitative predictability of the AFSW model is analyzed in comparison with another well validated turbulent flame speed closure (TFC) combustion model in order to help assess the behaviour of the present model and to further help improve the understanding of the flow and flame dynamics. Characterization of non-reacting (or cold) and reacting flows are performed using various subgrid scale models for consistent grid size variation with 300,000 (coarse), 1.2 million (intermediate) and 2.4 million (fine) grid cells. For non-reacting flows at inlet velocity of 17?m/s and inlet temperature 288?K, coarse grid leads to over prediction of turbulence quantities due to low dissipation at the early stage of flow development behind the bluff body that convects downstream eventually polluting the resulting solution. The simulated results with the intermediate (and fine) grid for mean flow and turbulence quantities, and the vortex shedding frequency (f_s) closely match experimental data. For combusting flows for lean propane/air mixtures at 35?m/s and 600?K, the vortex shedding frequency increase threefold compared with cold scenario. The predicted results of mean, rms velocities and reaction progress variable are generally in good agreement with experimental data. For the coarse grid the combustion predictions show a shorter recirculation region due to higher turbulent burning rate. Finally, both cold and reacting LES data are analyzed for uncertainty in the solution using two quality assessment techniques: two-grid estimator by Celik, and model and grid variation by Klein. For both approaches, the resolved turbulent kinetic energy is used to estimate the grid quality and error assessment. The quality assessment reveals that the cold flows are well resolved even on the intermediate mesh, while for the reacting flows even the fine mesh is locally not sufficient in the flamelet region. The Klein approach estimates that depending on the recirculation region in cold scenario both numerical and model errors rise near the bluff-body region, while in combusting flows these errors are significant behind the stabilizing point due to preheating of unburned mixture and reaction heat release. The total error mainly depends on the numerical error and the influence of model error is low for this configuration.     

A Novel One- and Zero-Dimensional Model for Turbulent Jet Ignition

Flow, Turbulence and Combustion - Turbulent jet ignition (TJI) is a promising combustion technology for burning highly diluted air-fuel mixtures. Computationally efficient models to assess the...     

Evaluations of SGS Combustion,Scalar Flux and Stress Models in a Turbulent Jet Premixed Flame

A newly developed fractal dynamic SGS (FDSGS) combustion model and a scale self-recognition mixed (SSRM) SGS stress model are evaluated along with other SGS combustion, scalar flux and stress models in a priori and a posteriori manners using DNS data of a hydrogen-air turbulent plane jet premixed flame. A posteriori tests reveal that the LES using the FDSGS combustion model can predict the combustion field well in terms of mean temperature distributions and peak positions in the transverse distributions of filtered reaction progress variable fluctuations. A priori and a posteriori tests of the scalar flux models show that a model proposed by Clark et al. accurately predicts the counter-gradient transport as well as the gradient diffusion, and introduction of the model of Clark et al. into the LES yields slightly better predictions of the filtered progress variable fluctuations than that of a gradient diffusion model. Evaluations of the stress models reveal that the LES with the SSRM model predicts the velocity fluctuations well compared to that with the Smagorinsky model.     

微型射流涡流器对湍流边界层控制的数值研究

运用数值方法,模拟出展向分布的同向倾斜微型射流列与平板湍流边界层相互作用形成流向涡列的流场结构,验证了利用其来对湍流边界层进行控制的可能性.随射流间距减小,流向涡列控制作用流向渗透能力增强,但作用区域减小;随射流速度提高,流向涡列控制作用增强,但过大的射流速度反而会导致流向涡列在局部区域内控制作用的下降;随射流俯仰角减小、倾斜角增大,流向涡列初始控制作用增强,但过小的俯仰角、过大的倾斜角会导致流向涡列流向控制区域明显缩小.要保证流向涡列具有较强的湍流边界层控制作用,必须通过合理配置射流列各主要参数,在保证各流向涡具有一定强度的同时,还要确保各流向涡在形成时部分嵌入边界层内部.     

Large-Eddy Simulation of Lean Premixed Turbulent Flames of Three Different Combustion Configurations using a Novel Reaction Closure

This large eddy simulation (LES) study is applied to three different premixed turbulent flames under lean conditions at atmospheric pressure. The hierarchy of complexity of these flames in ascending order are a simple Bunsen-like burner, a sudden-expansion dump combustor, and a typical swirl-stabilized gas turbine burner–combustor. The purpose of this paper is to examine numerically whether the chosen combination of the Smagorinsky turbulence model for sgs fluxes and a novel turbulent premixed reaction closure is applicable over all the three combustion configurations with varied degree of flow and turbulence. A quality assessment method for the LES calculations is applied. The cold flow data obtained with the Smagorinsky closure on the dump combustor are in close proximity with the experiments. It moderately predicts the vortex breakdown and bubble shape, which control the flame position on the double-cone burner. Here, the jet break-up at the root of the burner is premature and differs with the experiments by as much as half the burner exit diameter, attributing the discrepancy to poor grid resolution. With the first two combustion configurations, the applied subgrid reaction model is in good correspondence with the experiments. For the third case, a complex swirl-stabilized burner–combustor configuration, although the flow field inside the burner is only modestly numerically explored, the level of flame stabilization at the junction of the burner–combustor has been rather well captured. Furthermore, the critical flame drift from the combustor into the burner was possible to capture in the LES context (which was not possible with the RANS plus k–ɛ model), however, requiring tuning of a prefactor in the reaction closure.     

Effects of the Local Flow Topologies Upon the Structure of a Premixed Methane-air Turbulent Jet Flame

Local flow topologies have been identified and their interactions with the iso-scalar surfaces geometries have been investigated using the results of a three-dimensional direct numerical simulation (DNS) of a turbulent premixed methane-air flame in a piloted Bunsen burner configuration with tabulated chemistry. The universal teardrop shape of the joint probability density function (jpdf) of the second and third invariants of the velocity-gradient tensor disappears in the different flame regions under study. A ‘canonical’ vortex, which affects the fine-scale structure of the turbulent premixed flame, has been identified and analyzed at three times, differing by increments of the order of the Kolmogorov time micro-scale.     

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.     

Experimental and Numerical Study of the Scalar Turbulent/Non-Turbulent Interface Layer in a Jet Flow

Based on two large-eddy simulations (LES) of a non-reacting turbulent round jet with a nozzle based Reynolds number of 8,610 with the same configuration as the one that has recently been investigated experimentally (Gampert et al., 2012; J Fluid Mech, 2012; J Fluid Mech 724:337, 2013), we examine the scalar turbulent/non-turbulent (T/NT) interface layer in the mixture fraction field of the jet flow between ten and thirty nozzle diameters downstream. To this end, the LES—one with a coarse grid and one with a fine grid—are in a first step validated against the experimental data using the axial decay of the mean velocity and the mean mixture fraction as well as based on radial self-similar profiles of mean and root mean square values of these two quantities. Then, probability density functions (pdf) of the mixture fraction at various axial and radial positions are compared and the quality of the LES is discussed. In general, the LES results are consistent with the experimental data. However, in the flow region where the imprint of the T/NT interface layer is dominant in the mixture fraction pdf, discrepancies are observed. In a next step, statistics of the T/NT interface layer are studied, where a satisfactory agreement for the pdf of the location of the interface layer from the higher resolved LES with the experimental data is observed, while the one with the coarse grid exhibits considerable deviations. Finally, the mixture fraction profile across the interface is investigated where the same trend as for the pdf of the location is present. In particular, it is found that the sharp interface that is present in experimental studies (Gampert et al., J Fluid Mech, 2013; Westerweel et al., J Fluid Mech 631:199, 2009) is less distinct in the LES results and rather diffused in radial direction outside of the T/NT interface layer.     

A Numerical Study Promoting Algebraic Models for the Lewis Number Effect in Atmospheric Turbulent Premixed Bunsen Flames

This numerical investigation carried out on turbulent lean premixed flames accounts for two algebraic – the Lindstedt–Vaos (LV) and the classic Bray–Moss–Libby (BML) – reaction rate models. Computed data from these two models is compared with the experimental data of Kobayashi et al. on 40 different methane, ethylene and propane Bunsen flames at 1 bar, where the mean flame cone angle is used for comparison. Both models gave reasonable qualitative trend for the whole set of data, in overall. In order to characterize quantitatively, firstly, corrections are made by tuning the model parameters fitting to the experimental methane–air (of Le = 1.0) flame data. In case of the LV model, results obtained by adjusting the pre-constant, i.e., reaction rate parameter, C_R, from its original value 2.6 to 4.0, has proven to be in good agreement with the experiments. Similarly, for the BML model, with the tuning of the exponent n, in the wrinkling length scale, L_y = C_l⋅ l_x(s_L/u′)ⁿ from value unity to 1.2, the outcome is in accordance with the measured data. The deviation between the measured and calculated data sharply rises from methane to propane, i.e., with increasing Lewis number. It is deduced from the trends that the effect of Lewis number (for ethylene–air mixtures of Le = 1.2 and propane–air mixtures of Le = 1.62) is missing in both the models. The Lewis number of the fuel–air mixture is related to the laminar flame instabilities. Second, in order to quantify for its influence, the Lewis number effect is induced into both the models. It is found that by setting global reaction rate inversely proportional to the Lewis number in both the cases leads to a much better numerical prediction to this set of experimental flame data. Thus, by imparting an important phenomenon (the Lewis number effect) into the reaction rates, the generality of the two models is enhanced. However, functionality of the two models differs in predicting flame brush thickness, giving scope for further analysis.     

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.     

Large Eddy Simulation of a Propagating Turbulent Premixed Flame

A large eddy simulation of a turbulent premixed flame propagatingthrough a chamber containing a square obstruction is presented anddiscussed. The governing equations for compressible, reacting flowsare Favre filtered and turbulence closure is achieved using thedynamic Smagorinsky subgrid model. A simple flame surface densitymodel based on the flamelet concept is employed for the subgrid scalereaction rate. The simulation gives very good agreement with experimentalresults for the speed and the shape of the flame as it propagates throughthe chamber. The peak pressures, however, are underpredicted by20–30%. Reasons for this are discussed and it is concluded that amore sophisticated combustion model is required for complex flowssuch as this one, if a more accurate prediction of the pressureis to be achieved.     

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.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.     

Turbulent Structure of a Gas-Liquid Impinging Jet

The turbulent structure of a submerged axisymmetric impinging jet containing small gas bubbles is studied experimentally under conditions of periodic external excitation. On the basis of measuring the surface-friction pulsatory component in the jet impinging on an obstacle, the effect of the suppression of large-scale eddies at large gas volume fractions is registered. The conditions of resonant growth of coherent structures and the suppression of wide-band turbulence are determined for both the single-phase and the two-phase impinging jet. An analysis of the development of different pulsatory friction components in the impinging-jet gradient region is presented.     

Geometrical Properties and Turbulent Flame Speed Measurements in Stationary Premixed V-flames Using Direct Numerical Simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier?CStokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u??/s _L , of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17?? _th downstream of the flame holder, where ?? _th is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, s _T is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable $\tilde{c}$ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, s _C , calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518?C535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of $\tilde{c}$ is taken into account.